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Abstract 



We have observed extragalactic sources Mrk421 (2;=0.30) and H1426+428 (z=0.129) in very 
high energy 7-ray region using TACTIC telescope at Mount Abu (24.6°N, 72.7°E, 1300 m 
asl), Rajasthan, India. Both of these sources belong to the blazar sub-class of active galactic 
nuclei whose jets are closely aligned to our line of sight. Blazars are most extreme and pow- 
erful variable sources of radiation of energies ranging from the radio to the 7-ray regimes. 
Mrk421 was observed during 07 December 2005 to 18 April 2007 for a total of ~331 hours 
while H1426+428 was observed for ~244 hours between 22 March 2004 and 1 1 June 2007. 
Detailed analysis of Mrk421 data revealed the presence of a 7-ray signal with the excess 
of ~(951 ± 82) 7-ray Uke events corresponding to a statistical significance of ~ 12.0(7 in 
r^97 hours during 27 December 2005 to 07 February 2006. The combined average flux 
was measured to be ~(1.04 ± 0.14) Crab units. Obseved time averaged differential en- 
ergy spectrum was found to be consistent with both a pure power law d^/dE — (4.66 ± 
0.46) X 10~^^ ^(-3.ii±o.ii) s~^ TeV~^ and power law with an exponential cutoff func- 
tions d^/dE = foE-'^exp(-E/Eo) with /o = (4.88 ± 0.38) x IQ-^^ cm-'^ s'^ TeV-\ 
r = 2.51 ± 0.26 and Eq = (4.7 ± 2.1)^6^. The later function fits the observed spectrum 
well with lower value of reduced x^- However, during the observation period from 18 De- 
cember 2006 to 18 April 2007, the source was found to be in a low emission state and we 
have placed an upper limit of 1.44 x 10'^"^ photons cm~'^ s~^ at 3a level on the integrated 
7-ray flux above 1 TeV. In addition, the TACTIC light curves have also been compared 
with those obtained with the RXTE/ASM and Swift/BAI satellite -based experiments. Fur- 
ther, we have obtained the SED of Mrk421 in the VHE range using TACTIC telescope and 
compared it with those obtained by the HESS, MAGIC, VERITAS, Whipple, HEGRA and 
CAT telescopes. The detection of TeV photons from the distant object H1426-I-428 have 
important implications for estimating the EEL, in addition to the understanding of the par- 
ticle acceleration and 7-ray production mechanisms in the AGN jets. Detailed analysis of 
H1426-I-428 data does not indicate the presence of any statistically significant TeV 7-ray sig- 
nal. Accordingly we have placed an upper limit of < 1.18 x 10~^'^ photons cm~'^ s~^ on the 
integrated 7-ray flux above 1 TeV at 3a significance level. 
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Very high energy gamma-ray 
astrophysics 

1.1 Introduction 

Cosmic 7-rays represent the highest energy band of the electromagnetic spectrum and pro- 
vide insight into some of the most dynamic processes in the Universe. The window to the 
high energy universe was opened a few decades ago with the advent of X-ray astronomy and 
recent advances have given us access to observe photons of still higher energy. Although the 
distinction between X-rays and 7-rays is not very sharp, often X-rays are considered to be 
photons produced by atomic, thermal or non-thermal processes, while 7-rays are produced 
in nuclear or non-thermal processes. Above IMeV all photons are 7-rays but this very 
wide energy range needs some subdivision for clarity. Weekes [I294II defined "High Energy" 
(HE) gamma rays as those between 30 MeV to 10 GeV, "Very High Energy" (VHE) as pho- 
ton energies from 10 GeV to 100 TeV (10^° to 10^^ eV), "Ultra High Energy" (UHE) from 
100 TeV to 100 PeV (10^^ to 10^^ eV) and "Extremely High Energy" (EHE) from 100 PeV 
to lOOi^eV^ (10^^ to 10^° eV). These energy ranges overlap the regions of the spectrum 
historically covered by the different detection techniques: satellite based detectors (HE), at- 
mospheric Cherenkov detectors (VHE), air shower arrays (UHE) and nitrogen fluorescence 
detectors or very large air shower arrays (EHE). In this thesis I shall follow the convention of 
Very High Energy (VHE) 7-ray astronomy as the study of astrophysical sources of photons 
> WGeV. 

In this chapter, we discuss the scientific motivation for pursuing 7-ray astronomy, produc- 
tion and absorption processes of VHE 7-rays along with their propagation in the interstellar 
medium. We also discuss various types of galactic and extragalactic sources of 7-rays . 
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1.2 Brief history 

While 7-ray astronomy in the 100 MeK energy region started with the seminal paper of 
Morrison [I234L astronomy at TeV energies seems to date from the prescient prediction by 
Cocconi HI 2011 in 1959 of the detectability of TeV 7-rays from the Crab Nebula. Although 
the model for 7-ray emission adopted by him was not correct (as it over estimated the de- 
tected flux by orders of magnitude) he sowed the seeds for the first serious atmospheric 
Cherenkov experiments to detect VHE 7-rays from cosmic sources. Cocconi's prediction 
was based on the detection of air showers from 7-rays by existing particle arrays. The detec- 
tion of Atmospheric Cherenkov Radiation (ACR) from Cosmic Ray (CR) air showers H 15911 
opened the path for ground-based studies of 7-rays . The first systematic VHE 7-ray ob- 
servations using the Atmospheric Cherenkov Technique (ACT), were made in Crimea by a 
group at the Lebedev Institute in 1960 [I119II . Operating at an energy threshold of ~ ATeV, 
no statistically significant fluxes were detected from any source directions by this telescope. 
However, they laid the ground work for the future development of the ACT. This was fol- 
lowed by the development of Whipple 7-ray telescope in 1968, with a large optical reflector 
(10 m diameter), at Mount Hopkins in Southern Arizona. While larger collection area of the 
Whipple telescope led to the reduction in energy threshold, it did not immediately result in 
a significant improvement in the flux sensitivity. However, observations with the Whipple 
telescope resulted in the first detection of 7-rays from the Crab Nebula [I148L More details 
of the early history of 7-ray astronomy can be found in a review paper by Weekes [I294II . 

In the first couple of decades of observations with marginal detections, it was clear that 
the sensitivity of ACT would need to be substantially improved by an efficient means of re- 
jecting the CR background. Now it is known that even for the strongest 7-ray sources, the CR 
background events outnumber the 7-rays by a factor of about 1000 or so. With the advances 
in computers, rigorous Monte Carlo simulations of shower development revealed some char- 
acteristic differences in the electromagnetic (initiated by 7-rays ) and the hadronic showers. 
In order to get the best method for background rejection, a number of discriminants were 
attempted e.g. differences in the spectral contents exploited by the early Whipple and the 
Crimean groups, lateral distribution differences utilized by CELESTE lfT35ll . STACEE ll298]l 
and Pachmarhi group [|94l . It was the idea originally recognised by Jelley and Porter II191II 
that dramatically enhanced the performance of atmospheric Cherenkov telescopes. The im- 
age of ACR incident on the telescope could be used to reject the CR background by signif- 
icantly improving the angular resolution of the telescope and the shape of image could be 
used as an independent discriminant [I182II . Demonstrating the power of the imaging atmo- 
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spheric Cherenkov technique, the Whipple group firmly established the detection of TeV 
7-ray emission from the Crab Nebula with a statistical significance of 9.0a in 1989 [1296| . 

The stereoscopic approach (detection of ACR from the shower in different projections) 
allows precise determination of the arrival direction of primary 7-rays on an event by event 
basis and effective suppression of the background light due to CR showers and from local 
muons, therefore provides further improvement in the performance of the imaging atmo- 
spheric Cherenkov technique. The observations of the Crab Nebula with HEGRA [I131II con- 
firmed the predictions about the stereoscopic Imaging Atmospheric Cherenkov Telescope 
(lACT) arrays [35 1. The collective experience of all the historical experiments pointed to- 
wards a larger collection area of the lACT, small pixel size of the camera and the power of 
stereoscopic technique as the parameters with the potential for improving the sensitivity of 
the lACT 

1.3 Cosmic rays 

Ever since their discovery by the Austrian physicist Victor Hess, nearly a century ago in 
1912, CRs have been puzzling physicists and astrophysicists. He observed that an ionizing 
radiation impinges at the top of the Earth's atmosphere continuously. The origin of CRs 
(exact sites and processes of their production), their composition over the entire observed 
energy range and their propagation in space basically still remains unknown after so many 
years of research. The CR spectrum measured by different experiments from a few GeV 
energy up to 10^° eV is shown in Fig. |l.l[ Up to energies < 10^^ eV CRs can directly be 
measured by satellite based experiments, however at higher energies their small flux makes 
the direct detection by these experiments difficult. As shown in Fig. |l.l| the CR flux decreases 
from ~ 1 particle at energies around 100 GeV to below 1 particle year~^ for 

the energy 10^^ eV. For energies above 10^^ eV it is possible to detect the CR initiated at- 
mospheric cascade of secondary particles called Extensive Air Shower (EAS), with ground 
based detectors at high altitudes. The detection of charged secondary particles and the mea- 
surement of ACR, produced by an air shower are the two established techniques to observe 
EAS. The secondary particles produced by lower energy primaries do not survive till the 
ground level while ACR produced by them can still be measured. The energy spectrum in 
the region beyond 1 GeV is unaffected by the magnetic field of the Earth and the propagation 
of CR charged particles to the Earth through the solar wind. The differential flux dN/ dE, in 
this energy region can be described by a power law, dN/dE oc E~°'. Upto the knee point 
located at 5 x 10^^ eV the spectral index has a value of a = 2.7 while the spectrum steep- 
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Fig .1.1: The all-particle spectrum of CRs measured by different experiments. This figure has been 
adapted from S. S wordy [61. 



ens from the knee and the index assumes a value a = 3.3 upto a point called the ankle (at 
~ 3 X 10^® eV). From the ankle onwards and upto ~ 3 x 10^^ eV the differential spectral 
index has a value of a = 2.6. More recent experiments HiRes [[88l and Auger ETll have 
shown the suppression of CR spectrum by a factor of two in comparison to the power law 
extrapolation above 3 x 10^^ eV. 

The suppression at the highest energy end in CR spectrum is presumably due to the 
Greizen-Zatsepin-Kuz'min (GZK) cutoff which is a theoretical upper limit on the energy of 
CRs from distant sources. Independently computed in 1966 by Greisen [I170II . and Kuz'min 
and Zatsepin B03il . this limit is caused due to CR energy losses in their interaction with 
Cosmic Microwave Background Radiation(CMBR). They predicted that CRs with energies 
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over the threshold energy of 5 x 10^ eV would interact with CMBR photons to produce 
pions via the resonance 

ICMBR + p A* p + 7r° 

ICMBR + p A* ^ n + 7r+ 

This process continues until the CR energy falls below the pion production threshold. Due 
to the mean path associated with this interaction, extragalactic CRs traveling over distances 
larger than 50 Mpc and with energies greater than this threshold should never be observed 
on Earth. This distance is also known as GZK horizon. The existence of GZK cutoff has 
been a much debated subject. A number of observations have been made by the AGASA 
experiment that appeared to show CRs from distant sources with energies above this limit 
(called ultra-high-energy CRs, or UHECRs). The observed existence of these particles was 
the so-called GZK paradox or cosmic ray paradox. 

The chemical composition of CR has been measured directly up to the knee point and 
found to comprise mainly charged particles, about 98% of these particles are protons and 
nuclei, whilst 2% are electrons. Of the protons and nuclei, about 87% are protons, 12% are 
helium nuclei and the remaining 1% are heavier nuclei [I218L The spectrum of electrons ( 
and positrons) is steeper than the one of protons and nuclei. The 7-ray component is just a 
tiny fraction < 10^''. It is generally believed that CRs with the energy below the knee value 
mostly originate from galactic accelerators and they have been confined inside our Galaxy 
by its magnetic field which is of the order of 1 jiG. The region between knee and ankle is 
still under debate, while at higher energies above 10^^ eV particles are believed to be mostly 
of extragalactic origin, since the galactic magnetic field is unable to trap them within the 
galaxjQ More recently the Auger collaboration has announced that CRs of energy above 
10^^ eV appear to come from Active Galactic Nuclei (AGN) located within 75 Mpc of our 
Galaxy [I146II . The value of extragalactic magnetic field is of the order of 5 x 10^^*^ G within 
a radius of 75 Mpc from the solar system [|133|| . 

The energy spectrum of CR shows very few variations which can be interpreted as the 
evidence that the same mechanism and the same class of sources are responsible for most 
of their production. The best known possible sites of production of CR includes Supernova 
Remnants (SNRs of shell type), Plerions, Pulsars, Microquasars (X-ray binary systems) and 
Young stars ( and young open star clusters) within our galaxy while extragalactic ones are 
AGN, Gamma Ray Bursts (GRBs), Starburst galaxies and Clusters of galaxies. 

'Rigidity R = B x rg where B is Magnetic field and Vg is the gyroradius 
gyroradius Vg = niv±^ / qB 
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Despite being the most abundant component of CRs, protons do not retain the directional 
information about their origin as they get deflected by the magnetic fields. For a proton of 
energy E crossing a magnetic field B, the deflection radius is approximately given by [I132II : 

, E/lTeV 
R/lpc . (0.01) 

This means that a proton of energy below 300 PeV|^ coming from the galactic center (~ 
8 kpc from the Earth) will lose its directional information. 7-rays , neutrinos and neutrons 
are possible tracers for search and investigation of CR sources because they are not deflected 
by magnetic fields and thus point to the place of their creation. Mean decay time of a neu- 
tron being 886 s, implies that only the highest energy neutrons can be detected over large 
distance^ The detection of neutrinos from AGN allows to distinguish between the two ma- 
jor competing classes of models of 7-ray production mechanisms as they favor the hadronic 
one (discussed in next chapter) and therefore neutrino observations are complementary to 7- 
ray astronomy. However detection of neutrinos is very difficult as they have small interaction 
cross section, necessitating very large size detectors (> 1 km^). Currently 7-ray astronomy 
is the most suitable and efficient probe to study CR sources. 



1.4 Scientific motivation 

It is known that VHE astrophysical processes produce relativistic particles and the associated 
gamma radiation over a wide range of energy. 7-rays provide one of the best windows into 
the non-themal universe and a means of probing fundamental physics beyond the reach of 
terrestrial accelerators. The energy region of 7-rays has been explored so far rather inhomo- 
geneously. Neverthless, unprecedented progress has been made during the last decade in this 
field which has the following scientific motivations: 

1.4.1 Origin of cosmic rays 

As the cosmic 7-rays reach us without changing their original direction, they provide in- 
formation about their sources. CRs are remarkable for a variety of reasons. Some of their 
important characteristics are: 

^PeV = IQis eV 

^Neutrons from the center of the Galaxy can only be detected if it has energy > 10^* eV 
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• The energy range of CRs is enormous i.e. 10® eV to 10^° eV. The presence of VHE 
CRs indicates the evidence of powerful astrophysical accelerators. 

• The sources which power CRs have a total luminosity of 10^^ erg s^^. There is ev- 
idence that CRs have existed for as long as the galaxy itself (~ 10^° yr) but the es- 
timated galactic containment time is only ~ 10^|/r 016011 . Hence CRs must continu- 
ously be replenished during the lifetime of the galaxy. 

• CRs are abundant and play an important role in the energy balance of the galaxy. They 
have an energy density of ~ leV cm~^ comparable to that contained in the galactic 
magnetic field or in the CMBR. 

• CRs are the largest source of material reaching the Earth from outside the Solar Sys- 
tem. Their chemical composition (modulated by propagation effects) reflects the nu- 
cleosynthetic processes occurring at their origin. 

1.4.2 Astrophysical considerations 

While the typical energy of a visible photon is ~ 1 eV, a VHE 7-ray has an energy more 
than 10 orders of magnitude greater than this. Based on this energy difference the thermal 
radiation as the underlying production mechanism for high energy photons is not applicable 
for 7-rays. It is expected that 7-ray s are produced via particle interactions in very powerful 
accelerators. 7-ray astronomy therefore probes astrophysics not accessible by conventional 
astronomical techniques. 

As discussed later in this chapter, 7-rays can also be used as probes of Extragalactic 
Background Light (EBL), especially VHE photons which interact with optical to infrared 
background. The measured VHE spectra are used to set upper limits on EBL density. The 
TeV catalog of extragalactic sources contains a large fraction of nearby sources, which is 
possibly due to the absorption of 7-rays by the EBL. It is expected that the number of ex- 
tragalactic sources will increase significantly if the analysis threshold energy of a VHE 7- 
ray telescopes is reduced significantly below 100 GeV, as at these energies photons penetrate 
deeper into space. 
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1.5 Gamma rays 

1.5.1 Production mechanisms of 7-rays 

• Bremsstrahlung 

When a charged particle is deflected in the electric field of a nucleus it emits elec- 
tromagnetic radiation called Bremsstrahlung. Classically, the acceleration produced 
by a nucleus (Ze) on a particle of charge e and mass m is proportional to Ze^/m. 
The complete quantum mechanical treatment of electron bremsstrahlung by an atom 
is complex because of the effects of screening by the atomic electrons and the finite 
nuclear radius II147II . This process is particularly important for relativistic electrons in 
the presence of atomic and molecular material wherein electrons will be deflected and 
emit 7-rays via the bremsstrahlung [117511 mechanism. In astrophysical situations, the 
production of 7-rays depends on the electron energy distribution and the gas density. 
The spectrum of emitted photons has a power law form with the same spectral index as 
that of the parent accelerated particles. To produce a 7-ray photon of a given energy, 
the electron should have more than twice that energy [I257II . Bremsstrahlung in the 
field of atomic electrons is also important when the numbers of e^ are much higher 
than those of the nuclei as in regions where temperatures are more than 10^° K. 

• Synchrotron radiation 

Synchrotron radiation is the process in which radiation is emitted by relativistic charged 
particles in the presence of a magnetic field. Li a uniform magnetic field a non- 
relativistic electron moves in a spiral path at a constant pitch angle 9. Its velocity 
along the field lines is constant while it gyrates about the magnetic field direction at 
the angular frequency of Larmour precession 

eB 

col = (1.1) 

me c 

where mg is rest mass of electron and B is magnetic field intensity normal to the 
velocity vector of the electron. At relativistic energies, radiation is beamed into a 
cone of angle 9 = nieC^/E [I 295II and occurs as a continuum spectrum. The critical 
frequency at which the maximum power is emitted is 



3 eB 




(1.2) 



2 m-e c 
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where 7 is Lorentz factor of electron. The energy loss rate is given by 

dE 1 dE , 2e^ 



dx c dt 3 c* 



17^5^ ergs cm ^ (1.3) 



where E is in ergs and B in Gauss^T^. 

Because of higher mass of proton {mp/rrie — 1836) the energy loss rate appears 
{mp/rrieY — 10^^ times slower, characteristic frequency of synchrotron radiation is 
{mp/rrie)^ ^ 6 x 10^ times smaller and time needed for radiation of synchrotron 7-ray 
photon is (mp/meY^^ — 1.5 x 10^ times longer as compared to an electron of the 
same energy and therefore proton synchrotron radiation is considered as an inefficient 
process [|37l . 

If the radiating electrons have a power law energy spectrum of the form 

N,{E,)=KEf (1.4) 

the synchrotron photon number spectrum will be of the form 

G,{E,) oc (1.5) 

which is flatter than that of the parent electron spectrum for values of F > 1 02571 . 
It should be noted here that in other windows of electromagnetic radiation, spectra 
are generally expressed in terms of the power radiated. In order to find the radiated 
power spectrum(iym^ifz~^) one has to multiply the equation (1.5) by E^ leading to 
synchrotron power spectrum II257B 

J,[E,) oc (1.6) 

Synchrotron radiation of accelerated electrons is one of the most important processes 
in the non-thermal universe and in the context of VHE 7-ray s, this is the usual process 
for generating the seed photons for Inverse Compton (IC) scattering (discussed below). 
However, UHE CRs (electrons and protons) can directly emit synchrotron radiation at 
VHE range. Quantitatively, to produce 1.0 TeV 7-ray in the presence of a magnetic 
field of 1 G, the electron should have an energy of ~ 7.2 x 10^^ eV . 

Inverse Compton (IC) scattering 

In the IC scattering a relativistic electron or positron of energy E^ scatters off a low 
energy photon of energy e and transfers part of its energy to the photon. The energy of 



Chapter 1 . Very high energy gamma-ray astrophysics 



10 



Compton boosted photon is given by [I257II 



and 



i?^ ~ e 7^ for e <^ rrieC^ 



~ Ee for 7 e me 



(1.7) 



(1.8) 



The corresponding cross sections are given by ^25711 

27e 



CTt 1 



nip 



and 



nip c 



7e 



In 



27e 



nip 



1/2 



(1.9) 



(1.10) 



where ax = Snrl/S is the Thomson cross section and has a value of ~ 6.65 x 10^^^ 
cm^. In order to produce a 1.0 TeV photon against microwave background which has 
mean energy of 7 x 10^'* eV requires an electron of energy ~ 1.7 x 10^^ eV, while 
against starlight (e ~ 1.0 eV) and X-rays (e ~ 10 keV) nearly 1.0 TeV electron can 
produce the photon of similar energy 025711 . 

Considering the differential energy spectrum of electrons given by equation (1.4) the 
exponent of the 7-ray energy spectrum will be given by — (F + l)/2. 

Curvature radiation 

Very high energy particles (e^) radiate curvature radiation in their direction of mo- 
tion 022311 when forced to move along the curved magnetic field lines in very intense 
magnetoshere of pulsars. The characteristic energy of the radiation is given by 



Ec (eV) 



5 ^,3 



3hcY _ 2.96 X 10-^7 
4 7rpc 



Pc (cm) 



(1.11) 



where is the radius of curvature of the magnetic field line and 7 = Ep/rripC^. A 
10^^ eV electron moving along a field line with a curvature of 10^ cm, typical for a 
pulsar, emits photons of energy ^2.5 GeV . 

7r° decay 

Protons can produce tt mesons either in inelastic collisions with matter or interacting 
with ambient radiation, e.g. the microwave background. Both charged (7r+, tt^) and 
neutral pions (tt^) are produced with the same probability in the nucleonic cascades 
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and are the lightest (771,^0 = 135 MeV, ■m^± = 140 MeV) mesons. Neutral pions with 
a mean life of 0.83 x 10"^^ 5x7^ (where 7^ is the Lorentz factor of the pion) decay 
into 7-rays . 

While interacting with matter, the threshold kinetic energy, for producing a particle of 
mass m in a proton-proton collision is given by 

771 

Tth = 2mc\l + —) (1.12) 
nip 

where rup is the proton mass. To produce a 7r° meson, the threshold kinetic energy is 
279.6 MeV. This process produces a very weak 7-ray flux due to the steep spectrum 
of the CRs (^ E'^-'^) and its cross section (app ^ 2.5 m barns). 

The second process of VHE 7-ray generation is the interaction of proton with radia- 
tion which has relatively small cross section (^ 250 /i barns) but is more promising. 
The threshold energy for the production of 7r° from photons of 2.7 microwave back- 
ground is ~ 10^° eV [I257II and is around 10^'^ eV for 1 keV photons. In the vast 
intergalactic distances this process leads to sufficient energy degradation and hence to 
a steepening of the primary CR proton spectrum at energies > 10^° eV. The resulting 
7-rays have energies in excess of ^ 10^^ eV and interact with the same microwave 
background photons to produce pairs. These particles further produce high energy 
7-rays by IC scattering process. The resulting electromagnetic cascade gives 7-rays of 
lower energies [l293] , |27l]| • 

1.5.2 Absorption of VHE 7-rays 

The absorption mechanisms of 7-rays can be divided into two parts: 7-matter interaction and 
7-7 interaction. 

• Gamma-matter interaction 

There are three processes by which 7-rays can interact with matter. (?) Photoelectric 
effect: this process has importance at low energy (hard X-ray region). (//) Compton 
scattering: it is the dominant process in the 7-ray energy range from few hundred keV 
to 10 MeV. {Hi) Pair production: this is the most important energy loss mechanism 
for HE and VHE 7-ray with matter. The incident 7-ray is completely annihilated, i.e. 
hu —7- e+ + e^. The interaction takes place in the electric field of a nucleus and thresh- 
old for the interaction is > 1.02 MeV (i.e. 2 nie c^). It can also take place in the field 
of an electron (a rare phenomenon) but the cross section is much less and the threshold 
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is higher (4 mg = 2.044 MeV). Note that the kinetic energies are not equally shared 
between electron and positron and the initial trajectory is not necessarily the mean of 
their emission angles. The pair production becomes a dominant process above en- 
ergy of about 30 MeV. It is the key process in HE space telescopes and also plays a 
vital role in the development of the atmospheric electromagnetic cascades that make 
VHE 7-ray astronomy possible. The mean distance that a 7-ray photon travels before 
undergoing pair production is 

V = ^ (1.13) 



where N is the number of target nuclei per unit volume, app is the cross section for 
pair production whose value is given by 



^0 



28 / 183 \ 2 

In 



9 V y 27 



cm^ atom ^ (1-14) 



where ctq = (1/137) {e^/{mlc^)) = 5.8 x 10'^^ cm"^ / nucleus = 0.58 mbarn [|295l. 
At very high energies the cross section is independent of energy. The mean free path 
for pair production is related to the radiation length, Xq [discussed later] by 



9 
7X0 



Kp = ^ (1-15) 



Gamma-gamma interaction 

This process requires very unusual combination of high energy photons and a high 
density of low energy photons. 7-ray s are absorbed by photon-photon pair production 
(7 + 7 — e+ + e^) on background photon fields provided the center of mass energy 
of the two photon system is more than the square of rest mass energy of electron. A 
7-ray photon of energy ei collides with another photon of lower energy 62 and gives 
rise to a pair of particles, each of mass m, if ei is greater than a threshold value e*, 
given by [|257]| 

et = ^l m (1-16) 

€2 (1 — COSU) 

where 9 is the angle between the photon trajectories. Taking a case of head on col- 
lisions, e2et = 0.26 x 10^^ (all energies are in eV). Following are the cases with 
astronomically important radiation fields, with their mean energy and required thresh- 
old energy of incident photon [I257II 
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2.7 KCMBR (6 x 10'^ eV) 
starlight (2 eV) 
X-ray (1 keV) 



The cross section for this process is not very small but such collisions are not very fre- 
quent because of the low densities of target photons. Nevertheless, even with low pho- 
ton densities these collisions lead to significant attenuation of photons traveling over 
typical extragalactic distances. The effect is particularly important where the presence 
of EBL, limits the distance to which VHE 7-ray telescopes can detect sources. 



et ^4 X W^eV 
et ^ 10^^ eV 
et^3x 10^1 eV. 



1.6 Extragalactic background light (EBL) and 7-ray 
horizon 

Diffuse extragalactic background radiation has been observed from radio to high energy 7- 
rays . The faint point sources are the main contributors for this radiation and hence it is 
a good tool to study global parameters of source populations and universal physics. The 
optical to infrared diffuse radiation, also called EBL, carries important information about 
galaxy and star formation history of the universe, and is produced by direct starlight (UV 
and visible range) and starlight that has been absorbed and reemitted by the interstellar dust 
(IR to sub-mm range) [I197II . Accordingly, a two peak structure (first peak around 1 — 2/^m 
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Fig. 1 .3: Illustration of absorption of extragalactic VHE photons by EBL traveling over large dis- 
tances. This picture has been taken from [ 191 . 



and second in the wavelength range ~ 100 — 200 fxm) is commonly expected in the spectral 
energy distribution (SED) of EBL [|37l (Fig. 1.2). The SED of EBL is not known accurately, 
especially in the near and mid infrared range. Upper limits and absolute measurements come 
from direct observations which might be polluted by foreground emission, while indirect 
upper limits can also be set by observations of high energy 7-ray sources. Galaxy number 
counts integrations of observable galaxies and missing possible faint sources, give strict 
lower limit II197II . While traveling from the cosmic sources towards the Earth, 7-rays can 
be absorbed in the interstellar or intergalactic space. Attenuation of VHE 7-rays by matter 
is negligible, while they are absorbed (Fig. 1.3| ) by low energy photons of EBL via pair 
production (explained in the last section) and leaves an imprint on the measured spectra 
of different sources. The cross section for this process is described by the Bethe-Heitler 
formula [Ell 



a{E{z), e{z), x) = 1.25 x 10-^5(1 - 



2(3{I3^ - 2) + (3 - (3^) In 
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(1.17) 

(1.18) 



Eex{l + zf _ 

me being the mass of electron, E is the energy of a 7-ray photon, e is the energy of EBL 
photon, z is the redshift of the source and x = (1 — cos9). The cross section is maximum 
when /3 = 0.7 which corresponds to e x ~ 0.5 x 10^^ V, where e and E both are in units of 
eV. Note that while computing the optical depth the energy of a photon scales with redshift 
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Fig. 1 .4: Illustration of 7-ray horizon which can be constrained by observations upto z=0.536. The 
prediction range of EBL models is illustrated by thick solid black line and dashed-dotted 
blue line. Dashed blue line represents an upper limit of EBL based on 3C279 data taken 
with MAGIC telescope. Grey area indicates r > 1 i.e. flux of 7-rays is strongly suppressed. 
The image has been adapted from [i74il . 



z as (1 + 2;). For E = 1 TeV, the cross section is maximum at e ~ 0.5 eV, which means 
VHE 7-rays then interact mainly with optical or infrared photons. The cross section with 
CMBR peaks when the energy of 7-ray is ~ 1 PeV. The attenuation suffered by observed 
VHE spectra of sources, can be used to derive constraints on the EBL density [I272L First 
limits on the EBL were given by Stecker et al. n273ll . while the recent measurements through 
the detection of distant AGNs can be found in [|57] . 12291 . 

The probability for a photon of observed energy E to survive absorption while traveling 
from its source at a redshift z is expressed as 013211 

e-^^^''\ (1.19) 



The coefficient t{E, z) is called optical depth. The energy dependence of r leads to sig- 
nificant deviation of the observed source spectrum as compared to spectrum at the source 
even for small differences in r, due to the exponential dependence. Since the optical depth 
increases with energy, the observed SED becomes steeper than the emitted one. 

The 7-ray horizon niOOM also called attenuation edge for a photon of energy E is defined 
as the distance corresponding to a redshift z for which t{E, z) = 1, i.e. the path length 



which provides the attenuation of the photon flux by a factor of 1/e (Fig. 1.4). Once the 
opacity of the Universe for 7-rays is known it is possible to calculate the intrinsic emission 
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Fig. 1 .5: Third EGRET catalog (1999) at energy more than 100 MeV. This figure has been adapted 
from ||7J. 



spectrum of TeV sources and test the emission models. 

Recently detected blazar 3C279 (z=0.536) by the MAGIC collaboration [|74l has raised 
the question of transparency of the Universe to 7-rays . The redshift region over which the 
7-ray horizon can be constrained by observations has been extended up to z=0.536. The grey 



area in Fig. 1 .4 indicates an optical depth r > 1, i.e. the flux of 7-rays is strongly suppressed. 



1.7 Gamma ray sources 

Almost a century after the discovery of CRs, currently we are in the search for the sources 
of these energetic radiations. More than one half of the 271 sources detected at energies 



> 100 MeV in the third EGRET catalog of 1999 (see Fig. O) are still unidentified. Cur 



rently, Cherenkov telescopes have made an unprecedented progress in the field of VHE 7- 
ray astronomy which has resulted in the detection and identification of a large number of 



sources (see Fig. 1.6). To date more than 100 sources (61 galactic and 42 extragalactic) 
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Fig. 1 .6: Known sources in the VHE 7-ray sky as of August 2010. This figure has been adapted 
from d. 



have been detected in VHE 7-rays . The Large Area Telescope (LAT) which is the primary 
science instrument on the Fermi Gamma-ray Space Telescope (Fermi) started surveying the 
sky in 2008. After the first 11 months of the science phase of the mission, the first source 
catalog has been released. This catalog contains 1451 sources (see Fig. |1.7| ) detected and 
characterised in the 100 MeV to 100 GeV energy range [|26l . 



1.7.1 Galactic sources 
Supernova remnants 

When a star reaches its final stage of evolution by burning all its fuel necessary for the fusion 
reactions that counteract the gravitational pressure, the core collapses followed by ejection of 
material into outer space. This catastrophic explosion of a star known as a supernova, emits 
~ 10^^ ergs of energy in a very short time interval [|295|| . The remains of the destroyed 
star can be detected for thousands of years over a wide range of energies. The 7-rays that 
result from supernova explosions may be detected instantaneously (within few seconds of 
the explosion) as a GRB, as the steady but periodic emission from pulsars (rotating core) or 
as the expanding outer shell of the star known as the supernova remnant (SNR). 
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Fig. 1 .7: Fermi LAT first catalog contains 1451 sources detected in the energy range 100 MeV to 
100 GeV. This figure has been adapted from [t26il . 

Shell type SNR 

This is the most common type of SNR, characterized by the shell of interstellar material 
swept up by the expanding ejected material and shock wave developed along the way, which 
is clearly visible in X-rays. Particles are accelerated in these shocks by their coUisionless in- 
teractions with magnetic clouds in the interstellar medium (Fermi acceleration). For decades 
these sources have been believed to be the principal objects for the production of galactic 
CRs, based on the kinetic energy available in supernova events. Nearly 10% of the total 
mechanical energy released by the supernova explosions in our Galaxy is required to main- 
tain the observed density (10~^^ ergs cm^^) of the galactic CRs ll63ll . However it should be 
noted that some other sources (e.g. pulsars or X-ray binaries) may also meet this energy re- 
quirement. The theory of diffusive shock acceleration II222II provides acceleration efficiency 
(the fraction of the mechanical energy of the shock transferred to the non-thermal particles) 
as high as 10% and can explain very hard power law source spectrum (a ^ 2. 1) inferred from 
CR measurements. The high acceleration efficiency and hard energy spectra of protons (up 



Chapter 1 . Very high energy gamma-ray astrophysics 



19 



to multi-TeV energies) should lead to production of VHE 7-rays of hadronic origin (from 
7r° decay). Almost all young SNRs are emitters of non-thermal X-rays, presumably origi- 
nating from the synchrotron radiation by electrons of energy ^ 10 TeV. These electrons 
can also radiate TeV 7-rays , therefore SNRs are expected to radiate in TeV energy by two 
competing processes. Recently deep survey of the galactic plane conducted by HESS firmly 
established the young shell type SNRs as copious TeV 7-ray emitters. Examples of shell 
SNRs (emitting TeV 7-rays ) are Cassiopeia A, RX J1713.7-3946, RX J0852.0-4622 etc. 

Pulsar Wind Nebulae (PWN) or Plerions 

PWN or plerions are the SNRs which consist of swept up material of supernova explosion 
and have a rotating magnetized neutron star (pulsar) at their center. The swept up material 
called nebula has strong non-thermal radio to 7-ray emission. Generally accepted magneto- 
hydrodynamic (MHD) model presumes that spin down energy of the neutron star is carried 
away by the stream of particles (mechanism is not completely understood), mostly electron 
and positron which are produced by the pulsar. The particles in the pulsar wind interact 
with magnetic field of the nebula and emit radiations from radio to X-ray via synchrotron 
process. The pulsar wind encounters the interstellar matter at the edges of the nebula which 
ends in a shock and can accelerate the particles up to PeV energies. These accelerated par- 
ticles interact with their own synchrotron photons or CMBR and can efficiently emit VHE 
7-rays through IC process. They currently constitute the most populated class of identified 
galactic VHE sources. As compared to the shell type SNRs the upscattering process in these 
objects seems to be very efficient which explains their abundance in the Galaxy. The Crab 
Nebula a prototype object of this class is the brighest known steady VHE 7-ray emitter (the 
standard candle source) and was the first source discovered in the VHE regime II296II . 

Pulsars 

Pulsars are fast spining magnetized neutron stars produced in supernova explosions. They are 
divided into two categories: rotation powered (detectable at radio wavelengths) and accretion 
powered (detectable at X-ray energies - generally thermal). Another category is known as 
millisecond pulsars. It is generally believed that rotation-powered pulsars are formed when 
cores of massive stars (> 8 Mq) collapse [I220II . Pulsars have their masses ranging from 
1.4 M0 to 3 Mq, diameter of about 10 km, core density reaching upto 10^^ g/cm^ and can 
produce the strongest magnetic field in the universe, of the order of 10^^ The rotational 
^1 G= lO-'^T 
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energy is given by 

K=^Icu^ = 27r^^. (1.20) 

Here symbols have their usual meanings. For M = 1.4 Mq and i? = 10 km, moment of 
inertia I = (2/5) M i?^ = 1.1 x 10^^ g cm^. The rate of rotational energy loss (called spin 
down luminosity or spin down power) can be written as 

dK ^ , I dP 

The period P can range from a millisecond to a few seconds. Most of the pulsars show 
gradual increase in the spin period, for the Crab pulsar, P = 33.3 ms, dP/dt = 4.21 x 
10"^^ s and dK/dt = 5 x 10^^ ergs s"^ 

The strong co-rotating magnetic field causes a strong electric field of the order of 6 x 
10^° V m'^ HI 6711 at the poles of the neutron star, which accelerates the particles near the 
pulsar. The rotation and magnetic axes of pulsars are misaligned, we see the emission only 
when the beam crosses our line of sight. The emission mechanism is not yet completely 
understood, however it is generally accepted that the primary radiation mechanism is the 
synchrotron-curvature radiation taking place near the polar caps of the pulsars. The sec- 
ondary possible mechanisms are ordinary synchrotron and IC scattering in the outer region 
of the magnetoshere. 

It is not known whether the emission of radiation takes place near the magnetic poles 
of the rotating neutron star (the polar-cap scenario 1117211 ) or farther out in the magnetoshere 
(slot- gap or outergap scenario [[80] , I116B and references there in). The high energy part 
of the 7-ray spectrum differs substantially between near and far cases. Again the slot gap 
and the outer gap have different predicted 7-ray spectrum. The detection of 7-rays above 
lOGeV would allow to differentiate between the different pulsar emission models. The 
energy spectra and the structure of light curves in the region from several GeV to 30 GeV 
carry information about the location of 7-ray production region in the pulsar magnetospere. 
Generally due to pair production in the strong magnetic field, no TeV 7-ray emission is 
expected from the pulsar magnetoshere [1173! , 12631 , 1186L Recently MAGIC collaboration 
reported the detection of 7-rays from the Crab pulsar above 25 GeV [fTTil which favors the 
outer gap model. 

Only six pulsars were firmly established as the sources of high energy 7-rays by the 
EGRET on board CGRO launched in 1991 [|295l . Recent resuks from the Fermi-LAT col- 
laboration reported the detection of 30 7-ray pulsars, out of which 15 were not known earlier 
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at other wavelength Il24l . 
Compact binary systems 

Most of the stars in our galaxy are associated with other stars making multiple star systems 
which are orbiting around each other. Half of the stars in the steller population occur in 
some kind of binary association [I295II . The compact binary system consists of a compact 
star generally a neutron star or black hole and a normal massive star called the companion 
star. The compact star accretes mass from the companion star. They are known for the 
powerfuU emission of temporally and spectrally variable X-ray radiation. At the moment 
four binaries are recognised VHE 7-ray emitter: PSR B 1259-63 LS 5039 LS I-1-6I 
303 [[67ll and Cygnus X-1 [TtTI . Particle acceleration is possible in these systems either due 
to the termination of pulsar wind or through the internal shocks in the jets formed in the 
vicinity of the black hole. The 7-rays are then produced due to the presence of the dense 
target material in the form of optical photons or gas provided by the companion star. 

X-ray binary systems 

X-ray binaries are a class of binary stars that are luminous in X-rays and are generally treated 
as thermal sources. The X-rays are produced by matter falling from the companion star to the 
compact object which effectively transforms the gravitational energy of the compact object 
into X-ray emission radiated by hot accretion plasma. The most likely emission process is 
the injection of relativistic electrons into the surrounding medium and further acceleration of 
these electrons by the shock region created between winds of the two objects. VHE 7-rays 
are the result of IC scattering of the electrons with the ambient photons. 

Presently, four compact binary systems are found to emit VHE 7-rays . Three of these bi- 
nary systems: PSR B1259-63 153, LS 5039 ^ and LS I-1-6I 303 [|67l are clearly associated 
to X-ray binaries [I245II . These three Binary TeV sources (BTV) have a bright high-mass pri- 
mary star, which provides intense UV seed photons for IC scattering of particles accelerated 
around the compact object. All of them have been detected at TeV energies in several parts 
of their orbits and show variable emission and hard spectrum. The emission is periodic in LS 
5039 and LS I -1-6I 303 with periods of 3.9 and 26.8 days respectively [l56|,|75l and consistent 
with their orbital periods [[79] , lllll , 11121 . These two sources also share the distinction of 
being the only two known high energy emitting X-ray binaries that are spatially coincident 
with sources above 100 MeV listed in the Third EGRET catalog [I174L Both these sources 
have also been recently detected by the Fermi observatory ||23] . l25l . PSR B 1259-63 was not 
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detected by EGRET and its periodicity has also not been determined yet, as the orbital pe- 
riod is quite large (~3.4 years) which requires the monitoring of the source for an extended 
period (few years). 

It is believed that the BTV Cygnus X-1 has a black hole as a compact object where as in 
the case of PSR B 1259-63 it is a neutron star. The nature of compact objects in other two 
binary systems LS 5039, LS I-I-61 303 is still not completely understood as there is no strong 
evidence as yet in support of either a black hole or a neutron star. HESS telescope array 
discovered a new BTV candiadate HESS J0632-I-057 which appears to be a point like source 
within experimental resolution f6T\. It is possibly associated with the massive Be star MWC 
148 and has been suggested to resemble known TeV binary systems like LS I -1-61 303 or 
LS 5039. No evidence for flux variability from this source was found by the HESS group. 
The VERITAS collaboration has also observed this source for 31 hours in 2006, 2008 and 
2009 with no significant 7-ray signal above ITeV being detected from it [|30ll . This non- 
detection excludes with a probability of 99.993% that the HESS J0632-I-057 is a steady 7- 
ray emitter. Contemporaneous X-ray observations with the Swift X-Ray telescope reveal a 
factor of 1.8±0.4 higher flux in the 1-10 keV range than earlier X-ray observations of HESS 
J0632-I-057. Future multiwavelength observations combined with results from ground-based 
and space-based 7-ray observatories will provide a deeper understanding of the true nature 
of HESS 10632-^057. 

Microquasars 

A microquasar (or radio emitting X-ray binary) is a smaller cousin of a quasar. Microquasars 
are named after quasars, as they have some common characteristics: strong and variable 
radio emission, often resolvable as a pair of radio jets, and an accretion disk surrounding 
a compact object. In quasars, the black hole is supermassive (millions of solar masses); in 
microquasars, the mass of the compact object is only a few solar masses. In microquasars, 
the accreted mass comes from a normal star, and the accretion disk is very luminous in 
the optical and X-ray regions. Microquasars are sometimes called radio-jet X-ray binaries 
to distinguish them from other X-ray binaries. A part of the radio emission comes from 
relativistic jets, often showing apparent superluminal motion. 

Microquasars are very important for the study of relativistic jets. The jets are formed 
close to the compact object, and timescales near the compact object are proportional to the 
mass of the compact object. Therefore, ordinary quasars take centuries to go through varia- 
tions which a microquasar experiences in one day. MAGIC observed Cygnus X-1 recently 
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during a short lived flaring episode and reported strong TeV emission [Jlj . This high mass 
X-ray binary is a microquasar [I129II . 

Young open star clusters 

Star clusters which are groups of stars are of two types. One is the globular clusters which 
are tight groups of hundreds of thousands of very old gravitationally bound stars, while other 
is the open clusters, a more loosely clustered group of stars, generally contain less than a few 
hundred members, and are often very young. The particles can be accelerated to extreme 
energies by means of shock created by the collision of winds of these massive young stars. 
HEGRA reported the first evidence of this new class of galactic VHE source in the vicinity of 
Cygnus 0B2 region which was confirmed by MAGIC. HESS detected VHE emission from 
a massive star called Wolf-Rayet (WR) star ( 1112911 and references therein). 

Galactic center 

The galactic center (1=0°, b=0°), located at a distance of 8.33±0.35 kpc from the Earth [I164II , 
is a region of interest and it is significant that there is an EGRET source 3EG J 1744-3039, po- 
sitionally coincident with it. Fermi, which is currently surveying the sky in the energy range 
between 20 MeV to more than 300 GeV with unprecedented sensitivity and resolution has 
detected a source at the galactic center with high significance after 3 months of observation, 
which is listed in the published LAT Bright Source List [|2l as OFGL 1746.0-2900. It is be- 
lieved that there is a Supermassive Black Hole (SMBH) at the center of our Galaxy (as most 
if not all galaxies are believed to contain SMBHs at their centers). The galactic center emits 
TeV 7-rays as reported by the Whipple BM, CANGAROO [|28l, HESS Il5n,[58],l60ll and 
MAGIC [|65ll collaborations. The 7-ray source shows steady emission even during the X-ray 
flares and behaves like a point source within the angular resolution of the current instruments. 
Since the region around the galactic center contains many potential 7-ray sources (the young 
sheU type SNR Sag A East, the newly discovered PWN G 359.95-0.04, the SMBH Sag A*), 
the exact interpretation of the emission process is rather difficult [I132L The 7-ray production 
in the vicinity of SMBH is expected to be variable on time scales as short as 1 hour. The 
lack of variations in TeV flux therefore interpretated as 7-ray production in extended re- 
gions however it can not be used as a decisive argument against the black hole origin of TeV 
emission [|63l . Other possible emissions because of SMBH could be acceleration of elec- 
trons at termination shock of its hypothetical wind (similar to PWN but powered by black 
hole). Plausible radiation mechanism therefore includes IC scattering of energetic electrons. 
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the decay of pions produced in the interactions of energetic hadrons with the interstellar 
medium or radiation field and curvature radiations of UHE protons close to SMBH. These 
considerations disfavour dark matter annihilation as the main origin of the detected flux. 

The possibility of indirect dark matter detection through its annihilation into VHE 7- 
rays has aroused interest to observe the center of our galaxy. Recently HESS and MAGIC 
observed the galactic center, measuring a steady flux consistent with a differential power 
law (spectral index of about 2.2) upto energy ~ 20TeV with no apparent cutoff [ISTI . [65l. 
The detection of TeV 7-rays from the galactic center (within ^1 arc minute around Sgr 
A* (a hypothetical SMBH) is the first model-independent evidence of high energy particle 
acceleration in the central 3 pc region of the galaxy. There are several objects located in this 
region which can be responsible for VHE emission, in particular the dark matter halo, the 
young sheU type SNR Sag A East, the newly discovered PWN G 359.95-0.04, the SMBH 
Sag A* itself and the entire central 1 pc diffuse region filled by dense molecular clouds and 
CRs, are possible counterparts of HESS J1745-290. 

Unidentified sources 

There exists a population of TeV 7-ray sources which does not have counterparts at longer 
wavelength i.e. in the X-ray or radio energy bands which spatially coincide with the location 
of the observed excess. These sources cluster around the galactic plane, indicating that they 
are likely of galactic origin. Most of these sources are spatially extended sources, further 
suggesting that they might be associated with unseen PWNe or SNRs which constitute the 
largest population of galactic TeV 7-ray sources. Many previously unidentified TeV gamma- 
ray sources have subsequently been found to be probably associated with PWNe or SNRs 
( 1129711 and references there in). 

The situation for being unidentified is the likely outcome of experimental and physical 
considerations. As an example of the former, extended sources may have different mor- 
phologies at different wavelengths while for the later, sources may emit substantially in the 
GeV — TeV range and are too dim to be detected at other spectral bands ("dark accelera- 
tors") e.g TeV J2032-I-4130. The explanation for these unidentified sources may either be 
exotic where these sources originate from the annihilation of dark matter or conventionally 
from p-p interaction which entails emission mostly in the VHE, contrary to electrons which 
produce comparable fluxes by synchrotron emission. 

TeV J2032-I-4130 is the first unidentified TeV 7-ray source and has to date remained 
as such. This source was first discovered by HEGRA [,43 J with marginal statistical signif- 
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icance initially but latter on confirmed by a number of experiments [|55] , 11991 , |72l|- The 
HEGRA results indicated that this source was an extended TeV source and the measured 
7-ray spectrum was hard with a power law index of 1.9. These results have been confirmed 
by the MAGIC collaboration. The MILAGRO group had discovered three unidentified VHE 
7-ray sources MGRO J2031+41, MGRO J2019+37 and MGRO J1908+06 ||22|,|2ll, among 
them the last one is still unidentified though the other two have now been identified at other 
energies ( 1113211 and references therein). A catalog of bright 7-ray sources which has been 
released recently contains all the sources detected by Fermi LAT with a statistical signifi- 
cance of more than lOcr [|24ll . The Fermi LAT source OFGL J2032. 2-1-4122 lies very close to 
CoG of TeV J2032-I-4130 (about 8 arc min) and therefore is a promising candidate for being 
the GeV counterpart of the latter, based on spatial coincidence alone. Interestingly OFGL 
J2032. 2-1-4122 is also one of the 29 7-ray pulsars detected by the Fermi LAT. This source is 
named as LAT PSR J2032-I-41, because it falls in a special category of pulsars which seems 
to show pulsation only in the 7-ray band. Being "dark" at other wavelengths, these sources 
could have easily escaped previous pulsar surveys. 

Diffuse emission 

The diffuse galactic 7-ray emission is produced by CR protons, nuclei and electrons inter- 
acting with the interstellar gas and the radiation fields. It is generally believed (still not 
experimently proved) that CRs upto the knee point (< 5 x 10^^ eV) are produced in SNRs. 
The spectrum of these 7-rays arising from the interaction of the CRs with matter in the galaxy 
can provide important information to understand the origin and propagation of CRs in our 
galaxy. The measured differential spectral index of CRs is a ~ 2.75 11821 . while shock accel- 
eration theory predicts a spectrum close to a ~ 2 112221 . [T45l . The difference in the spectral 
index is attributed to diffuse energy loss of higher energy CRs from our galaxy I1184II . Mea- 
surements by EGRET indicated an excess of 7-ray emission at energies > 1 GeV relative 
to diffuse galactic 7-ray emission models that were consistent with directly measured CR 
spectra in the Earth's neighbourhood. However, the measurements of galacic diffuse emis- 
sion with Fermi LAT for energies 100 MeV to 10 GeV do not show the GeV bump II285II 
discovered by EGRET. 
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1.7.2 Extragalactic sources: 
Active Galactic Nuclei 

AGN are compact regions at the center of galaxies which have a much higher than normal 
luminosity over some bands or all of the electromagnetic spectrum. A galaxy hosting an 
active nucleus is called an active galaxy. AGN host a very compact SMBH with masses 
between 10^ to 10^*^ solar masses. In the standard model of AGN, cold material close to the 
central black hole forms an accretion disc. Dissipative processes in the accretion disc trans- 
port matter inwards and angular momentum outwards, while causing the accretion disc to 
heat up. At least some accretion discs produce jets, twin highly coUimated and fast outflows 
that emerge in opposite directions from close to the disc and perpendicular to it. The jet 
production mechanism and indeed the jet composition on very small scales are not known at 
present, as observations cannot distinguish between the various theoretical models that exist. 
The jets extend up to several kpc distance scale from the central SMBH and have the most 
obvious observational effects in the radio waveband. The jets contain excellent conditions 
for efficient particle acceleration and subsequently highly variable emission of non-thermal 
radiations which has been measured over all the wavelengths. Currently, the AGN classi- 
fication schemes are based on the differences due to anisotropic radiation pattern i.e. the 
pointing directions are more important than the intrinsic physical differences. AGN with 
relativistic jets close to the line of sight of the observer are classified as blazars. More details 
about AGN along with their emission mechanism are discussed in the next chapter. 

Gamma Ray Bursts 

GRBs are flashes of 7-rays associated with the extremely energetic explosions in distant 
galaxies. They are the most luminous electromagnetic emissions occurring in the Universe 
in any wavelength band and perhaps the brightest phenomenon since the big bang. Bursts 
can last from milliseconds to nearly an hour, although a typical burst lasts a few seconds and 
typically (10^^ - 10^^) ergs s^^ are released within this time . The short duration of GRBs 
hints at very compact progenitors. The initial burst mostly in 7-rays is usually followed by 
a longer-lived "afterglow" emitting at longer wavelengths. In spite of many successes in the 
interpretation of the various phases of a GRB |304l|, there are still fundamental unanswered 
questions involving essentially all aspects of the GRB phenomenology. Currently, the most 
accepted picture of GRBs is the cosmological fireball model [13041 , |250l , |230ll . In the fireball 
model the various observed phenomenology from a GRB are due to an ultra-relativistic out- 
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flow with initial Lorentz factor between (10^ - 10^) 023311 generated during the final collapse 
of a high mass star or the coalescence of a binary system made of compact objects. The for- 
mer case is usually associated with long duration GRBs, while the later to the short duration 
class. The VHE 7-rays are expected to be produced in internal shocks due to synchrotron 
and IC processes. The mean distance of GRBs is about z = 2.8, which essentially means that 
in most cases the VHE 7-rays will not be able to reach the Earth because of EEL absorption. 
In the past several detections of GRBs were made in the MeV — GeV range by EGRET 
and more recently by AGILE and Fermi. On the contrary, in spite of continuous attempts, 
so far no confirmed detection at VHE range has been obtained. First upper limit on VHE 
7-ray emission of GRBs has been published by MAGIC collaboration ll64l . 

Star-forming galaxies 

These young galaxies are characterized with an exceptionally high star formation rate as 
compared to the normal galaxies. High star formation and supernova rate in starbust galax- 
ies enhances the energy density of energetic nonthermal particles - mostly electrons and 
protons which are Fermi accelerated in the sites of SNRs. Coulomb, synchrotron and Comp- 
ton energy losses by the electrons and pions decay (after their production by the interaction 
of energetic protons with ambient gas) results in emission from radio to high energy 7- 
rays . An estimation of VHE emission has also been obtained by analytical calculations of 
hadronic process. However, only a weak level of emission is expected making the nearby 
and extremely powerful starburst galaxies like M82, NGC253 and Arp220 the obvious tar- 
gets for observations. The flux predicted from M82 by Persic et al. [I246II could be de- 
tectable with long exposure by MAGIC II and VERITAS, whereas Cherenkov Telescopes 
Array (CTA) would detect it in nearly 10 hours of observations. Indeed, VERITAS ob- 
served the M82 galaxy during 2007-2009 for about 137 hours and detected the source ||29ll 
at energies > 700 GeV. This is one of the weakest VHE sources observed, with a luminosity 
about 0.9% Crab. Similarly NGC253 has been reported by HESS to be detected in a long 
campaign during 2005-08 with 119 hours of observation time [f34ll . The source was detected 
at the 0.3% Crab level, which appears to be lower than what models have predicted [I140L 
Observations of Arp220 by MAGIC has placed a loose upper limit to its VHE flux f,68ll . 
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2.1 Introduction 

The fundamental building blocks of the Universe are stars which can also be found in clus- 
ters. The most common star system is the galaxy which has 10^ to 10^^ stars. The "Milky 
Way", is a typical Galaxy with ~ 10^^ stars in a disk like structure. Galaxies themselves are 
found in clusters which are members of so-called superclusters [1295 L A Galaxy is called a 
normal galaxy if it is not ultra-luminous and most of the observed optical radiation seems to 
be the sum of all the steller emissions. The Andromeda galaxy is our closest normal galaxy 
at a distance of about 570 kpc. On the other hand, active galaxies are those whose core lumi- 
nosity exceeds that of the normal galaxies and points to the evidence of relativistic particle 
acceleration II295II . 

In this chapter, we discuss the classification of AGN and models for observed 7-ray emis- 
sion in blazars (a particular class of AGN). We also present a list of extragalactic VHE 7- 
ray sources detected till date i.e. August 2010. 

2.2 Active Galactic Nuclei 

AGNs [I195II unusually have a bright, starlike central core called a nucleus. The nucleus emits 
more radiation as compared to the rest of the galaxy and the emitted radiation shows high 
variability of flux. Because of the enormous energies coming from rather small volumes, it 
is generally accepted that these nuclei have supermassive black holes (usual mass of SMBH 
is of the order of 10® to 10^° solar masses) at their centers as a source of energy for the 
extraordinary luminosity. The observed energy spectra of AGN are bright at all wavelengths, 
indicating the non-thermal nature of emission, unhke stars which show a typical black body 
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Fig. 2.1 : Empirical classification scheme of AGNs adapted from 112901 . 



radiation. Though the 7-ray luminosity also varies at all wavelengths spanning from low 
energy to VHE range, variations at lower energies (radio and optical) are usually weaker and 
slower than at high energies. These objects are characterized by low emission quiescent state 
with occasional outbursts referred to as "flares", which may last for hours to months. 



2.2.1 Classification of AGN 

Detailed description of classes of AGN can be found in [12891 . 12401 . 1241| . Earlier these objects 
were empirically discerned by their radio emission, the properties of their optical emission 
lines, morphological considerations, 7-ray emission and some other properties (Fig. |2.1[ ). 
AGN are primarily divided into two groups of spiral and elliptical galaxies. The spiral galax- 
ies are subdivided into Seyfert galaxies (type I and II) based on their optical emission lines. 
On the other hand, one class of elliptical galaxies with weak radio emission are called Radio 
Quite Quasars (RQQ), while the other with strong radio emitting objects are further classi- 
fied on the basis of their optical emission lines. The objects with strong optical emission 
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Fig. 2.2: An artistic view of the unified AGN scheme. When the viewing angle is quite narrow with 
respect to the jet, system is called a blazar. This figure has been adapted from [,15.1 . 

lines are Radio Quasars (they have two groups: Flat Spectrum Radio Quasars (FSRQs) and 
Steep Spectrum Radio Quasars (SSRQs)). The second category in this class i.e. objects 
having weak or no optical emission lines are further classified as Radio Galaxies having no 
7-ray emission and the BL Lacerate (BL Lac) objects which have 7-ray emission and also 
show strong flux variability. The FSRQ which also have 7-ray emission, strong variability 
and optical polarization, together with BL Lacs form the class of Blazars. All blazars emit 
7-rays . 

2.2.2 The unified sclieme of AGN 

A unified theory of all AGN has been proposed [12891 . which reduces the differences in the 
observed properties to mainly the orientation of the sources relative to our line of sight. In 
this theory the central SMBH 112 1911 accretes surrounding material forming a rotating plasma 
accretion disk, which is heated up and emits thermal radiations from optical to X-ray wave- 
length. The resulting spectrum peaks in X-ray band. The gravitational energy of infalling 
matter is transferred into thermal radiation with a conversion efficiency of about 10% [11291 . 
The rapidly moving gas clouds close to the central engine (SMBH) or above the disk also 
heat up and form a broad-line emission region, observed by the Doppler-broadened emission 
lines. This strong optical and ultraviolet radiation is obscured along some lines of sight par- 
ticularly along the equitorial plane by a dusty torus that surrounds the whole system of AGN. 
Beyond the torus i.e. farther from the central source, slower moving molecular gas clouds 
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Fig. 2.3: Another view of the AGN unified scheme wherein the classification has been made on the 
basis of viewing angle. The figure illustrates the different regions around the central black 
hole. The image has been adapted from ||96l . 



produce narrow line absorption and emission spectra. The broad line region has a particle 
density of > 10^ particles cm^^ and line widths of a few hundred thousand km/ s, while 
the narrow line region has a lesser density of 10^ to 10^ particles cm~^ and line widths of 
a few hundred km/s [|165l . l290ll . The spin of the black hole induces twisted magnetic field 
lines as jets at the poles of the rotating source fil03i . Particles can be accelerated by strong 
shock waves in the jets to ultra-relativistic energies. The jets may extend up to 1 Mpc for 
FR-II type AGN [fTTl . 

For a 10® Mq black hole, the radius is ~ 3 x 10^^ cm, the broad line clouds are located 
within ~ 2 — 20 X 10^^ cm of the black hole and the inner radius of dusty torus is perhaps 
~ 10^^ cm. The narrow line region extends approximately from 10^® — 10^° cm and radio 
jets have been detected on scales from 10^^ to several times 10^° cm [I289L 



A schematic drawing of the unified model can be seen in (Fig. 2.2 and 2.3 1. This model 
suggests that the wide variety of AGN types seen is largely a function of viewing angle , and 
hence of geometry rather than physics. It is generally believed that the jets are filled with 
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relativistic particles and that much of the observed radio emission is synchrotron radiation 
from these particles. The thick torus which surrounds the accretion disk may completely 
obscure it, if the system is viewed edge-on. The radio lobes are then clearly seen and the 
AGN is called a radio galaxy (with radio emission dominant). 

If the line of sight is far away from the jet (almost perpendicular), the optical and radio 
emission from the core can be seen along with broad emission lines. Such AGN systems 
are known as lobe dominated quasars (type II quasars and Seyfert II galaxies are generally 
equitorial views). The core becomes more and more apparent with decreasing viewing an- 
gle and the objects are classified as core dominated quasars (type I quasars and Seyfert I 
galaxies are intermediate views). The synchrotron emission from jets becomes the strongest 
mechanism and radiation is variable and polarized [I295II . When the viewing angle becomes 
quite narrow, the jet is the most obvious feature and the system is called a blazar. If the ob- 
server is directly looking into the jet (essentially zero viewing angle), it is an extreme blazar. 
Blazars with no emission lines are classified as BL Lacs, while the ones with emission lines 
are FSRQ. 

2.3 Blazars 

Currently the blazars constitute the largest class of known Te\^7-ray emitters. They are 
characterized by a close alignment of the jet axis with the line of sight of the observer. The 
blazar types have differences in the optical emission, as Optically Violently Variable (OVV) 
quasars and FSRQs have strong emission lines while BL Lacs show weak / no optical lines. 
These differences are probably due to the amount of cold molecular clouds encountered 
along the line of sight. The relativistic jets which are believed to be the key elements of 
the observed emission of blazars, produce highly variable component of non-thermal emis- 
sion. The dramatically enhanced fluxes of Doppler-boosted radiation (cx 6^) [|63l . resulting 
from the fortuitous orientation of the jets towards us, make these objects ideal laboratories 
for studying the underlying physics of AGN jets through multiwavelength studies of tem- 
poral and spectral characteristics of radiation from radio to VHE 7-rays . First of all this 
concerns the BL Lac objects, several of which are already established as TeV 7-ray emit- 
ters. The rapid variability of these objects in brightness and in polarization is explained by 
relativistic synchrotron jets which indicates the presence of ultrarelativistic electrons. The 
TeV radiation not only tells us that the particles (electrons and/or protons) in these objects 
are accelerated to VHE, but also provides the strongest evidence in favor of the commonly 
accepted paradigm that the non-thermal radiation is produced in relativistic outflows (jets) 
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with Doppler factor 5 > 10 [l37ll. 

Indeed, the enormous apparent VHE 7-ray luminosity of TeV blazars, reaching 

Lapp = A7i£f^ ~ 10^^ ergs/s (2.1) 

(where d = cz/ Ho is the distance to the source and is the integrated 7-ray energy flux) 
during the strongest flares of Mrk421 and MrkSOl can be reduced to more reasonable level 
assuming that the radiation is produced in a source moving with a large Doppler factor, 

Lint ~ S '^Lapp. (2.2) 

This assumption also allows a more comfortable (larger) size of the emitter based on the 
observed variability timescales, 

R < cAUarS (2.3) 

In fact the same arguments apply to other components of non-thermal variable radiation as 
well, but the uniqueness of 7-rays comes from their "fragility" due to interaction with ambi- 
ent photon fields, and the assumption about the relativistic bulk motion becomes unavoidable 
in order to overcome the problem of catastrophic internal 7 — 7 absorption. 

Observations of multi-Tey 7-rays from MrkSOl, the spectrum of which in the high state 
extends to ~ 20 TeV, provide a very robust constraint on the Doppler factor of the jets. 
Assuming that the optical/UV flux of this source in a high state ~ 5 x 10^^^ erg / cm? s 
II249II . is produced in the jet, the catastrophic absorption of 20TeV 7-rays can be avoided 
only when 5 > 10. 

The flux variability on different time scales is a remarkable feature of TeV radiation 
from BL Lac objects. Since the discovery of TeV 7 radiation from Mrk421 {z = 0.030) 
and MrkSOl (z = 0.034) by the Whipple collaboration [12551 . l256ll . these sources have been 
the subject of intensive studies through multi-wavelength observations. Both sources are 
variable with typical average TeV fluxes between 10% to 50% that of the Crab Nebula, but 
with rapid variability on time scales of less than 30 minutes [|156Ll45ll . Recently the MAGIC 
collaboration detected variability of Mrk 501 on timescales less than a few minutes ITTOI . 

lACTs are well suited to search for short duration signals from blazars in the flaring 
states [|63l . This is a key condition which makes the simultaneous observations of X-ray 
and 7-ray wavebands particularly meaningful [11221 . This has been demonstrated in many 
observations of Mrk421, Mrk501, lES 1959-1-650, PKS 2155-304 and some other BL Lacs 
simultaneously with ground based 7-ray and satellite based X-ray telescopes. One of the 
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best results in this area has been obtained during the spectacular flares of Mrk421 in 2001. 
A nice sample of such events, the 19 March 2001 flare detected by RXTE and VERITAS, 
shows a clear keV/TeV correlation on sub-hour timescales H 15411 . The correlated flux vari- 
ability is often interpreted as a strong argument in favor of the synchrotron-Compton jet 
emission models in which the same population of ultra-relativistic electrons is responsible 
for production of both X-rays and TeV 7-rays via synchrotron radiation and IC scattering, 
respectively. However, the very fact of correlation between X-ray and TeV 7-ray fluxes does 
not yet rule out other possibilities, including the so called hadronic models which assumes 
that the observed 7-ray emission is initiated by the accelerated protons interacting with the 
ambient gas, low-frequency radiation and magnetic fields [|37l . 

Since the intrinsic source 7-ray spectrum can be modified by both internal and extra- 
galactic absorption of TeV 7-rays , the identification of radiation processes becomes very 
difficult. As long as the measurements of CMBR remain highly uncertain at all relevant 
wavelengths, the information contained in the observed 7-ray spectrum can not be reliable. 
On the other hand, the spectral variability of TeV emission as well as the correlation of ob- 
served TeV fluxes with other energy bands do not depend on the extragalactic absorption. 
Therefore, only detailed studies of spectral evolution of TeV 7-rays and their correlation 
with X-rays from several X-ray selected BL Lac objects in different states of activity, and lo- 
cated at different distances within 1 Gpc, will allow us to follow and resolve simultaneously 
the fluctuations predicted by different models on timesscales close to the shortest one likely 
in these objects [|37ll . 

The spectral studies of MrkSOl and Mrk421 performed during quiescent or moderately 
high states generally show steep energy spectra with photon indices close to F = 3. However 
the observations of these objects by Whipple, HEGRA and CAT groups during the excep- 
tionally bright and long lasting activity of MrkSOl in 1997 and Mrk421 in 2001 revealed 
very hard energy spectra fl63|- The high state time-averaged 7-ray spectra of both sources 
can be described by a canonical "power law with an exponential cutoff" function of the form 

dN/dE = KE^ e^-^'^""^ (2.4) 

with r ~ 2. The high states of these sources are described by different [F, Eq\ combinations: 
r = 1.92 and Eq = 6.2TeV for MrkSOl 1997 flare and T = 2.19 and Eq = 3.6TeV for 
Mrk421 2001 flare ( ||63l and references therein). The exponential cutoff in the spectrum of 
Mrk421 starts earlier. Since both sources are located at approximately the same distances, 
the difference in cutoff energies can be interpreted as an indication against the hypothesis 
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Fig. 2.4: The average SED of the blazars studied by Fossati et al. II153II luminosity, including the 
average values of the hard X-ray spectra. The thin solid lines are the spectra constructed 
following the parameterization proposed in II142II . This figure has been adapted from [1421. 

that attributes the cutoffs to the pure extragalactic absorption. 

Therefore one may conclude that for unambiguous identification and deep understanding 
of acceleration and radiation mechanisms in the jets, the correlation of absolute 7-ray and 
X-ray fluxes is very important, but not yet sufficient. The key information seems to be 
contained in the spectral variability in both energy bands on time scales comparable to the 
characteristic dynamical times of about an hour or less. 

2.3.1 The blazar sequence 

The observed SED of blazars exhibits a two-bump structure with the first peak called syn- 
chrotron peak in the infrared to keV energy range and the second peak or Compton peak in 
the range of MeV to GeV energies. The statistics of spectral characteristics of the detected 
blazars allows us to make several phenomenological studies about their properties. A popular 
scheme to unify the blazar population is the so called blazar sequence proposed by Fossati 
et al. II153II . The authors analyzed the SEDs of 126 blazars and identifying a remarkable 
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continuity in their sample and drew the following conclusions: 

• The first peak occurs in different frequency ranges for different samples/luminosity 
classes, with most luminous sources peaking at lower frequency (i.e. synchrotron peak 
is anti-correlated with the source luminosity). 

• The X-ray spectrum becomes harder while the 7-ray spectrum softens with increasing 
luminosity, indicating that Compton peak also moves to lower frequencies for most 
luminous sources. 

• The peak frequency of the 7-ray component correlates with the peak frequency of the 
lower energy one (the smaller the i^peak,sync, smaller the peak frequency of the high 
energy component). 

• The luminosity ratio between the high and low frequency components increases with 
bolometric luminosity. 

In this blazar sequence a single parameter, related to luminosity, seems to govern the physical 
properties and radiation mechanisms in the relativistic jets present in the BL Lac objects as 
well as in FSRQs II153II . In the simplest form of blazar sequence (Fig. |2.4[ ), one distinguishes 



low frequency-peaked BL Lac objects (LBLs) and high frequency-peaked BL Lac objects 
(HBLs). LBLs are more luminous and have their first peak in the infrared to optical bands 
and the second one is mostly in MeV region. The energy emitted in the second peak domi- 
nates the energy released in the first one. In case of HBLs, the first peak occurs in the energy 
range of UV to X-ray energies and dominates over the second peak in GeV — TeV energies. 
A simple one-zone leptonic model HI 621 describes the results satisfactorily by attributing 
the observed differences between these objects to the maximum energy of the accelerated 
electrons. 

The model given in HI 6211 was used to predict sources with detectable VHE emission 
by II123L Indeed, after this publication, ten new extragalactic TeV 7-ray emitters were dis- 
covered by 2007 and almost all of them fuUfiUed the criteria proposed in [I123L that observe 
the sources at moderate redshift of highest combination between the radio flux (at 5 GHz) 
and X-ray flux (at 1 KeV). This supports the major predictions of the blazar sequence, 
however there is no proof that there are no VHE 7-ray emitters outside this selection cri- 
teria I1129II . It was recently shown that the anti-correlation between the luminosity and the 
synchrotron peak frequency does not hold and it was an artifact of a selection effect [12421 . 
More recently few intermediate frequency peaked BL Lac objects (IBLs) have been found 



(see Table 2.1 ) bridging the gap between LBLs and HBLs. 
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Fig. 2.5: Sky map of total 42 extragalactic VHE 7-ray sources detected till date (August 2010). 
Redshift of sources indicated in different colors. The image has been adapted from |[T4]| . 



2.3.2 Blazars detected in TeV energy 



All extragalactic sources detected till August 2010 are shown in Fig. 2.5 and also listed in 



Table 2.1 along with their discovery references. Most of these detected sources are blazars 
except two radio galaxies (M87 and Centaurus A), two starbust galaxies (NGC 253 and 
M82) and three unidentified sources (MAGIC J0223+430, RGB 0648+152 and IC 310). 



Majority of the detected blazars as is clear from the Table 2.1 belong to HBL type, while as 
to date only three sources each of LBL and FSRQ have been detected. Four recently detected 
sources belong to type IBL. A giant elliptical radio galaxy M87 (jet is viewed side on) is 
in our immediate cosmic vicinity, just 55 million light years away (z=0.0044). While the 
Whipple group has reported only an upper limit on the TeV 7-ray flux from this source [I214L 
the HEGRA collaboration found a weak TeV 7-ray signal [|47]| . which was confirmed by the 
HESS [|5l, VERITAS [fTHl and MAGIC [El observations separately. M87 was also jointly 
observed by three ground-based VHE 7-ray telescopes VERITAS, MAGIC and HESS along 
with radio observations at 43 GHz during 2008. Simultaneous observations at the lowest and 
highest ends of the electromagnetic spectrum indicate that this AGN accelerates particles to 
very high energies in the immediate vicinity of its central black hole [31 J. This source is 
currently being observed in a coordinated effort by VERITAS and MAGIC and an increasing 
VHE 7-ray flux level has been measured reaching a historic high state of about 20% of the 
Crab Nebula flux during the observations on 09 April, 2010 UJ. 
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Table 2.1 : Extragalactic TeV sources detected so far (August 2010). 



Sr No 


Source 


Type 


Redshift 

z 


Discovery reference 


1 


Mkn421 


HBL 


0.030 


Punch et al., Nature 358 (1992) 477 


2 


Mkn 501 


HBL 


0.034 


Quinn et al., ApJL 456 (1996) L83 


3 


lES 2344+514 


HBL 


0.044 


Catanese et al., ApJ 501 (1998) 616 


4 


lES 1959+650 


HBL 


0.047 


Nishiyama, Proc. 26th ICRC (1999) 3, 370 


5 


PKS 2155-304 


HBL 


0.116 


Chadwick et al., ApJ 513 (1999) 161 


6 


H 1426+428 


HBL 


0.129 


Horan et al., ApJ 571 (2002) 753 


7 


M87 


radio 


0.0044 


Aharonian et al., A&A 421 (2004) 529 


8 


PKS 2005-489 


HBL 


0.071 


Aharonian et al., A&A 436 (2005) L17 


9 


lES 1218+304 


HBL 


0.182 


Albert et al., ApJL 642 (2006) LI 19 


10 


H 2356-309 


HBL 


0.165 


Aharonian et al.. Nature 440 (2006) 1018 


11 


IBS 1101-232 


HBL 


0.186 


Aharonian et al.. Nature 440 (2006) 1018 


12 


PG 1553+113 


HBL 


0.09 


Aharonian et al., A&A 448 (2006) LI 9; 
Albert et al., ApJL 654 (2007) LI 19 


13 


Mkn 180 


HBL 


0.045 


Albert et al., ApJL 648 (2006) L105 


14 


PKS 0548-322 


HBL 


0.069 


Superina et al., 30th ICRC, Mrida, 2007 


15 


BL Lacertae 


LBL 


0.069 


Albert et al., ApJL 666 (2007) L17 


16 


lES 0229+200 


HBL 


0.140 


Aharonian et al., A&A 475 (2007) L9 


17 


lES 0347-121 


HBL 


0.185 


Aharonian et al., A&A 473 (2007) L25 


18 


IBS 1011+496 


HBL 


0.212 


Albert et al., ApJL 667 (2007) L21 


19 


3C 279 


FSRQ 


0.536 


Albert et al.. Science 320 (2008) 1752 


20 


RGB J0152+017 


HBL 


0.080 


Aharonian et al., A&A 481 (2008) L103 


21 


lES 0806+524 


HBL 


0.138 


Acciari et al., ApJ 690 (2009) L126 


22 


W Comae 


IBL 


0.102 


Acciari et al., ApJ 685 (2008) L73 


23 


S5 0716+71 


LBL 


0.31 


Teshima et al., ATel 1500, 1 


24 


3C 66A 


IBL 


0.444 


Acciari et al., ApJ 693 (2009) L104 


25 


MAGIC J0223+430 


UNID 




Ahu et al., ApJ 692 (2009) L29 


26 


Centaurus A 


radio 




Aharonian et al., ApJ 695 (2009) L40 


27 


RGB J07 10+591 


HBL 


0.125 


Ong etal., ATel 1941, 1 


28 


PKS 1424+240 


IBL 




Ong et al., ATel 2084, 1 


29 


NGC 253 


starburst 


0.0008 


Acero et al.. Science 326 (2009) 1080 


30 


M82 


starburst 


0.0007 


Karlsson, N., arXiv:09 12.3807 


31 


VER J0521+211 


others 




Ong et al., ATel 2260, 1 


32 


RBS 0413 


HBL 


0.19 


Ong et al., ATel 2272, 1 


33 


lES 0414+009 


HBL 


0.287 


Hofmann et al., ATel 2293, 1 


34 


lES 0502+675 


HBL 


0.341 


Ong etal., ATel 2301, 1 


35 


PKS 0447-439 


HBL 


0.2 


Raue et al., ATel 2350, 1 


36 


PKS 1510-089 


FSRQ 


0.36 


Wagner, HEAD 2010 


37 


RGB 0648+152 


UNID 




Ong et al., ATel 2486, 1 


38 


IC310 


UNID 


0.019 


Neronov et al., arXiv: 1003.4615, 
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39 


4C +21.35 


FSRQ 


0.432 


Mariotti et al., ATel 2684, 1 


40 


AP Lib 


LBL 


0.049 


Hofmann et al., ATel 2743, 1 


41 


MAGIC J200 1+435 


HBL 




Mariotti et al., ATel 2753, 1 


42 


lES 1440+122 


IBL 


0.162 


Ong et al., ATel 2786, 1 



The MAGIC collaboration has discovered VHE 7-ray emission from a new unidentified 
source MAGIC J0317+413 above 100 Ge^ in March 2010 IS. The emission position is 
consistent with the head-tail radio galaxy IC310 (z=0.0189) located in the outer region of 
the Perseus cluster of galaxies. Their preliminary analysis indicates emission at the level 
of ~ 2.5% of the Crab Nebula flux above 300 Gel^. Further, the source was also recently 
detected by the Fermi LAT satellite based detector at a statistical significance of ~ 6a above 
30 GeV^l^. 

The VERITAS collaboration has discovered VHE 7-ray emission from four IBLs (a new 
class of AGN) W Comae (z=0.102) [|2l, PKS 1424 +240 (unkown distance) [3], 3C 66A 
(z = 0.444) lEa and 1ES1440+222 (z=0.162, last source in Table [O]) above 100 GeV. 
The detection of PKS 1424 +240 was also confirmed by the MAGIC group S and Fermi- 
LAT ll24l in VHE and HE respectively. It may be noted that the radio galaxy 3C 66B lies in 
the same field of view as 3C 66A at a separation of 0.12° and is also a plausible source of 
VHE radiation II278L The MAGIC collaboration has detected VHE 7-ray emission from 3C 
66A/B region [|78| favoring 3C 66B as the source. Advanced telescopes with better angular 
resolution will be able to resolve these sources clearly. In addition, the MAGIC telescope 
has also detected VHE 7-ray emission from another interesting object 3C279 which falls in 
the FSRQ class of AGN and is the most distant (z = 0.536) source detected so far fTM . Such 
detections indeed constrains current theories about the density of the EBL implying that the 
Universe is possibly more transparent at cosmological distances than eaher believed to be 
in the 7-ray band. Further, the addition of new blazar classes namely IBL, HBL and FSRQ 
to the VHE catalog is expected to play an important role in our understanding of blazar 
populations, their kinematics and VHE 7-ray production mechanisms. 



2.4 Production of 7-rays in AGN 

Although the SMBH at the center of the AGN, is the putative source of energy for the entire 
system, it is not intuitively obvious how energy can be extracted from the black hole. Par- 
ticularly it is not obvious why it should emerge in the form of a relativistic jet. Black holes 
are generally conceived as energy sinks -a gravitational sink hole that pulls in every thing 
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within the event horizon. While the formation of the hot accretion disk around the black hole 
is understandable in terms of conservation of angular momentum and hence, the reason that 
these objects are detectable as bright X-ray sources, it is less clear why beams of relativistic 
particles should emerge from the vicinity of the black hole 029511 . The challenges for under- 
standing the relativistic jets that are the heart of blazars are: what is the source of energy; 
how these relativistic energies can be achieved; what is the reason that two jets are formed 
pointing in exactly opposite directions and what is the mechanism by which the jets remain 
coUimated as they emerge through the turbulent surroundings of the black hole [I295II . 

Though there is no general consensus on the origin of the emission components seen in 
7-ray blazars, it is generally agreed that the low energy component arises from incoherent 
synchrotron emission by relativistic electrons within the jets. This is supported most strongly 
by the high-level variable polarization observed in these objects at radio and optical wave- 
lengths. The observation of HE and VHE 7-rays from blazars provides independent evidence 
that the production of radiation takes place in the jets and also the jets make a small angle 
with the line of sight. This conclusion is also drawn from the inferred compact emission 
region (based on short time variability): given the observed luminosity (Eq. 2.1), it would 
not be possible for 7-rays to escape 7 — 7 pair production unless there is relativistic beaming. 

It is generally agreed that the relativistic jets are caused by bulk relativistic motion, char- 
acterized by bulk Lorenz factor F, defined as 

r = (1 - /32)-i/2 (2.5) 

The Doppler factor of an object moving at (3 = v/c, making an angle 9 with the line of sight 
is 

(5 = l/(r(l - /3 cos ^)) (2.6) 

The introduction of Doppler factor alleviates the explanation of observed properties of blazars 
in a variety of ways [I102II . Some of these are explained by equation [2 . 2 1 and [23| In addition 



^source- 



the problem of accelerating particles to such high energies is alleviated by Vobs = S Vg. 
In practice, generally 5 ~ 10 — 20, so that the enhancement in energy can be as much as 
103-4 [|295J . 

In the following sub- section, we briefly summarize the models generally used to explain 
7-ray emission from blazars. The models are usually aimed not only to explain the obser- 
vational data but also to constain the jet parameters like the magnetic field, geometry of the 
emission region and the Lorenz factor of the jet. Since the responsible particles for VHE 
7-ray emission are for example accelerated electrons and protons, accordingly the models 
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are divided into lep tonic and hadronic categories. 
2.4.1 Leptonic models 

The charateristic double peak shape of SED of blazars is immediately suggestive of the 
synchrotron Compton model at work in these sources. Electrons are accelerated beyond 
the velocity of the bulk Lorenz outflow and the probable mechanism of acceleration are the 
shocks propagating down the jet which are caused by colliding inhomogeneities in the jet 
i.e blobs of material moving down the jet with different velocities. The electrons radiate 
synchrotron radiation in the magnetic field associated with the jet and produce the first peak 
in the SED. The position of synchrotron peak is determined by the efficiency of the shock 
acceleration and the cooling processes (i.e. synchrotron and Compton scattering) [I295II . The 
VHE 7-ray emission is explained by IC scattering of optical to X-ray photons off relativistic 
electrons or positrons in the jet. The following two sub-classes differ in the origin of seed 
photons for IC. 

Synchrotron self-Compton (SSC) 

In SSC models [I105l . ll63l . l227l . l268l . ll94ll the seed photons are synchrotron photons emitted 
by the same electron population. Experimental evidence for the SSC mechanism has been 
provided by the observed correlation of the X-ray and VHE 7-ray flux levels during large 
flares of VHE 7-ray emitting blazars [|228l . I276l. I207l. The SSC model is widely accepted 
to describe VHE 7-ray emission from the HBL objects. 

External radiation Compton (ERC) 

The second class of models involve external ERC [l269] . [T36l . [T0T] . [T37] . [T6T] . l292ll . wherein 
the seed photons are ambient infrared or optical photons, photons of CMBR or thermal 
radiation either directly from the accretion disk or first scattered by the surrounding gas and 
dust clouds. Hence in these models the 7-ray emission takes place outside the jet. The 
existence of such sources of soft radiation is deduced from observations of non-blazar AGN, 
i.e. in which the jet is not seen edge-on. The SED of these sources shows additional thermal 
bumps in UV, infrared and X-ray energies [12951 . 

One characteristic feature of BL Lac objects, i.e. lack of strong emission lines, is com- 
monly interpreted as an evidence that ambient photon fields are not as important as syn- 
chrotron photons. Therefore SSC models are more likely to explain the data than ERC mod- 
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els. However, often the predicted emission is below the observed emission ( e.g. in FSRQs), 
making ERC contribution necessary [I290II . 

Presently, the leptonic models independent of the specific assumption concerning the 
acceleration and radiation scenarios, represent the preferred concept for TeV blazars. Two 
important features characterize these models [|37l : 

• Electrons can be readily accelerated to TeV energies e.g. through the shock accelera- 
tion mechanism, 

• They can radiate X-rays and TeV 7-rays with very high efficiency via synchrotron and 
IC channels. 

2.4.2 Hadronic models 

Although the associated acceleration of protons is expected with at least the same efficiency 
as that of electrons (for most acceleration mechanisms), the hadronic models require proton 
acceleration up to energies of 10^° eV, otherwise they can not offer efficient 7-ray production 
mechanisms in the jets [[371 . These models assume that the observed 7-ray emission is 
initiated by accelerated protons interacting with; 

• Ambient matter: the so-called mass-loaded hadronic models [[252L 

• Photon fields: called photo-pion hadronic models [12241 ■ [225[ . [2261 . 

• Magnetic fields: "pure" proton synchrotron model [[36l[ . 

• Both magnetic and photon fields [[235[ . l236l[ . 

While the synchrotron and photo-pion hadronic jet models require acceleration of protons to 
extremely high energies, the mass-loaded hadronic models need protons of relatively modest 
energies. 

Hadron models are attractive as they are characterized by particle acceleration up to 
energies lO'^^eV, and explain CR acceleration to extreme energies. The leptonic models 
generally lack this virtue. The hadron models have no problem explaining the highest 7- 
ray energies but they do have problems accounting for the rapid cooling necessary to account 
for the short time variations observed. Additionally, they have difficulties in explaining the 
observed correlations of X-ray and VHE 7-ray emissions of blazars during flares. This is 
because proton cooling by synchrotron radiation is less than that by electrons by the factor 
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{nie/mpY. However cooling can also come from collisions with photons or ions [I295L 
These models also require extreme parameters of the acceleration region, particularly high 
plasma densities and high magnetic fields B ^ IOC 

A by-product of proton models is the production of energetic neutrinos (from decay of 
charged pions) which might be detectable with next generation of neutrino telescopes. The 
detection of such flux would effectively eliminate the lepton models. This possibility is also 
quoted as the justification for the construction of large neutrino telescopes. Orphan VHE 
7-ray flares i.e. unaccompanied by outbursts in other wavelengths e.g. observed for BL 
Lac object lES 1959+650 in 2002 [I128II can be interpreted as an indication of hadronic 
acceleration at work. 

The HE and VHE 7-ray observations strain both the leptonic and hadronic models but do 
not, at present, rule either out. However, the observations generally seem to favor leptonic 
models [|295ll . Note that the leptonic and hadronic acceleration may both be present in blazar 
jets at the same time e.g. hadronic acceleration (or a leptonic/hadronic admixture) might ac- 
count for a baseline flux and fast flare could be governed by additional electron acceleration 
( 1129011 and references therein). 

2.5 SSC models: Parameters and observables 

SSC models describe the multi- wavelength spectra of HBLs quite successfully, wherein 
synchrotron and IC emissions are responsible for low and high energy components respec- 
tively [|2761.I280H. 

In its simplified one-zone homogeneous version, SSC model assumes a spherical accel- 
eration region (blob of radius R), filled with an isotropic electron population and a random 
magnetic field (B), moving at relativistic speed (3 = v/c along the jet within a small angle 
to the line of sight. The observed radiation is strongly affected by the relativistic Doppler 
beaming factor (5). 

Considering the diffusive shock as the acceleration mechanism, the electron injection 
spectrum can be taken as a power law dN/dE oc i?^"^ with n\ c:^ 2. The spectrum is mod- 
eled to show a break point where radiative cooling through synchrotron radiation becomes 
the dominating process. At this energy it steepens to ri2 ~ 3 and thus can be described by 
three parameters, Emm, E^ax and Ehreak (left plot of Fig. |2.6[ ). as observed in BL Lacs, 
such an electron spectrum is able to produce the synchrotron photon spectrum in randomly 
oriented magnetic fields of approximately constant strength [11961 . Usually nine parameters 
are required to completely specify the SSC models, namely size of emission region (R), 
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Fig. 2.6: Schematic of the relativistic electron energy spectrum (left) and corresponding SED pro- 
duced by this electron injection (right). For high Lorentz factors, radiative cooling decreases 
the number density of electrons and leads to a break in the electron spectrum. On the right 
side, the breaks stem from radio self absorption of created synchrotron photon, radiative 
electron cooling and the cutoff in the electron injection spectrum. This figure has been 
adapted from 12901 . 

minimum electron energy (Emm), maximum electron energy (Emax), break energy (Ebreak), 
spectral index of electron spectrum before break (nl), index after break (n2), normalization 
factor of electron density (K), magnetic field (B) and the Doppler beaming factor (5). 

The measured SED can be used to drive seven "observable" quantities i.e. photon index 
of synchrotron radiation before break (al), after break (a2), peak synchrotron frequency 
(Usfi), synchrotron frequency at high energy cutoff (Us^max), peak IC frequency (z/c), total 
measured energy flux for synchrotron component (Is) and for the IC component (/c) [I176L 

Using the standard formula for radiation one can obtain the relation between the model 
parameters and observables. Observed synchrotron frequencies from a single electron of 
Lorenz factor Ebreak and 'jmax are respectively given by I1196II 

z/,,, = 1.2 X 10^5 7L,fc (2.7) 
iy„ = 1.2 X 10^5 7^,, -— (2.8) 

1 ~r Z 

The maximum value of the observed SSC energy in the Klein-Nishina regime is 



hl^c C\ '^max C 

1 + Z 



(2.9) 
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The Constant Ci < 1 represents the uncertainty of the Klein-Nishina effect. 

One more relation can be obtained from the ratio of synchrotron to Compton luminosities 
( 1117611 and references therein). 

where di is the luminosity distance and mb(= -B^/Stt) is the energy density of magnetic 
field. One more observable quantity is the variability time scale which is related to the size 
of emission region (_R < ctyarS). 

When the synchrotron photons are upscattered to VHE energies via the IC process by 



the electron population, the second peak in the SED of blazars is produced (Fig. 2.6). The 
position of the IC peak is determined by Ebreak, as the peak synchrotron power is emitted by 
the electrons with this energy. An approximate mapping between synchrotron and IC photon 
energies produced by electrons of given energy is represented as follows [I209II 



5/10 



Ejc/TeV^^-^^j {EjkeVf' (2.11) 

within the present uncertainties, the left side of both the peaks can be described by the same 
spectral index [I279II . 

The Doppler boosting factor is an important parameter which describes the effect of time 
dilation, boosting and jet focusing. Note that the observed source frequency appears as u.S, 
the emission solid angle is reduced by 5'^, adding to an increase in brightness by ~ 5^. An 
additional constant to the model is given by the requirement of source transparency which 
can be translated into a hmit on the optical depth of the source and thus on energy density 
of the soft radiation ( [I290II and references therein). Time dependent modeling has to be used 
whenever a flare duration is very short. To fit the observed data consistently from several 
Mrk501 flares, time dependent modeling and the introduction of the second emission zone 
is necessary [I208II . 



2.6 Implications of VHE 7-ray observations 

• VHE 7-ray observations of blazars have helped in resolving the nature of differences 
between LBLs and HBLs. It had been proposed that these objects are same except that 
the jets of LBLs are aligned more closely with our line of sight, based on their smaller 
numbers and higher luminosities. However, the observation of rapid variability and 
spectra that extends upto TeV energies in HBLs point to the differences between these 
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two sub-classes being more fundamental: the HBLs have higher maximum electron 
energies and lower intrinsic luminosities. 

• Simultaneous measurements of the synchrotron and VHE 7-ray spectra constrain the 
magnetic field (B) and Doppler factor (5) of the jet. In case of Mrk421, the observed 
correlation between VHE and optical/UV photons indicates that if both sets of photons 
are produced in the same region of the jet, 5 > 5 is required for VHE photons to 
escape pair-production losses significantly. If SSC mechanism produces the VHE 7- 
rays , 5 = 15 - 40 and = 0.03 - 0.09G for Mrk 421 and 6 ^ 1.5 - 20 and 
B = 0.08 — 0.2G for MrkSOl [I295II . Proton models which utilize synchrotron cooling 
as the primary means for proton energy losses require -B = 30 — 90G for 6 ~ 10, in 
order to match the timescales of the correlated emission ( [I295| and references therein). 
The Mrk421 values of 5 and B are extreme for blazars but still within allowable range 
and are consistent with its extreme variability. 

• VHE observations have constrained the types of models that are likely to produce 
the 7-ray emission e.g. the X-ray and VHE flares correlation is consistent with SSC 
models. The relative absence of flaring at EGRET energies can be explained because 
the lower energy electrons (responsible for 7-ray emission at EGRET range) radiate 
away their energy more slowly. 

• In ERC models, where the target seed photons are the ambient photons, should have 
energies < O.leV (i.e infrared band) in order to avoid significant attenuation of VHE 
photons by pair production. 

• Models in which 7-ray emission is due to proton proginators through cascades 
originating close to the base of the AGN jet, have difficulty explaining TeV emission 
observed in Mrk421. The high densities of unbeamed photons near the core, cause 
high pair opacity to TeV photons. These models predict that the radius at which the 
optical depth for 7-7 production drops below unity, increases with 7-ray energy. Hence 
VHE 7-ray s should vary either later or more slowly than MeV — GeV 7-rays HIOIL 
However, this contradicts the Mrk421 observations [I295II . 

• Both Mrk421 and MrkSOl are often considered to be very similar, but their 7-ray ob- 
servations have shown them to be quit different. Mrk421 always has a higher mean 
flux, shorter time variation and peak in the synchrotron spectrum {Usync) always in soft 
X-ray (does not show dramatic shift) as compared to MrkSOl [I295II . 
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2.7 Flux variability 

Most of the blazars are found to be highly variable at all timescales (from minutes to years) 
with change of flux level by more than an order of magnitude. The first clear detection of 
flaring activity in the VHE emission of an AGN came in 1994 observations of Mrk421 by 
Whipple telescope where a 10-fold increase in the flux, from an average level that year of 
15% to approximately 150% of the Crab flux, was observed. Many flaring states detected 
in VHE 7-ray energies are usually correlated with X-ray and some times also in the optical 
band [I113II . The correlation between X-ray and 7-ray flux has been seen in Mrk421, MrkSOl 
and PKS 2155-304. However, a campaign on the TeV blazar 1ES1959-I-650 has shown an 
orphan 7-ray flare without X-ray or optical counterpart. An orphan flare for Mrk421 was 
also found in X-ray without a TeV counterpart ( 1112911 and references therein). Spectral 
shape variability has been established for the strong TeV blazars like Mrk421, Mrk501, 
H1426-I-428 [SI and lES 1959-1-650 [l46],|66l. However during a very strong outburst of 
PKS 2155-304 no spectral change was observed [|6TI . 

The rapid variability found in TeV blazars provides additional constraints on emission 
models. The flux doubling time for Mrk421 was found to be ~ 30 minutes [I156L while 
for Mrk501 and PKS 2155-304 the variability timescales were shorter than 3 minutes fTO', 
|6TI . Very large amplitude flux variability at short timescales implies that TeV emission 
originates from a small region very close to SMBH and therefore VHE 7-ray observations 
provide an unique opportunity to study the close vicinity of the central object [I129II . Several 
mechanisms have been proposed to explain the observed fast flare [12901 : 

• Additional population of more energetic electrons or altering the electron injec- 
tion spectra: Large variations in the flux are achieved in X-ray and 7-ray energies 
if Imax is increased by a factor ~ 5, with no noticable change in flux levels at other 
energies. Alternatively flare can also be produced by increasing the luminosity of the 
injected electrons. 

• Modulation in the seed photon field for the IC process: Assuming a small X-ray 
hot spot rotates on the surface of the accretion disk whose emission may be coUimated 
by the hot spot geometry by an ordered magnetic field. Although strong absorption 
does not allow generation of VHE 7-rays in this region, the energy release of such a 
region can explain fast flaring e.g. as observed in Mrk421. 

• Sub-shocks inside the jet e.g. unsteady jets, fed by discrete ejecta: Fast flares can 
also be explained by assuming unsteady jets, fed by discrete ejecta. For example M87 
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(often interpreted as misaligned blazar) shows bright knots in the jets in optical region 
which can be interpreted as internal shocks. 

• Change in the acceleration environment: Fast flaring implies small acceleration 
regions or high Doppler factors. Considering small blobs with enhanced magnetic 
field strength inside the acceleration region can explain the flares. 

The observed correlation of X-ray and VHE 7-ray emissions can help to discern the models, 
e.g. the correlation between these two bands points to SSC origin and hence to sub-shock 
scenarios while orphan VHE flares favor ERC and hence modulation models. 

In this thesis we study the blazars Mrk421 and HI 426+428 in VHE energy range above 
1 TeV. The observations were carried out during 2005-2007 and 2004-2007 on Mrk421 and 
HI 426+428 respectively. We present the results of our data analysis for these sources in 
chapter 6 and 7 of this thesis. 



Chapter 3 



Detection technique of VHE gamma-rays 



3.1 Introduction 

Since 7-rays encompass a very wide range of energies (1 MeV to 10^° eV), we can not 
expect that a single type of detector will work at all 7-ray energies. As already mentioned 
in chapter 1 that depending on the energy regions the satellite experiments, atmospheric 
Cherenkov telescopes, air shower arrays and fluorescence detectors are employed to detect 
7-rays . At MeV — GeV energies, 7-rays have been typically observed with space-based 
instruments but at VHEs these instruments are completely unusable because of the rapidly 
falhng 7-ray fluxes. Therefore astronomy at energies more than few tens of GeV, can only 
be done by ground-based instruments using the Earth's atmosphere as the detection medium. 
With the advent of lACTs in late 1980's, ground-based observations of TeV 7-rays became 
possible and since then the number of 7-ray sources has rapidly grown up to over a hundred. 
It is worth mentioning here, the recently launched Fermi LAT covers the energy range from 
about 20 MeV to more than 300 GeV. 

For most wavelengths, the atmosphere poses significant challenges for astronomy. The 
quality of optical and infrared astronomy is compromised by observations made through the 
atmosphere. For UV, X-rays and 7-rays , direct detection by ground-based instruments is 
impossible. A cascade of secondary particles called "EAS" is initiated by VHE 7-rays , 
which propagate down through the atmosphere. The secondary charged particles and pho- 
tons produced in the EAS can be detected by the ground- level instruments. Thus the possible 
disadvantage of the atmospheric absorption is converted into an advantage by the ground- 
based 7-ray telescopes [I239II . An lACT detects the ACR emitted in the atmosphere by sec- 
ondary electrons which are produced in the cascade initiated by a primary 7-ray photon. 
CRs (mostly protons, nuclei, electrons) also produce the ACR via hadronic interactions and 
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Fig . 3 . 1 : Variation of atmospheric depth or air mass grammage (g cm~'^) with the atmospheric height 
from sea level. Corresponding variation in number of radiation lengths is also shown. 
Mount Abu (1300 m asl) corresponds to atmospheric depth of ~ 890 5 cm~^ and radia- 
tion length of ~ 24.2. 

essentially become background for observations of 7-rays . These two types of primaries 
differing in development of EASs, can be distinguished by lACTs. 

In this chapter, we briefly describe the model of the atmosphere which is an integral part 
of detectors at VHEs, development of an EAS in the atmosphere, production of ACR and the 
imaging ACT which has been proven to be the most sensitive technique for the detection of 
VHE 7-rays . 



The Earth's atmosphere is a layer of gases surrounding the planet that is retained by its 
gravity. Dry air contains roughly (by volume) 78.09% nitrogen, 20.95% oxygen, 0.93% 
argon, 0.038% carbon dioxide, and small amounts of other gases. Air also contains a variable 
amount of water vapor, on an average around 1%. The atmosphere does not technically end 

at any given height, but becomes progressively thinner with altitude. 

It is a common practice to remove the effect of the density of the medium in which an 
EAS develops by defining the pathlengths in units of g cm~'^. The atmospheric depth X [in 



3.2 The atmosphere 
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g cm ^] is related to the density profile of the atmosphere by ( 1117611 and references therein) 

X= p{h')dh' (3.1) 



Jh 

where p{h) is the density of the atmosphere at altitude h. 

For an isothermal atmosphere, the above equation can be described by means of a constant 
called scale height ho as 

X = Xgexp{-h/ho) (3.2) 

where Xg ^ 1030 (7 cm~^ and ^ 8.5 km. However, for a more realistic model of the 
atmosphere the relationship between X and h I1259II is given by 

The change of atmospheric density with altitude is given by [|63l 

p = Poexp{-h/ho) (3.4) 

Where po ~ 1.205 x 10"^ g cm~^ is the typical atmospheric density at sea level. 

The index of refractiorj^of the atmosphere n is a critical parameter for understanding the 
production of ACR. If we define r] = n — 1 then r] is proportional to the density and is given 
by 

r] = r]oexp{-h/ho) (3.5) 

where r]o = 2.9 x 10""* and ho = 7.25 km, differing slightly from the standard scale height 
of the atmosphere which is 8.5 km. 

The atmosphere is far from an ideal detector. Poisson fluctuations in the Night Sky 
Background (NSB) light limit the energy threshold of an lACT. Atmospheric absorption and 
scattering significantly modifies the spectrum of the detected ACR, especially in the UV part. 
The dominant contributions to absorption or scattering of ACR are [63 ]: 

(i) Rayleigh scattering: where photons get scattered by particles of much smaller size than 
the wavelength of the photon. The absorption length (A/j) depends on the fourth power 
of photon wavelength i.e. oc A'^. 



?/ = v^e/i where e is dielectric constant and /i is the magnetic permeability of the medium. 
For air /i = 1 and refractive index rj = 1.00029 



Chapter 3. Detection technique of VHE gamma-rays 



52 



(ii) Mie scattering: where photons are scattered by particles with dimensions comparable 
to the wavelength (e.g. aerosoles), here Am oc A. 

(iii) Absorption by ozone: O3 + photon — )■ O2 + O. 

This strongly absorbs near UV region ( 200 to 300 nm). 

(iv) Scattering by water vapor in clouds or from local humidity. 

A rough estimate for the magnitude of these effects under good weather conditions can be 
obtained by considering the empirical model for the absorption length derived by Hillas II181II 

A = 1385 - 0.3798 a: - 0.000556 {gcm-^). (3.6) 

Radiation length 

HE electrons predominantly lose energy in matter by bremsstrahlung while high-energy pho- 
tons by pair production. The characteristic amount of matter traversed for these related in- 
teractions is called the radiation length Xq. It is both the mean distance over which a HE 
electron loses all but 1/e of its energy by bremsstrahlung, and 7/9 of the mean free path for 
pair production by a high-energy photon. For example an electron with initial energy Eq 
travels a distance X (g cm^^) in a dissipative medium, its energy can be written as 

Ex = Eoexpi-X/Xo). (3.7) 

The radiation length is given, to good approximation, by the expression ( 1128111 and references 
therein) 

Xoigcm-') = ( ^-^^ —] (3.8) 

\Z{Z + l)ln{287/^/Z)J 

where Z and A are the effective atomic number and atomic mass of the medium. For air 
at STP, Z = 7.34 and A = 14.37: therefore using above formula, the radiation length for 
electron Xq ~ 36.S gcm~'^ in air. Another term called shower unit is defined as i? = 

Xq ln2. 

Note that total vertical thickness of the atmosphere is 1030 g cm^'^ up to sea level and 
since the radiation length is 36.8 (7 cm~^, this amounts to ~28 radiation lengths of atmo- 
sphere. The relationship between atmospheric height and the air mass grammage is shown 



in Fig. 3.1 For electrons at lower energies (below few tens of MeV?,), the energy loss by 
ionization is predominant. Therefore one more parameter called critical energy is defined 
which characterizes the effectiveness of a medium as a radiator. 
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Critical energy 

Critical energy (Ec) is the energy at which the rate of energy loss by ionization becomes 
comparable with the rate of energy loss by radiation (e.g. Bremsstrahlung in case of elec- 
tron). The Ec for electrons (positrons) in gases can be computed by the following expression 
( 1128111 and references therein): 

The electron has a critical energy MeV in the atmosphere (air). 



3.3 Extensive air shower 

When VHE 7-ray photons or other CR particles enter the atmosphere, they interact with 
nuclei of atmosphere molecules typically at an altitude of about 20 — 25 km (the point of 
first interaction) and generate the so-called EAS. The interaction of electrons, muons and 7- 
rays with nuclei is electromagnetic i.e. they generate secondary particles by pair production 
and Bremsstrahlung, while protons and ionized nuclei interact via hadronic interactions i.e. 
they produce pions, muons and kaons as secondary particles. As shown in Fig. X2]the sec- 



ondary particles further interact with atmospheric nuclei starting a cascade of particles i.e. 
an EAS. The number of secondary particles grow up rapidly, while their energy decreases 
continuously in each interaction. Since the energy of the primary particle is distributed over 
all the secondary particles, at some point of shower development the energy of the shower 
particles is too small to produce further particles i.e. it falls below Ec. At this point shower 
reaches to its maximum size called "shower maximum" beyond which particles suffer en- 
ergy loss due to ionization, as a result the number of shower particles decay exponentially 
and finally the shower dies out. Depending on the type of primary particle we deal with two 
types of air showers: 



3.3.1 Electromagnetic EAS 

The Earth's atmosphere is opaque to high energy 7-rays . At sufficiently high energies (more 
than twice the electron rest mass) the predominant interaction of 7-rays with atmospheric 
nuclei is pair production. For a 7-ray of energy > 10 GeV, the photon is effectively replaced 
by two charged particles (e^) sharing its energy and traveling in almost the same direction 
as that of the 7-ray photon. Basically two processes mainly contribute to the shower devel- 
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Fig. 3.2: Schematic description of the extensive air showers. Left sketch shows the development of 
an electromagnetic shower while right one is for the hadronic shower. This figure has been 



adapted from 1129011 . 



opment: pair production and Bremsstrahlung. Pair produced electrons and positrons further 
produce photons via Bremsstrahlung; such photons in turn produce pair, and so on, giv- 
ing rise to a shower of charged particles (mainly e^*") and photons called an electromagnetic 
EAS, as shown in the left of Fig.|3]2l The process is described in Il264l . l238ll . 



At shower maximum the energy is ^ 300 MeV. The energy loss of an electron due 
to Bremsstrahlung is proportional to its energy: 

- dE^/dx = Ee/Xo (3.10) 

Once the energy of reaches Ec = 86 MeV, shower gets absorbed. As already mentioned 
the radiation length of for Bremsstrahlung is 36.8 g cm~'^ and the interaction length or 
mean free path of photons for pair production is 9/7Xq i.e. 47.3 (^cm^^. Therefore the 
same order of interaction lengths in both the processes, indicates a quite simple and sym- 
metric structure of an electromagnetic shower with respect to the original 7-ray trajectory 
(also called shower axis). The average opening angle of emission in all these processes 
is 1/7, where 7 is the Lorenz factor of electron. The resulting electromagnetic cascade is 
remarkably tightly bunched along the shower axis. However, tansverse broadening (lateral 
development) of the shower takes place which is predominantly due to multiple coulomb 
scattering process of electrons traversing in the air. Although small, the lateral spread of 
secondary particles also occurs due to the Earth's magnetic field. Bremsstrahlung and pair 
production also contribute to some extent in the transverse development of the shower. The 
shower particles spread out laterally to tens or even hundreds of meters, moving downwards 



Chapter 3. Detection technique of VHE gamma-rays 



55 



Priiiiai'v 

Panicle 




R E. 



2R E / 2 



3R Eo/4 



4R E / 1 



MSML 



iiR E i> / 2 



Fig. 3.3: A simple model of an electromagnetic shower development proposed by Heitler. This figure 
has been adapted from [,189,1 . 



with relativistic speeds. 

The distribution of Coulomb scattering is described by the theory of Moliere [|89ll . It has 
a roughly Gaussian shape for small scattering angles, while at larger angles it approximates 
Rutherford scattering with longer tail than a Gaussian distribution. The lateral spread of 
electromagnetic showers in different materials scales quite accurately with Moliere radiu:^ 
(Rm), which in turn depends on radiation length and critical energy in a given medium. 
Accounting that the distance corresponding to Xq varies with density of atmosphere (equa- 
tion (|3.1[)) the expression for Rm is ( 111551 and references therein): 



9.6(9 cm-^) _ (3.11) 

Pair 

At sea level Rm is approximately 80 m. 

The basic features of the longitudinal development of an electromagnetic shower, illus- 
trated by a toy model in Fig. |3.3|has been proposed by Heitler [I178B . The model considers 



only pair production and Bremsstrahlung process with assumption that both radiation length 
for electron and interaction length for photon have same value i.e. Xq. Instead of Xq the 



^The Moliere radius is a characteristic constant of a material giving the scale of the transverse dimension 
of the fully contained electromagnetic showers initiated by an incident high energy electron or photon. By 
definition, it is the radius of a cylinder containing on average 90% of the shower's energy deposition. It is 
related to the radiation length Xq by the following approximate relation; 

Rm = 0.0265Xo(Z + 1.2) 

where Z is the atomic number. 
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quantity defined earlier called shower unit (R = Xq ln2) is used so that the probability of 
an electron or photon to interact along a distance R is therefore, exp(—R/Xo) = 1/2. It is 
also assumed that in each interaction the energy of the interacting particle is equally divided 
between the products of the interactions. 

As shown in Fig.|33, a 7-ray photon of initial energy (Eq) which enters the atmosphere, 



undergoes pair production after traveling a distance R. The resulting each has on an 
average half of the initial energy i.e. Eo/2. After another distance R, each particle produces 
a photon of an average energy Eq/A via Bremsstrahlung, while retaining Eq/A energy with 
itself. At a distance nR from the first interaction, there are 2" particles each with an average 
energy of Eo/2"'. The particle multiplication continues till the average energy of particles 
drops below critical energy (Ec = 86 MeV for air) and thus halts the cascade. 

While for 7-rays , the critical energy corresponds to a point where cross section for 
Compton scattering is approximately equal to pair production cross section {crcompton ~ 
Cpair), for electron it is the energy at which energy loss due to ionization is comparable to 
that by Bremsstrahlung {(dE/dX) radiation ~ (dE/dX) ionization)- At shower maximum total 
number of particles (e^ and photons) are 

N,na. = Eo/E, = 2\ (3.12) 

Assuming that the shower stops abruptly at shower maximum, the maximum penetration 
depth of the shower can be expressed as 

Xma. = nR= ^^^^^R = Xo ln{Eo/E,). (3.13) 

A more realistic approximation for the number of particles (e^) produced in an electromag- 
netic shower as a function of atmospheric depth t{= X/Xq) expressed in units of radiation 
lengths and primary energy Eq is given by Greisen HI 6911 

iVe(t,Eo) = ^=2^=exp[t{l-l.hlns)] (3.14) 
^ln{EQlEc) 

Where s is a dimensionless quantity called shower age and defined as 

?>t 

(3.15) 



t + 2ln{EQ/Ec) 

The parameter s, basically indicates the degree of development of the shower whose meaning 
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Atmospheric depth (radiation lengths) 

Fig. 3.4: Longitudinal development of an electromagnetic shower for different primary energies us- 
ing Greisen approximation. Different blue curves represent different energies (£'o) of pri- 
mary photons. The green lines characterize the shower development by the shower age s. 
s = 1 represents the maximum development of the shower. The altitude of the TACTIC 
telescope (1300 m asl) is also indicated. 



is given by the derivative dNe(t, Eo)/ds, 

= N^{t,E) X ^—^^ (3.16) 

As one can see from the above derivative ^ > in the range < s < 1, i.e. number 
of grows as s increases. At s = 1, ^ = 0, i.e. A^e becomes maximum, while in the 



range s > 1, ^jf < 0, shower dies out. Therefore the parameter related to the shower such 
that [12641 : 

• s = — )■ first interaction point, 

• 5 = 1—)- the maximum development of the shower, and 

• s = 2 — the point at which shower begins to extinguish. 



By using this definition of shower age, the depth of shower maximum can be computed 
in terms of number of radiation lengths as tmax = ln(Eo/Ec) i.e. logarithmic dependence 
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on initial energy, same as obtained with the simple model of Heitler. Fig. |3.4| shows the 



longitudinal development of the shower as described by equation (3.14), for different primary 
energies. The higher the energy of the primary, the deeper the air shower penetrates into the 
atmosphere. The shower maximum is at about 13 to 7 km above sea level for primary 7- 
ray energies between 50 GeV to 10 TeV. 

The lateral distribution of electrons is modeled by the NKG-formula (Nishimura-Kamata- 
Greisen) [11691 . 11921 . which describes the electron density as a function of the distance r from 
the shower axis 

, NJt,E) f r r Y'^-' 1(4.5 - s) 

where T is the Gamma function and Rm = 21.2 MeV.^Xo/ Ec) ~ 80 m in air at sea level. 
Nearly 99% of shower energy is contained within 3.5Rm- 

Air showers induced by primaries, pose an irreducible background for 7-ray showers, 
as they develop in the same way as 7-ray showers do. Due to the steep electron spectrum 
of index ~ 3.3, these showers can normally be ignored above energies 500 GeV ( [12901 and 
references therein). 

Roughly 10% of the energy of primary 7-rays reaches an observational level of 5200 m 
above sea level in the form of electromagnetic particles, while only 5% of primary particle's 
energy reaches to an altitude 4300 m. The fluctuations in the development of an EAS are 
due to fluctuations in the depth of the first interaction and the probability that the primary 
particle will survive to an atmospheric depth of radiation lengths is given by [63i 

P{N) = exp{-^N) (3.18) 

The shower fluctuations excluding the point of first interaction, are described by [I290B 

ANe{s) ^{s - 1 - 3 In s) Ne{s) (3.19) 

Due to their statistical nature, the fluctuations from shower to shower are quite large and 
have important consequences particularly for the energy estimation of the primary particles. 



3.3.2 Hadronic EAS 



When CR particles, most of which are protons enter the atmosphere, the inelastic scattering 
with atmospheric nuclei causes nuclear disintegration and production of mesons e.g pions 
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Table 3.1 : Mesons produced in hadronic EAS. 



particle 


rest mass (MeV/c^) 


mean life (s) 


decay mode 


probability 




135 


8.4 X IQ-i^ 


7 + 7 
7 + e~ + 


0.988 
0.012 


7r± 


139.6 


2.6 X IQ-*^ 


7r+ /i+ + Uf, 

7r+ — )■ e+ + Ue 
TT" + z/- 


0.999 
0.999 
0.0001 
0.0001 




493.7 


1.24 X IQ-s 


^ /i± + z/^ 
-> 7r± + vrO 


0.635 
0.212 




105.7 


2.2 X 10-6 


/i+ e+ + z/g + z/^- 
/i- e- + z/- + z/^ 





and kaons as well as nucleons. The secondary particles in turn themselves produce disin- 
tegrations thus starting a cascade of particles called hadronic EAS, as shown at the right of 



Fig. 3.2 This interaction can be represented by the following equation 



CR + AtmosphericN uclei{AN) =^ CR' + AN' + mvr^ + nvr^ + other mesons 

CR', which is a fragment of the primary CR can keep repeating the above interaction or 
even reach the ground if the original CR has sufficient energy. 

Nearly 90% of the secondary particles in a hadonic shower are pions and remaining 10% 
are kaons and antiprotons. One third of the total pions are the neutral pions(7r°). The decay 
modes of mesons along with their rest mass, mean life and corresponding probability of 



decay are summarized in Table 3.1 It is clear from the Table that 7r° has very short life 
and decays into two photons before it undergoes any hadonic interactions. Even a tt" with 
100 GeV energy will decay in the laboratory frame in 6 x 10~^^ s. Photons from 7r° decay 
initiate the electromagnetic EAS. 

In the early phase of shower development, high Lorenz factor of mesons results in 
hadronic interactions while at a later stage their decay becomes more important. On the other 
hand, due to their longer lifetime and smaller cross section most of the produced muons are 
able to retain their energy and reach the ground. Due to high energy, their life is further 
extended by a factor E^/m^c^, where and are energy and rest mass of the muon. The 
major energy loss of muon's energy is through ionization, while the decay of low energy 
muons can start the electromagnetic sub-shower. Muons and neutrinos are the main source 
of energy loss in hadronic EAS 015511 . 

The interaction length or mean free path (A) for a proton is more than two times larger 
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than the radiation length in an electromagnetic shower. The interaction lengths of protons 
and mesons in air dXlTeV energy are HI 5711 

\p = 83 g cm~'^, 

= 107 gcm~'^, 
Xk = 138 gcm-^. 

As a result, hadronic showers penetrate deeper into the atmosphere (i.e. depth of shower 
maximum is larger) as compared to a pure electromagnetic EAS of similar primary energy. 
A hadronic shower grows until the energy per nucleon is below the pion production threshold 
of about 1 GeV. A hadron-induced EAS therefore consists mainly of three components: 

(i) A hadronic core, containing strongly interacting high energy nucleons, 

(ii) Mesons constantly feeding electromagnetic sub-showers (mainly from 7r° decay), and 

(iii) Non- interacting muons and neutrinos, carrying most of the primary energy. 

Study of hadonic-induced EAS can be performed by detailed Monte Carlo simulations. How- 
ever incomplete understanding of hadronic interactions at VHE leads to rather large uncer- 
tainties in the results. 

3.3.3 Differences between photon and hadron initiated showers 

7-rays are observed in the unavoidable huge background (orders of magnitude more) created 
by hadronic CRs, making it necessary to search for background rejection methods. These 
methods rely on the characteristic differences between the two shower types, which are as 
follows: 

• As already mentioned, because of the smaller cross sections of interaction of pro- 
ton with air, its interaction length is much larger (~83 g cm^'^) compared to radiation 
length of 7-ray photon (~37 (7 cm^^). The longitudinal development of an EAS is gov- 
erned by the interaction length of primary particle (also on secondary particles), the 
hadronic shower develops deep into the atmosphere as compared to that of an electro- 
magnetic EAS. 

• The strong interactions are different from the electromagnetic ones and as an effect the 
mean transverse momentum carried by pions and kaons is higher than electrons and 
positrons. The hadronic showers therefore have more lateral extent. 
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Fig. 3.5: Monte-Carlo simulated extensive air showers of vertically incident ITeV 7-ray photon 
and proton. Courtesy F. Schmidt, "CORSIKA Shower Images" [HI. 
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Generally, hadronic showers have several sub-cascades, resulting in much more wider 
and irregular shower than electromagnetic showers. This feature is used to separate 
7-ray induced showers from hadronic ones. However if only the sub-cascade of the 
hadronic EAS is detected, it can mimic 7-ray induced EAS and becomes an irreducible 
background. 

In addition to these intrinsic differences arising during shower development, one char- 
acteristic difference comes from the nature of primary particle. The 7-ray showers 
point back to the source of origin of the primary 7-ray , while the arrival direction 
of proton-induced-showers is isotropized due to the deflection in galactic and extra- 
galactic magnetic fields. This feature is the back bone for the background rejection in 
lACTs particularly stand alone systems like TACTIC. 

In general, the hadron induced showers produce less ACR (discussed in the following 
section) than the 7-ray induced showers. 

The proton showers are subjected to larger fluctuations in the shower development, 
both within the shower and between individual showers. These fluctuations are caused 
by the larger interaction length, larger transverse momentum of secondary particles 
and variety of secondary particles. 

Photoproduction or photomeson is another process occurring in both electromagnetic 
and hadronic showers. Here 7 -1- nucleus — hadrons: Hadronic states, mostly pions are 
produced [I127L The cross section for this process increases with increasing photon 
energy, but has a small value of 1 — 2 mb for 0.1 — 20 TeV photons [EHJ. This is much 
smaller than the cross section for pair production which is 520 mb. By measuring the 
muon component, it is possible to distinguish between 7-ray and proton induced air 
showers. As these muons originate from decay of pions, their rate is at least an order of 
magnitude less in 7-ray showers. This procedure of 7 and hadron separation is limited 
by the smaller number of muons at lower energies i.e. < 10 TeV and by the increased 
occurrence of photoproduction processes at energy > 1 PeV ( H 12711 and references 
therein). 



Fig. 3.5 illustrates the simulated EAS for 7-ray and proton primaries of same energy 
(1 TeV). As an EAS develops through the atmosphere, the charged particles mainly electrons 
produce ACR by polarizing the medium and fluorescent radiations by exciting the nitrogen 
molecules of the atmosphere. The atmospheric Cherenkov radiation is the main topic of 
concern here. 



Chapter 3. Detection technique of VHE gamma-rays 



63 



3.4 Atmospheric Cherenkov Radiation 

When a charged particle moves through a medium with the relativistic speed which exceeds 
the local speed of light in that medium, the emitted radiation is called Cherenkov radiatioij^ 
This radiation was discovered by P. A. Cherenkov in 1934 [I117II . while the theoretical ex- 
planation for the origin of ACR was given by Frank and Tamm in 1937 [I277II . Considering 
the case of an electron moving in a transparent medium (i.e. zero absorption coefficient thus 
refractive index is real). As shown in Fig 3. 6(a) [ electron is moving slower than the speed of 



light. The atoms of the medium are represented by circles near the track AB of the electron. 
At a position P on the track, the electric field of the particle distorts the atoms such that 
positive charges of nuclei are displaced towards the moving electron. The medium is thus 
polarized and radiates electromagnetic pulses, as atoms return to their normal position once 
the particle moves away. Because of the complete symmetry of polarization, the resultant 
field at a large distance is zero hence no emission is observed. However, when the electron 
moves faster than the velocity of light, the symmetry in the plane perpendicular to the track 
of particle still exists but there is a resultant dipole field along the particle track (Fig. |3. 6(b)] ). 
The waves from each small track interfere constructively and the radiation is observed. When 
V < c/n, the electromagnetic waves travel ahead of the particle, setting instant dipoles in the 
medium resulting in no net polarization. In the other case, no instant dipoles are set ahead of 
particle and net polarization exists along the track of the particle [I189II . 



Fig. 3.6(c) illustrates the geometrical interpretation of the ACR. The emission is observed 
only at an angle with respect to the particle track and is called Cherenkov angle (9). This 
angle represents the direction in which waves from arbitrary points on a particle track are 
coherent and forming a plane wavefront. The coherence takes place when time taken by the 
particle to move from A to B is same as the time taken by the wave from A to C i.e. 

cos 9='-^ = -^ (3.20) 
pet pn 

where v = (3 cis the particle's velocity and n is the refractive index of the medium. It follows 



from equation (3.20) that for a given medium there is minimum velocity for a particle called 



threshold or critical velocity 

Pmin = - (3.21) 

n 



^In this thesis, we call it Atmospheric Cherenkov Radiation (ACR) as in the case of the VHE 7-ray detection 
the Earth's atmosphere acts as a detector wherein the Cherenkov radiation is produced whenever a VHE photon 
or proton impinges at the top of the atmosphere. 
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track of charged particle 




(a) Particle moves with (b) Particle moves with (c) Propagation of ACR in the 
low velocity (v < c/n) high velocity {v > c/n) medium 



Fig. 3.6: Illustration of the polarization produced in a dielectric medium by a charged particle as it 
passes through it. 



below which no radiation takes place. Therefore the minimum velocity required for emission 
of ACR with angular frequency u depends on the refractive index of the material for the 
frequency uj. At (3min , = ^° i.e. radiations are emitted along the particle track. On the 



other hand, if particle moves slower than the radiation {v < c/n), the equation ( 3.20 1 has no 
solution and hence no radiation takes place. Moreover, angle of emission is maximum when 
[5 = 1 i.e. particle is ultra relativistic, therefore 

Omax = COS~^ (3.22) 

The emitted radiation is in the visible and near visible wavelength band for which n(uj) > 1 . 
A real medium is always dispersive, so radiations are restricted to those frequencies for 
which n(u) > ^. For X-ray region n{uj) < 1, i.e. radiation is forbidden since equation 
(1122]) can not be satisfied [fT89ll . 



There are two more conditions, which need to be fulfilled to achieve coherence HI 8911 

(i) The length of the track of the particle in the medium should be larger than the wave- 
length of the radiation, otherwise diffraction would be dominant. 



(ii) The velocity of the particle should be constant during its passage through the medium 
i.e differences in the time for particle to travel successive distances A should be less 



Chapter 3. Detection technique of VHE gamma-rays 



65 



than the period (A/c) of the emitted radiation. 

The quantum theory of Cherenkov radiation includes the effects of emitted radiation on the 
motion of the particle. These effects are rather small as the energy of the radiated quanta 
is much smaller than the energy of the particle. Quantum treatment modifies the equation 
(|3.20|) to the following expression ( 1128 11 and references therein) 



1 , Ao (n^ - 1) v/(r^ 
cose = — + - — — (3.23) 



where Aq = 0.0024 nm is the Compton wavelength. 



3.4.1 Spectral distribution of Cherenkov radiation 

Considering a particle of charge ze, the number of Cherenkov photons emitted per unit path 
length and per unit wavelength interval are given by [|189l 

d?N 2W / 1 \ 

X 1 - T^^T-JTTT (3.24) 



dxdX A2 V /32n2(A) 

where a = /he = 1/137 is the structure constant. For electron moving along a track of 
length /, the number of emitted Cherenkov photon in the wavelength range Ai and A2 is then 
given by [fT90l 

Most of the photons are therefore produced in a short wavelength range i.e. in UV and blue 
waveband. No ACR is produced below 180 nm. Considering the frequency v of Cherenkov 
photons, the total energy radiated is W = N hu. The spectral distribution may be repre- 
sented in terms of energy per unit path length per unit frequency interval i.e. {(PW/dl dcu) oc 
u or energy per unit path length per unit wavelength interval i.e.{d'^W/dl dX) oc 1/A^ or 
number of photons per unit path length per unit frequency interval i.Q.{d^N / dl doj) oc const 
or number of photons per unit path length per unit wavelength interval i.e.{d'^N/dl d\) oc 



l/A^ [I190L The differential photon spectrum of Cherenkov photons is shown in Fig. 3.7 for 
photons emitted at 10 km and observed at 2 km altitude. The observed Cherenkov photon 
spectrum differs considerably from the produced spectrum due to the transmission charac- 



teristics of the atmosphere, explained in section 3.2 

The energy carried by ACR is given by the Frank-Tamm formula which yields the amount 
of ACR emitted as a charged particle moves through a medium at superluminal velocity. The 
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Fig. 3.7: Spectra of Cherenkov light emitted by vertical 7-ray showers of different primary energy 
at 10 km height is represented by solid curves while the corresponding detected spectra 
(effect of absorption by ozone and scattering) at observational level is represented by dashed 
curves. This figure has been adapted from II290I . 



energy emitted per unit length per unit frequency is [I277II 



dE /i(w)g' 



dxdoj 



An 



ui 1 



(3.26) 



where ^ is the permeability of the medium, q is the electric charge of the particle. ACR 
does not have characteristic spectral peaks, as typical for fluorescence or emission spectra. 
Higher frequencies are more intense in ACR. This is why visible ACR is observed to be 
brilliant blue. In fact, most ACR is in the ultraviolet spectrum. The total amount of energy 
radiated per unit length in air by an electron is given by 



dE 

dx 



e 

An 



v>c/n{ui) 



1 



duj. 



(3.27) 



This integral is done over the frequencies for which the particle's speed v is greater than 
speed of light in the media The integral is non-divergent because at high frequencies 
the refractive index becomes less than unity. 



3.4.2 Cherenkov radiation in EAS 



If the energy of a primary particle is less than a few TeV, the shower gets absorbed in the 
atmosphere, while the ACR produced by the secondary charged particles can be detected at 
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the ground level. It was Blackett [|99!|. who suggested that there should be a small contri- 
bution of the order of 10~^ from the ACR produced by CR particles in the atmosphere, to 
the Light Of Night Sky (LONS). Galbraith and Jelley [[I59],[l5l detected the ACR pulses 
associated with EAS which are subjected to large fluctuations in the point of first interaction 
and development of cascade subsequently. The ACR carries information about the point of 
origin (direction) and energy of the primary 7-ray photon. There are some important char- 
acteristics of ACR emitted in an EAS: 



(1) Threshold energy for Cherenkov emission: Equation (3.21 ), can be re-written in terms 
of threshold energy {Emin), a particle must have in order to emit ACR 

Emin = Imin = = (3.28) 



where mo is the rest mass of the particle and n = 0.00029 at sea level. Equation (3.28) 
implies that the minimum energies required to produce ACR at sea level are 21.3 MeV, 
AAGeV and 39.1 GeV for electrons, muons and protons respectively. Thus the light parti- 
cles like electrons are the most efficient in producing the ACR in an EAS. Since the refractive 
index decreases and hence the (3min increases with height, the threshold energy for the pro- 
duction of ACR also increases with height. For electron Emin is 35 MeV at an altitude of 
7.5 km. However, the increase in thereshold energy with altitude does not practically affect 
the fraction of electrons emitting ACR due to the fact that in the early phase of shower de- 
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Fig. 3.9: Radius of ACR pool and Cherenkov pulse duration. 



velopment the mean energy of electrons is significantly higher than in the shower tail. The 
smaller number of electrons at high altitudes produce only a small amount of ACR, while 
most of this is produced in the region of shower maximum [I127II . 

(2) Maximum Cherenkov emission angle: The maximum Cherenkov emission angle 



given by equation (3.22) also depends on height of emission in the following manner 



cos Omax = TT" \ , ^ = 1 " ^70 exp{-h/ho) (3.29) 

i + ?7o exp[ — ii/ iiq ) 

On expanding cos O^ax iX - %^ + %^ - %^ H ) we get 

Omax = ^y2r]^exp{-h/2ho). (3.30) 

Thus the maximum emission angle (or emission angle in general for (3^1) for ACR with 
respect to shower angle depends on refractive index (and height of emission). The angle of 
emission increases with decreasing height. Note that the average scattering angle 6coi due 
to multiple Coulomb scattering which is described by ( 017611 and references therein) 
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Fig. 3.10: Radius of Cherenkov light cone as a function of height of emission. This is shown for sea 
level and Mount Abu altitude (1300 m asl). 



also increases with shower age. Both these effects lead to somewhat increased Cherenkov 



emission angle with decreasing height. Fig. 3.8 shows the variation of Emin for electron and 



maximum angle for ACR as a function of atmospheric height above sea level. 

(3) Cherenkov light pool radius: Let us consider a case wherein an energetic particle 



emits ACR from a height h as shown in Fig. 3.9 The emission is in a narrow cone with 



maximum opening angle 9max around its trajectory resulting in a nearly circular ring (also 



called Cherenkov ring) of radius r which can be determined by using the equation (3.30) 



r{h) = htanOmax = \/'^Vo h exp{—h/2ho). (3.32) 



The value of r{h) is maximum when h = 2ho which is rmax — 128 m. Fig. 3.10 shows the 
dependency of the radius of Cherenkov ring on the height of emission at two observational 
levels (at sea level and at Mount Abu altitude). It is clear purely from geometrical consider- 
ation that the radius of Cherenkov ring shrinks with increasing observational height for fixed 
height of emission. This radius for an observational height hobs can simply be written as 

r{h) = {h- hobs) tan 6max = ^/^{h - hobs) exp{-h/2ho). (3.33) 

It gives the maximum radius of the Cherenkov ring for Mount Abu altitude as Vmax — 
116.5 m. The ACR distribution of an EAS on the ground is the superposition of ACR cones 
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Fig. 3.1 1 : (a) Geometrical model for emission of Cherenkov radiation from 7 and hadron induced 
showers. The stripped region is the main region for the emission of Cherenkov light from 
7-ray showers, (b) Lateral Cherenkov photon distribution for both 7-ray photon (1 TeV) 
and proton (2 TeV) showers. The 7-ray induced shower shows a relatively flat lateral 
distribution of Cherenkov light upto nearly 125 m radial distance from the shower core. 
These figures have been adapted from 111891 . 



from individual shower particles (mostly electrons). In case of an electromagnetic shower, 
ACR illuminates typically a circle called ACR pool with an average maximum Cherenkov 
radius of 120 m at moderate altitudes like Mount Abu. Note that the circle is a special 
case for a vertically incident shower, while for an inclined EAS, the overlapping ACR forms 
an ellipse. Such ACR pools are typically produced by emissions taking place between 10 to 
20 km above sea level. 

(4) Geometrical Model: The ACR as seen at detector level, can be considered to come 



from three portions of a typical ITeV shower shown schematically in Fig. 3.11(a) The 
first portion which contains ~ 25% of total light comes from shower particles at elevations 
between the height of first interaction down to an elevation of 10 km. The mean value of the 
Cherenkov emission angle is about 1° in air. It increases from 0.66° at 10 km to 0.74° at 8 km 
and finally becomes close to 1.4° at sea level. This results in focusing of light on the ground 
into a blurry Cherenkov ring with radius of r ^ 10A;m x 0.012 = 120m for a typical 
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1 TeV 7-ray shower. Multiple Coulomb scattering gives rise to an exponential distribution of 
scattering angles with respect to the shower axis ~ e^^/^° where the characteristic scattering 
angle is about Oq = 0.83 -E~°„^'^ giving a typical value of ~ 5° ([63] and references therein). 
While this angle is much larger than the Cherenkov angle, the exponential distribution is still 
strongly peaked in the forward direction and the Cherenkov ring, while blurred by scattering, 
is still discernable. The Cherenkov ring is most prominent when the shower dies out far 
above the observational level, as in low energy 7-ray showers. The light from the highest 
altitude arrives at the same time as that from the lower altitudes as the particles in the shower 
travel close to velocity of light and compensate for the greater geometrical path travelled 
because light from higher altitudes goes at slower velocity in air. The result is that the light 
in the annulus is strongly bunched in time with a spread of ~ 1 ns ^29511 . 



The bulk of the light (~50%) comes from a cylinder (see Fig. 3. 1 1(a) I of length ~ 4 km 
(nearly between 6 to 10 km height) and radius 21m centered on the shower core; this cylin- 
der contains the shower maximum and hence bulk of the emitted light. The Cherenkov 
photons on reaching ground level form a ring like structure which is also called the hump 
of the Cherenkov pool. The light from this region is a good measure of the total energy and 
is best measured at a distance of nearly 100 m from the shower axis on the ground. The 
angular spread of this light will have a half-width of ~ 0.2° [I295II . For a proton cascade of 
same primary energy the similar column is represented by dashed line. The differences in 
the depth of shower maximum for a 7-ray photon and proton is due to the difference in their 
interaction lengths. 

The last ~25% of the light originates from the few surviving particles of the shower, 
radiating below an elevation of 6 km which is subjected to large fluctuations, and generally 
falls close to the shower axis. 

The intensity of ACR on the ground as a function of distance from the shower core is 



shown in Fig. 3. 1 1(a) In this lateral distribution, any particular radial distance on the ground 
from the core represents a particular height in the atmosphere. By comparing the lateral 
distribution of ACR for different energies it has been shown by Rao and Sinha [I259II that the 
hump feature is due to high energy electrons above 1 GeV. Electrons of lower energies are 
more scattered and do not contribute to the formation of the Cherenkov ring. The intensity 
near the shower axis rises when the tail of the shower of a high energy particle penetrates 
close to the observational level. This is because the light is emitted closer to the the ground 
and has a shorter distance to travel. The lateral distribution of ACR yield on the ground for 



a typical 1 TeV 7-ray shower is shown in Fig. 3.11(b) The 7-ray induced shower shows a 
relatively flat distribution upto the hump, in contrast to a proton shower of a similar energy. 
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(5) Cherenkov pulse duration: In a non-dispersive medium, the wavefront is infinitely thin 
and therefore the duration of ACR pulse is infinitely short. However, in a dispersive medium, 
the Cherenkov angle is different for different wavelengths and the wave trains will spread out. 



The observed duration along the particle track at a distance r' as shown in Fig. |3.9| HI 8911 is 



At = ^— ^ ■ {tan 9'^^^ - tan dmax)- (3.34) 
Another effect that increases the pulse duration is the different heights of emission of the 



Cherenkov photons. Again referring to Fig. 3.9 emissions of Cherenkov photons take place 



at two points A and B, assuming that photons are emitted along the track of the particle, the 
measured time difference between the detection of photons can be written as 

5t = '^^(exp-''/''' - exp-'/'^). (3.35) 
c 

Taking h' = 6 km and h = 10 km, we get 5t ~ 1.3 ns. The total pulse duration will 
be the sum of these two effects. When the development of a 1 TeV 7-ray photon induced 
EAS shower is simulated through the Monte Carlo procedure, the value of Cherenkov pulse 
duration turns out to be ~ 5 ns [I281L 

(6) Differences in ACR distribution of 7-ray and proton induced showers: 

(i) Because of the longer interaction length, the particles producing ACR in proton show- 
ers reach closer to the ground level than 7-ray showers, which gets reflected in Fig. 



3.1 1(b) as the intensity is more near the shower axis for protons. The Cherenkov ring 



also gets smeared out in proton showers. 

(ii) The ACR pool created by 7-rays is more regular in comparison to the protons since 
there are more sources for the Cherenkov emission other than electrons and positrons 
(e.g. muons) contributing to the light pool. Muons produce intense rings in ACR 
distribution on the ground. 

(iii) For the same primary energy, proton showers have a ACR pool less intense than the 
7-ray showers, as only 1/3 of the primary energy (to the first approximation) in proton 
shower contributes in the development of electromagnetic cascade. 

(iv) The Root Mean Square (RMS) of Cherenkov photon density measured over a 10 
area is calculated to be 16% for 7-ray showers and 36% for the proton showers at 
around 1 TeV energy and at 100 m from the shower core [I189L 



These differences are demonstrated in Fig. 3.12 
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Fig. 3.12: Comparision of the lateral distributions of ACR for vertically incident 7-ray and proton 
showers. The left shows the arrival positions on the ground of all Cherenkov photons 
produced in an air shower created by a single 50 GeV 7-ray . The right shows the same 
information for an air shower created by a single 200 GeV proton. This figure has been 
adapted from Ii239il . 

3.5 Detection techniques of extensive air showers 

The basic purpose of detection of an EAS experimentally, is to determine the type (though 
difficult) and properties of primary particle, such as direction, energy and chemical composi- 
tion in case of CRs. The detection methods are mainly based on the energy region of interest. 
In a 7-ray induced EAS, the particles are beamed to the ground and their trajectories retain 
the directionality. The shower front is narrow with a thickness of ~ 1 m at the center. Tight 
longitudinal structure (i.e. fast timing ~ few ns) can be used to measure the shower orienta- 
tion and hence the primary direction. Fast timing allows to separate the signal produced due 
to the shower front from random background processes. Finally the lateral distribution of air 
showers allows detector spread out on the ground to be sensitive to particles whose initial 
trajectory would not have intersected a detector in the absence of the atmosphere. Thus the 
effective area of an EAS detector is much larger than its physical area. 

At energies more than few TeV, sufficient number of secondary particles are able to 
reach the ground, which can be detected with a suitable detector e.g. array of detectors based 
on scintillation counters. Such a detector is normally called an air shower array. The fast tim- 
ing between detector elements is used to estimate the direction of the primary. The shower 
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core is approximated by the location of the highest particle density and the lateral distribu- 
tion is reconstructed from measurements of the particle density at different radii from the 
core [I239II . The advantage of this technique is the 100% duty cycle, while angular resolution 
is poor and has only modest possibility of 7-hadron separation. 

A small fraction of the primary energy is emitted as isotropically distributed fluorescent 
light (yield is ~ 10^^ of the total energy 1112711 ) from the excitation of nitrogen molecule^ 
Operating at rather high threshold energies (> 10^^ eV), fluorescent detectors can be oper- 
ated only during dark periods. Detection of neutrinos is also possible, but due to their weakly 
interacting nature with matter, these detectors are not very successful. 

On the other hand, the ACR produced by an EAS can be measured very efficiently above 
energies of more than few tens of GeV. The Cherenkov photon spectrum peaks near UV re- 
gion (1 / dependence) and the sky is also transparent to light between 300 to 600 nm, most 
of the produced Cherenkov photons reach the ground under good atmospheric condition^ 
Preserving the directional information contained in the shower particles, the Cherenkov pho- 
tons are strongly beamed in the forward direction and arrive at the ground with a tight time 
structure (~ 5 ns). This permits the detection of ACR by fast photosensitive detectors amidst 
the LONS. This technique is discussed in detail below. 

3.5.1 Atmospheric Cherenkov technique 

This technique has been used for more than forty years to search for point sources of VHE 
radiation. The details of the technical aspects of the ACT can be found in [11151 , 12941 , 
12391 . The first generation systems consisted of just a single light detector in the focal plane 
of a search-light mirror coupled to fast pulse counting electronics. These telescopes are 
characterized by the mirror collection area (A), the reflectivity (R), the solid angle (f2) and 
the signal integration time (r). Even with a simple detector like A = 2 m?, R = 85%, 
n = 10^^ and r = 10 ns, it is possible to detect the ACR signal from EAS (of few TeV 
energy) with high efficiency. 

3.5.2 Trigger threshold of the telescope 

An atmospheric Cherenkov telescope generates a trigger whenever an ACR signal hits the 
detector with a sufficient strength from both 7-ray photon or proton induced EAS amidst 
the noise of NSB photons. Characteristic features of an ACR pulse e.g. narrow pulse width 

"^during de-excitation, a distinct well known line spectrum 337 nm iV2-line is emitted near UV. 

^For example: 1 TeV 7-ray shower produces on an average ^ 3 x 10^ Cherenkov photons at ground. 
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(~5 ns), limited angular size (< 1°) on the ground and nature of its photon spectrum|^ facil- 
itate its detection. Within the ACR pool (circle of approximately 125 m radius), the number 
of Cherenkov photons is proportional to the energy of the primary 7-ray . Defining a scaling 
factor for Cherenkov photon yield (y^) as 

yj = P'y/E (3.36) 

where is photon density of Cherenkov photons on the ground and E is energy of the 
primary particle. The value of depends on energy e.g it is 0.033 photons /m'^ /GeV at 
25 GeV for wavelength between 300 and 550nm and increases to O.OQb photons / m? /GeV 
at 1 TeV [I239II . It also increases with height because of lesser atmospheric attenuation. 

The flux of night sky photons is highly location dependent and has a typical value of 
2 X 10^^ photons / m? / s/sr for a good quality site in the wavelength region of 300 to 600 nm 
( 1123911 and references therein). A telescope with 1° field of view and 10 ns trigger formation 
time would receive ~ 5 photons / from the LONS while a much higher expected value of 
~ 65 photons / m'^ from a 1 TeV 7-ray shower. Thus, a 1 TeV shower should be detectable 
above NSB with a telescope having 1 mirror area. 

The signal produced by Cherenkov photons in terms of a number of photoelectrons in a 
photosensitive instrument, which is usually a photomultiplier tube (PMT) can be written as 

S = p^AR7]pmt = y-fE AR7]pmt (3.37) 

where rfpmt is the quantum efficiency of the PMT. The ACR must be detected above the 
fluctuations in the NSB in the integration time r. The noise level can be expressed in terms 
of fluctuations as 

LONS ^ A RrjpmtT (3.38) 

where (pLONS is the NSB photon flux and is the solid angle subtended by a PMT. Note 
that, in actual case the wavelength dependence of both Cherenkov and LONS production and 
collection should also we considered. The energy threshold of an atmospheric Cherenkov 
telescope can be defined in a number of ways. At the moment the most appropriate defini- 
tion is the minimum 7-ray energy for which the signal to noise ratio (S/N) is sufficient to 
adequately trigger the telescope. The smallest detectable light pulse is therefore inversely 

^Cherenkov light peaks at short wavelength (blue/UV) whereas LONS peaks at longer wavelength 
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proportional to S/N, i.e. 



As is clear from equation (3.39), the energy threshold of a Cherenkov telescope can be 
reduced by working at darker sites, reducing field of view and integration time and at the 
same time increasing the mirror collection area, its reflectivity and quantum efficiency of the 
PMT. 



3.5.3 Flux sensitivity of the telescope 

The most troublesome factor in 7-ray astronomy is the huge background of CRs. The ACR 
distribution from these hadronic showers is similar to that from 7-ray showers so that first- 
generation Cherenkov telescopes were unable to distinguish between the two. In the field 
of view of a simple telescope, whose solid angle is optimized for 7-ray detection, the back- 
ground CR events are 10^ times as numerous as the strongest steady 7-ray discrete source 
thus far detected. Since the arrival directions of CRs are isotropic due to interstellar magnetic 
fields, a discrete source of 7-rays can stand only as an anisotropy in an otherwise isotropic 
distribution of air showers. However, there are very distinct differences in the cascade de- 
velopment in electromagnetic and hadronic showers, as already discussed and can also be 
found in [11831 , 11911 , that make the study of sources of VHE 7-rays with Cherenkov tele- 
scopes possible. 

The flux sensitivity basically defines the performance of a Cherenkov telescope and is 
determined by its ability to detect a 7-ray signal over the CR background. Defining the 1]^ 
and r]cr as the 7-ray and CR retention factors, after using some discriminating tools for two 
types of shower events, a "figure of merit" called quality factor (Q) of a telescope can be 
defined as: 




(3.40) 



The CR background has a power law type spectrum, which can be expressed in the form of 
integral flux as 

0er(>^) = C,,^-"=^ (3.41) 

In the range of interest acr = 1-7. Similarly, the 7-ray source energy distribution can be 
assumed to have the form 

(i)^{> E) = E-"^ (3.42) 
aj can have values from 1 to 3 and is generally assumed to increase with energy i.e. the 
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spectrum steepens [|295i 

If S is equal to the number of 7-rays detected from a given source in observation time 
T, and A-y is the collection area for 7-ray detection, then the 7-ray signal for some threshold 
energy E can be written as 

S = (f)y{E)Ay{E)T. (3.43) 
The background CRs recorded during the same time is given by 

B = (j)criE)Acr^T (3.44) 

where Acr is the collection area for the detection of CRs of energy E. Defining A^^r as the 
minimum number of standard deviations at which a source must be detected to be believ- 
able|3 we have [|239l 

AT^ OC ^ = i^(0.85-a,) ME) yl/2 



As is evident from the equation (3.45), to maximize i.e. optimizing the sensitivity, the 
product Acr^ should be minimized, while at the same time maintaining a reasonable col- 
lection area for 7-rays , which neccessitates the requirement of an efficient technique for 
background rejection. The source related information can also be used to optimize sensitiv- 
ity, e.g. steady sources like Crab Nebula should be observed for longer period. For typical 
sources, the telescope can be operated at lowest possible energy threshold, which can also 
help the detection of sources having cutoffs at higher energies. However, sources with ex- 
tremely hard spectra (a^ < 0.85, though rare), require observations at higher energies. 



3.5.4 Methods of rejecting cosmic ray background 

Since more than 99% of the air showers produced in the atmosphere are caused by the 
isotropic CRs, the study of 7-rays requires an efficient technique recognising and rejecting 
the hadronic showers and accepting only 7-ray showers. The characteristics of the ACR pro- 
duced in air showers induced by 7-rays and hadrons are similar to a large extent. However, 
as already pointed out Monte Carlo simulations show that there are perceptible differences 
between the two types of primaries, which makes their discrimination possible. Many tech- 
niques to achieve a good 7-hadron separation exist but all of them are based on geometrical 
considerations and involvement of different processes in the way the two kind of showers 

'^Ncr for a reliable source detection is generally taken as 5 
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develop in the Earth's atmosphere. Here are some of the characteristics of the ACR that can 
be used for their segregation: 

Lateral distribution: The ACR is sampled by an array of small telescopes. The arrival 
direction of the shower is determined through the arrival time information of the shower 
front at different detectors. The photon density is also sampled by digitizing the Cherenkov 
pulse. A 7-ray induced shower is expected to have a flatter distribution of photons as com- 
pared to hadrons. Examples of such telescopes are given in the next section. The flatness 
parameter is used for background rejection with a quality factor of ~ 1.8 [|93l . 

Pulse duration: The shape of the Cherenkov pulse is related to the development of the 
shower cascade, e.g. the rise time of the pulse reflects the longitudinal development and 
decay time is a measure of attenuation of the cascade below shower maxima [1281 II . Width of 
the Cherenkov pulse represents the production profile of Cherenkov photons. Some of these 
properties are exploited by the Pachmarhi group llil8L 

Spectral content: Hadron showers are comparatively richer in muon contents, which pro- 
duce ACR closer to the observational level. The UV part of ACR from a 7-ray shower 
gets almost absorbed in the atmosphere, thus making it possible to differentiate between 
7-ray and hadron based on their UV to visible light ratio. This approach was used by the 
Crimean group 113061 . 

Shower image parameters: In addition to other details (discussed in next section) the 
Cherenkov image also contains the important information about the nature of the primary 
particle e.g. the image that results from a typical 7-ray shower is elliptical and compact with 
an orientation that points towards the center of the field of view. In contrast, Cherenkov im- 
ages from hadrons are less regular and randomly oriented in the focal plane II182II . This is the 
imaging technique explained in next section and has been used in the work being presented 
in this thesis. This technique provides very high quality factor with the rejection of ~ 99.7% 
hadrons, while retaining ~ 50% 7-rays . 

3.6 The imaging technique 

The atmospheric Cherenkov telescopes fall into two broad classes: Wavefront sampling at- 
mospheric Cherenkov telescopes and lACTs. The former records the integrated while the 
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Fig. 3.1 3: Basic principle of the imaging technique, illustrating the detection of an EAS. This figure 
has been adapted from [.127.1 . 

later records the differential information of the ACR. Wavefront sampling instruments 
make use of a large array of mirrors that cover a large fraction of ACR pool and record 
the light intensity at a number of points on the ground. Examples of these instruments are 
THEMISTOCLE experiment il, AGS AT MM, Solar-II / CACTUS (2871, STACEE [fTTTl . 
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Fig. 3.14: Sketch of a shower and its image, (a) Longitudinal distribution of the shower related to 
the length of the image, (b) Lateral distribution of the shower and determination of the 
width of the image. This figure has been adapted from II155L 
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On the other hand, the development of atmospheric Cherenkov imaging technique gave 
the first effective discrimination between 7-ray and background hadron showers ( II295II and 
references therein). In the following, we explain the working principle of the imaging tech- 
nique. An array of PMTs in the focal plane of a large optical reflector constitutes a camera to 
record ACR image for each triggered shower. The camera is triggered when a preset number 
of PMTs detect a light level above a set threshold within a short integration time. The light 
level in all the pixels is then recorded digitally and the image is analyzed offline to determine 
whether it has the expected characteristics of a 7-ray shower with a point of origin at the 
center of the field of view. Discrimination against the background CR showers is based on 
geometry of the image e.g. 7-rays arriving parallel to the optic axis of the telescope have 
roughly elliptical images which appear to originate from the center of the camera. 



As depicted in Fig. 3.13 due to the mirror and camera geometry, Cherenkov photons 
viewed within the opening angle of a pixel are imaged onto the same point in the camera, and 
therefore the camera coordinates are usually defined in angular coordinate system. Since the 
height of emission of photons, governs the angle of emission, the different emitting regions 
of the shower are mapped at different positions in the camera plane. The image formation 



of different parts of the shower in the camera is shown in different colors in Fig. 3.13 The 
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upper part of the shower is mapped onto the camera closer to its center, whereas the lower 
part of the image is mapped further away. The shower axis which defines its orientation 
is mapped onto the camera plane as a straight line called image axis pointing towards the 
shower impact position on the ground. The image of the shower thus recorded in the camera 
contains information about the development of the shower and the properties of emitted 
ACR. It is possible therefore to infer the information of the primary particle, such as its 
energy, direction of incidence and also its nature. The energy of the primary particle is 
largely determined from the light content of the image, while the direction is determined 
through the orientation and shape of the image. Thus the procedure of inferring the properties 
of the primary particles based on the information from the image formed in the camera is 
called the imaging technique. 



Fig 3.14 demonstrates, how the geometry of an EAS is related to the shape of the light 



distribution imaged in the camera. Fig 3.14 ^a) represents the view of the plane containing 
the telescope and the shower axis, wherein length of the camera image determines the angle 
(02 — 0i) under which the shower is seen. If the shower is viewed over longer heights 
(points A and B are separated farther away) or shower is viewed from larger distance (large 



core distance), the camera image results in an elongation. Fig. 3.14;b) demonstrates the 
perpendicular view i.e. the width of the camera image determines the lateral extent of the 
shower. 

Based on the geometrical relations, explained above A. M. Hillas in 1985 HI 8211 con- 
structed a set of image parameters called Hillas parameters, which are very efficient tools 
for reconstruction of an EAS and improved the capabilities of imaging technique to distin- 
guish between 7-ray and hadron initiated showers. The Hillas parameters are still providing 
the best results in the field of VHE 7-ray astronomy. We will discuss the Hillas parameters 
in Chapter 5. 



Chapter 4 



TACTIC telescope 



4.1 Introduction 

In this chapter, we describe the instrumentation aspects of the TACTIC (TeV Atmospheric 
Cherenkov Telescope with Imaging Camera) ||9ll imaging element. We discuss the main 
features of its various subsystems and its overall performance with regard to (/) tracking ac- 
curacy of its 2-axes drive system, (ii) spot size of the light collector, (///) back-end signal 
processing electronics and topological trigger generation scheme, (zv) data acquisition and 
control system and (v) relative and absolute gain calibration methodology. The TACTIC 
has been in operation at Mount Abu (24.6° N, 72.7° E, 1300 m asl), a hill resort in West- 
em Rajasthan, India, for the last several years to study TeV 7-ray emission from celestial 
sources. The mean NSB brightness levels at Mount Abu is found to be ~ (1.34 ± 0.50) x 
10^'^ photons m^"^ s^^ sr^^ (pre-monsoon) and ~ (1.54±0.55) x 10^'^ photons m^"^ s^^ sr^^ 
(post-monsoon). The telescope uses a tessellated light-collector of area ~ 9.5 and de- 
ploys a 349-pixel imaging camera, with a uniform pixel resolution of ~ 0.3° and a total 
field-of-view ~ 6° x 6°. It records the images of the EAS via ACR which is produced by 
the secondary relativistic particles of the EAS. Using a trigger field of view of ~ 3.4° x 3.4° 
(corresponding matrix of 11 x 11 pixels) the telescope records a CR event rate of ~ 2.5 Hz 
at a typical zenith angle of 15°. The photograph of the TACTIC imaging telescope is shown 



in Fig. 4.1 



4.2 Mechanical structure of the telescope 

The main sub-assemblies of the TACTIC telescope are : mirror basket, mirror fixing/adjusting 
frames, zenithal and azimuthal gear assembly, encoders/motors for its two axes, camera sup- 
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Fig. 4.1 : Photograph of the TACTIC imaging element. 



port boom assembly and the PMT-based imaging camera. A three dimensional truss type 
structure has been used to support the mirror frame. The mirror basket is a 3 layer welded 
mildsteel tubular grid structure fabricated in 3 parts which are bolted together. The basket 
is provided with a central tie rod which extends into short stub shafts at the two ends which 
are housed in bearings. The individual mirror facets are supported on three levelling studs, 
so that desired inclination of the specific mirror with respect to the telescope axis can be 
achieved. The weight of the moving part of the telescope is around 6.5 ton. The complete 
basket assembly is held by the fork of the telescope which transmits the load through the cen- 
tral roller thrust bearing to the telescope foundation. The horizontal part of the fork frame 
is clamped on to a large gear of the azimuthal drive rotary table of the telescope. The thrust 
bearing using hardened rollers was specially designed and fabricated in-house. A central pin 
and a taper roller bearing, housed in the large gear is used to ensure proper rotation about 
its vertical axis. The base structure has been anchored to the ground by 4 foundation bolts. 
The zenithal motion to the telescope is given from only one end of the basket. A common 
shaft assembly fixed in self aligning ball bearing having a large gear fixed on it, provides the 
zenithal rotation of the mirror basket. Five stage gear box providing a total speed reduction 
of 5189.14 has been used in the zenithal drive of the telescope. Likewise, a four stage gear 
box leading to a speed reduction of 6348.94 has been used in the azimuth drive. A large ca- 
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pacity circular cable drag chain has been provided for easy and free movement of the nearly 
700 signal and high voltage cables. 

Deflection analysis using FEM software package has also been performed at various 
zenith angles to arrive at the shift in the focal point position due to structural deformations 
of the mirror basket frame, mirror holding attachment and support booms. The results of this 
study indicate that a shift of < 3 mm in the focal point position can be produced over the 
zenith angle range of 0° to 70°. 



4.3 Drive control system 

The large light collector aperture (~ 3.5 m) coupled with the attendant large telescope weight 
(~ 6.5 ton) have led to the choice of an altitude-azimuth (alt-azm) mounting for the TACTIC 
telescope, as against the comparatively simpler equatorial mounting. The main advantage of 
the alt-azm mount is that the telescope weight is supported uniformly on a horizontally- 
placed central thrust-bearing. For an alt-azm mount, the immobile axis is set towards the 
zenith and the telescope guidance is provided by the azimuth and zenith angles. As the 
source co-ordinates (right ascension, a and declination, 5) are known in the equatorial co- 
ordinate system, a co-ordinate transformation is required to find the position and speed of 
the celestial source along either axis of the telescope. The general equations for such a 
transformation are given by: 

cosS sink , , , , 

^ = — Z u X (4.1) 

smcp cosh coso — coscp smo 

cosz — sin4) sin5 + cos5 cos(f) cosh (4.2) 

dz/dh — coscf) sin A (4.3) 

dA/dh — sin(f) + coscf) cos A cotz (4.4) 

where (p is the latitude of the observatory (for Mount Abu (p = 24.63°), z is the zenith angle 
measured from local vertical, A is the azimuth angle measured from the geographical south 
with positive sign for the angle towards the west and h is the source hour angle. The hour 
angle of the source is related to the local sidereal time (LST) by the equation 



h — LST — a. 



(4.5) 
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The design of the TACTIC drive control system is based on the CAMAC standard. The tele- 
scope uses two 100 cm hybrid stepper motors (Pacific Scientific make, Model H 31NREB, 
NEMA Size 34) for driving its azimuthal and zenithal axes through multistage gear-trains. 
The maximum azimuthal speed of the TACTIC telescope is 100 rad/rad (equivalent to 
1500" hour^^), and the resulting blind spot has a diameter ~1.2°. A GPS based CAMAC- 
compatible digital clock (Hytec Electronics Ltd. make; GPS92) with a resolution of ~10ns 
and absolute time accuracy of ~100ris, is used to compute the source co-ordinates in real 
time. The new co-ordinates of the source are calculated after every second while tracking a 
candidate source. More details regarding various hardware components and related software 
of the telescope drive system are discussed in [12821 . 12811 . 120511 . 

The tracking accuracy of the telescope is also checked on a regular basis with so called 
"point runs", where a reasonably bright star, having a declination close to that of the candi- 
date 7-ray source is tracked continuously for about 5 hours. The point run calibration data 
(corrected zenith and azimuth angle of the telescope when the star image is centered) are 
then incorporated in the telescope drive system software so that appropriate corrections can 
be applied directly in real time while tracking a candidate 7-ray source. 



4.4 Light collector design 

The TACTIC light collector uses 34 front-face aluminium-coated, glass spherical mirrors 
of 60 cm diameter each with the following characteristics (/) focal length ~ 400 cm, (zz) 
surface figure ~ few A (zzz) reflection coefficient >80% at a wavelength of ~ 400 nm and 



(zv) thickness 20 mm to 40 mm. Fig. 4.2 ^a) shows the projection of the mirror layout onto 
a plane transverse to the symmetry axis of the basket. The shorter focal length facets are 
deployed close to the principal axis of the basket while the longer focal length facets are 
deployed around the periphery. The peripheral mirrors have the effect of increasing the 
overall spot size as they function in an off-axis mode. In order to make the design as close to 
the Davies-Cotton design as possible, we used longer studs on the mirror frame structure to 



raise the pole positions of the peripheral mirrors. Fig. 4.2 ^b) shows the measured reflection 
coefficients of individual mirrors. The reflection coefficient measurements of the mirrors 
were made using a reflectometer (Dyn-Optics make; Model 262). The alignment procedure 
of the various mirror facets is discussed in [12841 . 1205 [[ . After mirror aligment, a common 
focus with the minimum possible image spread was obtained at a focal plane distance of 
386 cm instead of at 400 cm, as would have been expected for a standard paraboloid or a 
Davies-Cotton design of the reflector. This value of focal plane distance is chosen on the 
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Fig. 4.2: (a) Projection of the 34 mirror facets of the tessellated light collector of the TACTIC tele- 
scope. The shaded circles represent the peripheral mirrors that have been raised along the 
z-axis of the light collector for obtaining a better point spread function (b) Reflection coeffi- 
cients for all the 34 mirrors used in the telescope. Seven randomly selected locations on the 
mirror surface have been used to quantify the mean reflection coefficient of the mirror. The 
error bars represent standard dispersion (Icr value) of the 7 reflection coefficient values. 



basis of the simulation results ^28411 . The alignment of the mirror facets is further confirmed 
by observing a bright star image at the focal plane. Gross misalignment in any of the facets 
is easily identified as it results in multiple images being seen on the focal plane. 

In order to evaluate the optical quality of the light collector experimentally, the telescope 
was pointed towards the bright star (^-Tauri and its image recorded by monitoring the anode 
current of the central pixel of the imaging camera. The anode current versus angular offset 



plot is shown in Fig. |4.3| ^a). The point-spread function shown has a FWHM of ~ 0.185° 
(= 12.5 mm) and -Dgo ~ 0.34° (= 22.8 mm). Here, Dgo is defined as the diameter of 
circle, concentric with the centroid of the image, within which 90% of reflected rays lie. An 
image of the star Sirius recorded at the focal plane of the telescope has also been shown 



in Fig. |4.3[b) and it has superimposed on it two circles which correspond to the diameter 
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Fig. 4.3: (a) Measured point spread function of the TACTIC telescope light collector, (b) Photograph 
of the image produced by Sirius. Circles superimposed on the image have diameters ~ 
0.31° and ~ 0.34° and represent the diameter of the camera pixels and Dgo, respectively. 



of the pixel and measured DgQ value calculated on the basis of Fig. 4.3 'a). The value of 
Dgo ~ 0.29° (= 19.3 mm), predicted on the basis of the simulation for an incidence angle 
of 0°, matches reasonably well with the measured value mentioned above. Other details 
regarding the ray-tracing simulation procedure and comparison of the measured point-spread 
function of the TACTIC light collector with the simulated performance of ideal Davies- 
Cotton and paraboloid designs are discussed in [I284II . 



4.5 The imaging camera 

The imaging camera with an angular resolution of 0.3° has been used in the TACTIC focal 
plane. The camera frame is made up of two 5 mm thick aluminium plates in which 19 mm 
diameter holes are drilled at a pitch of 22 mm corresponding to the locations of the 349- 
pixels. The PMT is held in place by a metallic collar fixed to its socket which in turn is held 
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to the rear plate by a specially designed fastner. The central pixel of the camera which is on 
the principal axis of the light collector is used for checking the alignment of the mirror facets 
and the tracking/pointing accuracy of the telescope. The pixels are numbered sequentially 
clockwise from the central pixel which is designated as pixel number 1. 

The camera uses 19 mm diameter PMT (ETL-9083 UVB), with bialkali photocathode 
having maximum quantum efficiency of ~27% at 340 nm and UV glass for the window to 
enhance its sensitivity in the 280 — 300 wavelength band as well. The 10 stage linear 
focussed PMT has a rise time of ~ 1.8 ns which is compatible with the time profile of the 
Cherenkov pulse. A low current zener diode-based voltage divider network (VDN) is used 
with the PMT. This VDN design fi92il has the advantage of ensuring stable voltages at the 
last two dynodes with VDN current of only about 240 fxA which is a factor of 5 less than the 
minimum current recommended for a resistive VDN. The VDN uses negative voltage and the 
photocathode is at a high voltage of 1000 — 1400 V while the anode is at the ground potential. 
The VDN of the PMT is permanently soldered to its socket and two RG174 coaxial cables 
from each VDN circuit board are terminated with coaxial connectors on the connector panels 
fixed to 4 sides of the camera. 

The Compound Parabolic Concentrator (CPC) shape was chosen for the light guides 
(made of SS-304) to ensure better light collection efficiency and reduction in the NSB light 
falhng on the PMTs. Some of the important geometrical parameters of the CPCs used in 
the TACTIC telescope camera are the following: entry aperture ~ 21.0 mm, exit aperture ~ 
15.0 mm, acceptance angle ~ 45.58° and height ~ 17.6 mm. The light collection efficiency 
of the CPC, which includes both the geometrical collection efficiency and the reflectivity of 
the surface was experimentally measured to be ~ 65%. 

4.6 Backend signal processing electronics and trigger 
generation 

4.6.1 Signal processing electronics 

The image of a typical Atmospheric Cherenkov Event (ACE) is registered in the form of 
varying amplitude pulses of~ 3 — 60mV^ produced by a group of 5 — 20 pixels of the camera. 
These voltage pulses are brought to the control room, using 55 m long high quality RG 58 
coaxial cables. In-house developed fast NIM Hex amplifier modules with a user selectable 
gain range of 2 — 50 and amplitude discriminator modules of 50 — 500 mV range are used for 
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amplification and threshold selection of the PMT signals. While the multichannel fast NIM- 
based amplifier and fixed threshold discriminator modules have front panel adjustments for 
gain and threshold respectively, the charge content and the scaler rates are directly read off 
the CAMAC bus. The complete back end instrumentation shown in Fig. |43 is based on 
inhouse developed medium channel density modules and is housed in seven 19 inch racks of 
36 U (1U=1.75 inch) height (Fig. |4.4| ). One of the two outputs of a discriminator channel is 
used for monitoring the Single Channel Rate (SCR) with a CAMAC scaler while the other 
output is connected to the TACTIC Trigger Generator (TTG). The outputs from each of the 
independently operating TTG modules are then collated in an Event Handler (EH) which 
generates the 22 ns duration gate pulse and interrupts the data acquisition for reading the 
charge ADC data from all the 349 pixels. The final trigger pulse is also used for latching the 
GPS -based clock. The scalers and charge-to-digital converters (CDC) for the 349 channels 
use 5 CAMAC crates each. Each of these crates is controlled using an in-house developed 
multi-crate CAMAC controller and five such controllers are daisy chained and connected 
to a data acquisition Personal Computer (PC). A similar strategy following a custom built 
standard has been used for the computer-programmable High Voltage (HV) units. 



4.6.2 Trigger generation 

The imaging camera uses a programmable topological trigger [|90l which can pick up events 
with a variety of trigger configurations. As the trigger scheme is not hard wired, a number 
of coincidence trigger options e.g Nearest Neighbour Pairs (NNP), Nearest Neighbour Non- 
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Fig. 4.5: Block diagram of the back-end signal processing electronics used in the TACTIC imaging 
element; TTG - TACTIC Trigger Generator; PCR-Prompt Coincidence Rate ; CCR- Chance 
Coincidence Rate. 



collinear Triplets (3NCT) and Nearest Neighbour Non-CoUinear Quadruplets (3NCQ) can 
be generated under software control. The layout of the 349-pixel imaging camera, which 
can use a maximum of 240 inner pixels for trigger generation is depicted in Fig.|43'a). The 
trigger criteria have been implemented by dividing these inner 240 pixels into 20 groups of 
12 (3 x4) pixels. A section of 5 such groups is connected to a TTG module and a total of 4 
TTGs are required for a maximum of 15 x 16 matrix of trigger pixels. The trigger scheme has 
been designed around 16k x 4 bit fast static RAM (Toshiba make TC55B417; access time of 
< 8 ns). Each of the 4 TTG modules uses 5 memory ICs (indicated by M1-M5 in Fig. 4.6[ b)) 
and has horizontal cascading built into it. A pictorial representation of the horizontal cas- 



cading followed within each TTG module is shown in Fig. 4.6 ^b). Vertical cascading has not 
been provided in the TTG due to non- availability of more than 3 outputs from the discrimi- 
nator. The loss of events as a result of the absence of vertical cascading has been estimated 
to be about 12% through Monte-Carlo simulations and actual on site experiment ni93L The 
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Fig. 4.6: (a) Layout of the 349-pixel imaging camera of the TACTIC telescope. Shaded pixels in 
the center and the same filled with circles at the periphery of the camera represent PMTs 

for which single channel rate is not stabilized, (b) Pictorial representation of the horizontal 
cascading followed within each trigger generator module. Each trigger generator module 
can handle 60 channels by using 5 memory ICs (M1-M5) of type TC55B417. 



memory address lines are connected to the CAMAC address lines and the front panel recep- 
tacles (for connection to the discriminator outputs) through two sets of tri-state buffers. The 
TTG operation starts with writing of the data, as per a user defined topology, from a disk 
file into each of its memories under CAMAC control. Once programmed, the TTG outputs 
follow the event topology as described earlier. Apart from generating the prompt trigger, the 
TTG has a provision for producing a chance coincidence output based on ^^C2 combinations 
from various groups of closely spaced 12 channels. This chance coincidence output is used 
as a system monitor for evaluating its overall functioning during an observation run. Mon- 
itoring of the Chance Coincidence Rate (CCR) has also helped in keeping a close check on 
the operation of the telescope and the quality of the collected data. Other details regarding 
the design, implementation and performance evaluation of the programmable TTG for the 



Chapter 4. TACTIC telescope 



92 



TACTIC telescope are discussed in [I193L 

4.6.3 Stabilization of single channel rates 

A cost-effective method for operating the imaging camera of the TACTIC 7-ray telescope 
at stable SCR and safe anode current values is being used despite variations in the LONS 
experienced by the individual pixels from time to time ll95ll . The camera operates 9 central 
pixels with fixed high voltages and the remaining pixels at different high voltages to ensure 
their operation within a pre-determined SCR range. The purpose behind using the central 9 
pixels of the camera at fixed high voltages is to facilitate the gain normalization (flat fielding) 
of the remaining 336 camera pixels, so that the event sizes (sum of CDC counts in the clean, 
flat-fielded image) recorded during a night's observations can be directly compared to one 
another. Operation of the pixels in a narrow SCR band has the advantage of ensuring a stable 
CCR which can be used as a system diagnostic parameter. An elaborate algorithm |95l 
has been developed to monitor the SCR rates of all pixels using the CAMAC front ends 
and ensure their operation within a narrow range despite changes in the NSB light level 
incident on them due to changes in the sky brightness and star-field rotation. The algorithm 
also ensures that all the pixels of the camera operate within safe anode current ranges. The 
feedback loop of the algorithm changes the HV to the various pixels using a multi-channel 
HV unit which has a resolution of 1 V. The decision of operating a pixel under enhanced light 
levels is solely based on the comparison of the SCR and applied high voltage with reference 
data generated under controlled light level conditions. The PMTs of TACTIC telescope are 
operated with a safe anode current limit of < 30 fiA. A detailed description of the SCR 
stabilization scheme can be found in [|95l . 

4.7 Data acquisition and control system 

The data acquisition and control (DAC) system of the telescope has been designed around 
a network of PCs running on the QNX (version 4.25 lfT6l ) Real-Time Operating System 
(RTOS). The software is designed for the real time acquisition of event and calibration data 
and on-line display of telescope status in terms of Prompt Coincidence Rates (PCR), CCR 
and the functional status of each of the 349 pixels of the camera. 
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Fig. 4.7: Block diagram of the multi-node PC-based data acquisition and control system. 

4.7.1 System architecture 

As the number of jobs to be handled by the DAC system is fairly large in addition to the 
high throughput of scaler and CDC data from the 349 channels, we have used a network of 
3 PCs equipped with QNX RTOS to perform the various functions of the data acquisition 
and control. The QNX RTOS was chosen for its multi-tasking, priority-driven scheduling 
and fast context switching capabilities. In addition, the QNS RTOS also provides a powerful 
set of Inter Process Communication (IPC) capabilities via messages, proxies and signals. 



The DAC of the TACTIC is handled by a network of three PCs as shown in Fig. 4.7 While 
one PC is used to monitor the SCRs and control the HV to the PMTs, the other PC handles 
the acquisition of the event and calibration data and the programming of the TTGs. These 
two front-end PCs, referred to as the rate stabilization node and the data acquisition node 
respectively, along with a master node form the multinode DAC network of the TACTIC 
Imaging telescope. All executable routines stored on the master node are spawned on to the 
other two front-end nodes as and when required. The same network is extended to two more 
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LINUX-based PCs which are used for on-line data analysis and archiving. An EH module 
described in the next section, controls the whole process of data acquisition. 

4.7.2 Event handler: Interface between CAMAC hardware and the ap- 
plication software 



The EH module shown in a block diagram form in Fig. 4.8 provides the link between the 
TACTIC hardware and the application software. The EH accepts the ACEs, calibration and 
chance trigger outputs from various TTG modules and interrupts the front end data acquisi- 
tion node. The interrupt driven system uses two interrupts at levels 3 and 4 for better data 
segregation of the various event types of the telescope as two different interrupt service rou- 
tines are used for the data acquisition. The system acquires the 349 channel CDC data for 
the trigger selected ACEs, relative calibration flashes generated by the calibration LED and 
sky pedestal events, in addition to CDC data for the 4 absolute calibration channels. The HV 
and scaler data are also logged continuously, though at a much lower frequency. 

The acquisition of CDC data for an ACE is handled at interrupt level 3 while the cali- 
bration and sky pedestal data are acquired using interrupt level 4. The switch over of these 
interrupt levels is managed by the EH hardware under the control of the data acquisition soft- 
ware and works in the following manner. The rate stabilization node detects the stabilization 
of the camera and informs the data acquisition node accordingly which in turn sets a flag in 
the EH resulting in the change in the interrupt level from 3 to 4. A logic pulse is also sent 
to the LED driver which produces a 20 ns duration light flash using the LED (Nichia, Japan 
make type SPB 500) for the purpose of Relative Gain Calibration (RGC) of the 349 pixels 
of the camera. This flash which is sensed by all the operating pixels of the camera produces 
a calibration output from the TTG modules and interrupts the DAC system at the interrupt 
level 4 selected earlier. At event occurrence the EH also generates a TTL output for latching 
the system clock and a 20 ns wide NIM pulse for gating the CDC modules. The hardware 
interrupt to the data acquisition node is generated after a delay of ~ 30 /is to account for the 
CDC conversion time. After every event the EH goes into a blocking mode until the CDC 
data of all the channels have been read out. After the generation of the preset number of 
calibration pulses, the EH sends out user programmed number of false triggers to acquire the 
sky pedestal data from all the channels. Even though the RGC data and the sky pedestal data 
are acquired using the same interrupt level, their segregation is possible because of the fixed 
number of these events. After the acquisition of the preset number of sky pedestal events 
the control software resets the EH flag leading to the selection of interrupt level 3 for the 
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Fig. 4.8: Block diagram of the event handler module which provides the interfacing between the 
CAMAC hardware and application software of the data acquisition and control system of 
the TACTIC Imaging Element. 



acquisition of the ACE data. This process of interrupt level switching continues during the 
course of observation. 



4.7.3 Data acquisition and control software 

The multinode DAC software is designed for the real-time acquisition of event and calibra- 
tion data and on-line display of telescope status in terms of PCR, CCR and the status of each 
of the 349 pixels of the camera. As each specific job is handled by an independent process, 
the software comprises of about 6000 lines of C-code distributed over 15 routines and de- 
vice drivers. Ten concurrent processes running across the network of 3 PCs under the QNX 
RTOS handle the various aspects of DAC. All executable routines stored on the master node 
have priorities (0 to 31) assigned to them and are spawned on to the other nodes as and when 



required. Some of the main processes of the data acquisition system are listed in Table 4.1 



The software is also provided with all necessary safeguard measures to ensure trouble-free 
operation of the telescope. Furthermore, the system also provides single-point monitoring 
and control from a designated node with elaborate GUI facilities for monitoring the func- 
tions of various sub-systems. Quasi on-online analysis of the data recorded by the telescope 
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Table 4.1 : List of processes involved in the TACTIC data aquisition and control system. 



Process name 


Description 


Priority 
(node) 


stJactic 


Spawns and controls the execution of all other processes 


23 
(2) 


gain_ctl 


Monitors the scaler rates and controls the high 
voltage of all the channels 


20 

(3) 


da_ctl 


Controls the execution of cal_rel, caLsky and da.acq pro- 
cesses and keeps the time information in the shared memory 


20 
(1) 


caljrel 


Acquires CDC data for all the channels corresponding 
to relative calibration 


21 
(1) 


cal_sky 


Acquires CDC data for all the channels corresponding 
to sky pedestal runs 


21 
(1) 


da_acq 


Acquires CDC data for all the channels corresponding 
to a Cherenkov event 


21 
(1) 


semaphore 


Regulates the access to a shared hardware 
amongst various peocesses 


17 

(3) 


abs-calib 


Acquires the absolute calibration data from 4 
specified channels and is part da-acq process 


21 
(1) 


pcr.ccr 


Acquires and displays the prompt and chance 
coincidence rates 


11 

(3) 


progjrncc 


Programs the RAM of the trigger generator modules 
as per user selected trigger configuration 


10 
(1) 


cdcjped 


Acquires the CDC pedestal data for all the channels 


10 

(1) 


sync-clk 


Synchronises the clocks of all the PCs on the network 
with the master node 


10 

(2) 



is also performed to evaluate the data quality and to look for the presence of a 7-ray signal 
in the data. More details regarding hardware and software features of the DAC system of the 
telescope are discussed in BOOI . 



Master process: 

The software is activated by the main process stdactic which runs on the master node with 
the highest priority and controls the execution of all other processes. On activation this pro- 
cess checks the network integrity and synchronizes the clocks of the other two PCs with the 
master node clock. Clock synchronization of the three nodes is important because after being 
spawned from the master node the various processes use the local PC clock for reference. 
Network integrity is checked by establishing a virtual circuit between processes on the local 
and remote nodes. The clock synchronization is followed by the spawning of a few other 
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Fig. 4.9: The pictorial representation of the 349-pixel imaging camera as created on the master node 
by the Gain Control Process, scrolling plots for PCR and CCR are also shown at the left of 
the pannel. 

basic processes like Semaphore, gain_ctl and daMcq on their respective nodes. A menubar 
of 16 icons each of which is associated with a specific task is also displayed on the master 
node. This master process communicates with other processes running on remote nodes to 
coordinate their execution by using messages and signals. The master process has a number 
of sequential interlocks built into it to ensure safe operation of the telescope. 

A semaphore has been developed under QNX to coordinate the functioning of multiple 
processes in the TACTIC software which are running on the same node and accessing the 
same CAMAC bus. This semaphore has been implemented using the IPC features of QNX. 
All processes accessing shared resources need to go through the semaphore to ensure that 
no other process is currently using the resource. In case of simultaneous access a semaphore 
list is prepared in which the highest priority process gets preference in execution. 



Gain control process: 

This process running on the Gain Control node controls the gain of the PMTs by changing 
the HVs applied to them on the basis of their SCR resulting from the ambient light-induced 
shot noise fluctuations f.95.1 . The process adjusts the HV of each of the PMTs to ensure that 
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its SCR is within the predefined range. The information about channel number, HV and SCR 
is sent to the remote master node for display. The pictorial representation of the 349-pixel 



imaging camera as depicted in Fig. 4.9 is also created by this process on the master node. 
This colour coded pictorial representation makes it convenient for the user to identify and 
switch off/on any malfunctioning pixel. An alarm is sounded whenever a pixel operates be- 
yond its safe limits and the HV to the camera is switched off whenever more than 30% of the 
pixels operate in an unsafe range. 

The gain control process apart from ensuring the operation of the PMT within safe anode 
current limits is also responsible for maintaining a stable value of SCR of individual pix- 
els despite being invariably exposed to different NSB light levels which can also include a 
bright star sometimes. These changes in the LONS background are further enhanced by the 
phenomenon of star field rotation [I282II which is inherent to altitude-azimuth mounted tele- 
scopes, like the TACTIC. This causes the image of an off-axis bright star in the field of view 
of the telescope to trace a circular path in the focal plane of the telescope thereby giving rise 
to large variations in the anode current and the SCR of concerned pixels over a time scale of 
tens of minutes. Since TACTIC imaging element hardware does not use any anode current 
monitoring ADCs, the same is measured in an indirect manner by keeping track of the SCR. 
If the SCR of a pixel under actual conditions is higher than the predesignated SCR range, 
the first attempt of the gain control process is to operate it at a lower HV to ensure that the 
PMT still remains ON, so that minimal information content is lost when the telescope gets 
triggered by an ACE. But if the recorded SCR is found to be more, indicating that the anode 
current of the PMT is still more than 30 fxA, the pixel is switched off. 



Data acquisition and control process: 

This process which runs on the data acquisition node controls the two processes on the local 
node which acquire event and calibration data. These two processes running with the same 
priority record the event data and the Relative Gain Calibration (RGC) data from all the 349 
channels along with time information on the occurrence of an event. A proxy is attached 
to each of these processes which is triggered by the interrupt handler. Operating in a mutu- 
ally exclusive mode the daMcq process records the ACE data while the caLrel and caLsky 
processes record the RGC and sky noise data, respectively. The Absolute Gain Calibration 
(AGC) data of the 4 designated channels is acquired by the abs.calib process which is also 
controlled by the da.acq process. This process also makes available the real-time event in- 
formation to all the processes which need it by using the shared memory. 
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Prompt and chance coincidence rates: 

The PCR and CCR processes running on the gain control node generate scrolling displays 
of the rates on the master node as shown in Fig. 4.9[ Any malfunctioning of the system is 



reflected immediately on these displays in terms of variations in these rates. This process 
provides facilities for changing online the integration time and the full scale count of the 
rate displays. These processes access the same CAMAC hardware as the gain control pro- 
cess through the Semaphore. The number of calibration runs taken by the system and the 
effective observation duration is also displayed on the console. 

Apart from the above listed main processes a number of stand alone routines for testing 
the various subsystems of the telescope are also available in the software. These routines are 
operated from the menu bar as and when required. 

Data filing system and data quality evaluation: 

In view of the large data volume recorded in an observation spell, the TACTIC data are stored 
on an hourly basis in a binary format on the master node. The system generates 15 files at 
the start up and adds 5 new files to it every hour. These 5 files which are referred to as cdc, 
cal, sky, scr and abs contain information regarding ACEs, RGC of the pixels, sky pedestal 
values, SCR values and AGC data, respectively. Apart from the AGC file which contains 
data for only 4 channels, all other files mentioned above have information about all the 349 
pixels used in the camera. These files are embedded with absolute time information with a 
resolution of 1 fis so that the data can be collated properly at the analysis stage. In addition 
to these files, several other files common to the whole observation spell are also generated 
by the system. These contain details regarding HV applied to the pixels, list of disabled pix- 
els and pixels operating at fixed HV, PCR, CCR and log of other system parameters. With 
each event consuming a space of ~ 708 bytes, the amount of data collected during a typical 
observation spell lasting for 5 hours (assuming an event rate of ~ 2.5 Hz with 7 calibration 
runs) turns out to be about 60 MB. 

Quasi on-line analysis of the data recorded by the telescope is also performed to evaluate 
the quality of the data recorded by the system and to identify hardware-related problems. The 
procedure involves transferring the hourly data files from the Master Node to Node4/Node5 
by using standard FTP protocol. Some of the checks which are routinely followed include 
analyzing the RGC and sky pedestal files so that pixels with either low or high gain can be 
identified and adjusted during the course of an observation spell. Checks like this have been 
found to be very useful for ensuring recording of good-quality data with minimum loss of 
observation time. The analysis code written in C-i-i- runs under the Linux operating system 
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and is equipped with elaborate data visualization tools using the interactive data analysis 
system ROOT fllO] , I109| . Apart from providing a graphical display of image parameter 
distributions, the software also contains suitable routines for extracting 7-ray signal from the 
recorded data by applying various event selection cuts. 

The data acquisition and control system of the TACTIC telescope has been working 
quite satisfactorily. The dead time of the system has been experimentally measured to be 
~ 2.5 ms by collecting the RGC data along with absolute time information at a trigger rate 
of ~ 400 Hz. After performing a time series analysis of the timing data it was found that 
any further increase in the event rate resulted in loss of some events by the data acquisition 
system. However, the loss of actual Cherenkov events corresponding to an event rate of 
~ 2.5 events s^^, by considering a non-paralizable dead time model turns out to be only 
0.63% which is safely acceptable. 

The gain control algorithm of the data acquisition and control system has been upgraded 
recently so that minimum possible time is spent on stabilizing the SCR of camera pixels. The 
modified approach being used now requires at least 3 pixels from the trigger zone to cross 
the cahbration limit for initiating a fresh calibration call. This has resulted in a substantial 
saving of the PMT gain control stabilization time and has allowed us to use about 80% of the 
actual observation time for collecting data on a candidate 7-ray source. 

4.8 Gain calibration scheme 

There are two schemes employed for calibrating the camera PMTs of the TACTIC telescope 
and are known as RGC and AGC. Currently the procedure followed by us for AGC, uses the 
same data which have been collected for the purpose of RGC. 

4.8.1 Relative gain calibration 

RGC scheme uses a high intensity blue LED placed at a distance of ~ 2 m from the camera 
surface to determine the relative gain of the camera pixels. The LED has been enclosed with 
a light-diffusing medium to ensure the uniformity of its photon field onto the imaging camera 
within ~ ±6%. The mean light intensity from the pulsed LED recorded by each pixel, in 
response to 2000 hght flashes is subsequently used for off-line RGC of the camera PMTs. 
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4.8.2 Absolute gain calibration 

In order to determine the AGC of PMTs, for the TACTIC telescope, two different procedures 
have been used: 

(/) The first procedure for determining AGC of the camera PMTs involves monitoring the ab- 
solute gain of a set of 4 gain calibrated pixels placed at the periphery of the camera as shown 



in Fig. 4.6 'a.). Since measurement of the absolute gains of these PMTs by determining their 
single photoelectron peaks a number of times during an observation run is rather time con- 
suming, we have instead used a relatively simpler method of measuring the light pulser yield 
of a calibrated source for the in-situ determination of the absolute gain of these calibration 
channels. The calibrated light sources used is the Am^''^ -based light pulsers (Scionix make) 
which produces fast optical flashes at an average rate of ~ 20 Hz with maximum emission 
at a wavelength of ~ 370 nm. The underlying principle for converting the charge content of 
an uncalibrated pixel from CDC counts to photoelectrons, uses the fact that the calibration 
pixels are also exposed to the light flashes from the LED during the relative calibration run 
and hence it becomes possible to obtain the conversion factors for all the remaining 345 pix- 
els of the camera [I283II . 

(//) The conversion for image size in CDC counts to number of photoelectrons has also been 
performed independently by using the excess noise factor method. Here we use the RGC 
data to estimate the conversion factor. After compensating for the gain normalisation be- 
tween PMTs, we calculate the variance in CDC counts between tubes for each RGC event 
and relate to the average number of photoelectron (pe). This relationship involves a larger 
variance than would be inferred from Poisson counting statistics alone due to fluctuations in 
the electron multiplication process [|97l . Considering Npe and cTpg are the average number 
of pe and correponding variance respectively at the photo cathode while Ncdc and acdc are 
measured average number of digital counts and measured variance in this count respectively. 
The following relations can be written between these quantities: 



a. 



pe 



(4.6) 



Ncdc = g CTpe (4.7) 

where g is the gain in units of digital counts per pe. The measured variance in digital counts 
can be written as 

C^cdc = g\/a y/Ope (4.8) 
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where a is a correction factor to take into account the additional fluctuations (caused due 
to electron multiplication process at the dynodes of the PMT). Using these relations the 
expression for the gain in terms of the measured quantities can be expressed as 

9--^ (4.9) 

The analysis of relative calibration data yields a value of Ipe = (6.5 ± 1.2) CDC for this 
conversion factor (g) when an average value of a ~ 1.7 is used for excess noise factor of the 
PMTs. Further work to determine the excess noise factor values more precisely for TACTIC 
PMTs is still underway. 
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5.1 Introduction 

lACTs do not directly detect the radiation emitted by the astrophysical sources being ob- 
served, but indirectly by monitoring the Cherenkov photons emitted by the secondary par- 
ticles (mostly electrons) produced in the EAS initiated by the incident 7-ray photon. The 
resulting ACR pulse is digitized by using the camera PMTs and associated electronics. The 
ACR images recorded in the camera are analyzed offline by subjecting them to various stan- 
dard analysis steps like noise cleaning, gain normalisation of PMTs and their parameter- 
ization. The statistical significance of the source detection is then evaluated by applying 
"gamma domain cuts" to the parameterized images (also called events). The energy es- 
timation of a primary 7-ray photon and subsequent determination of the observed source 
spectrum are the important steps in the data analysis chain. The data analysis software used 
to obtain results presented in this thesis has been developed on the LINUX plateform using C 
code integrated with the ROOT-library which is an object oriented data analysis framework 
developed at CERN [fTOll . 

Monte Carlo simulations have been performed to optimize gamma domain cuts, deter- 
mination of effective collection area and the quality factor of the ACT used for signal pro- 
cessing. In this chapter, we shall briefly discuss the simulation studies carried out for the 
TACTIC imaging element relevant to the present work and the detailed analysis procedure 
(starting with raw image to the spectrum determination), which have been followed. De- 
tails of the simulation studies for the TACTIC telescope like its energy threshold, trigger 
efficiency, optimization of imaging cuts and their energy dependence, estimation of the sen- 
sitivity can be found in ||203] , M , |266l , IM , [M] , IM. We begin with the discussion 
of Hillas parameters of an ACR image by describing its analytical formalism and physical 
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meaning. 



5.2 Parameterization of ACR images 
5.2.1 Hillas parameters 

The images of ACR observed by an individual lACT are typically analyzed by using moment 
analysis of the images to derive a set of quantities that characterize their roughly elliptical 
shape. These are referred to as the Hillas parameters after M. Hillas, who first proposed 
these parameters HI 8211 (see Fig |5.1 for their illustration). Moments are based on the CDC 



counts recorded in each pixel, together with its coordinates. Hence image parameters are re- 
constructed from the information contained in the pixels. Here we summarize the important 
image characteristics used in lACT image analysis 

• Size (S): The overall light content of the image is parameterized by the total number 
of pe contained in it. To the first approximation, the size is proportional to the energy 
of the primary particle. The size energy relation also depends on the impact parameter 
and and zenith angle of the observation. 

• Length (L): The RMS angular size along the major axis of the ellipse and related to 
the longitudinal development of the shower This parameter gives a measure of the 
parallax angle to the shower maximum and increases with the impact parameter 

• Width (W): The RMS angular size along the minor axis of the ellipse and related to 
the lateral development of the shower 

• Distance (D): The distance from the centroid of the image to the center of the field of 
view of the camera. This parameter gives information about the impact point of the 
shower i.e. a crude measure of core distance for a 7-ray shower It can also be viewed 
as the angle between the shower axis and the line joining the shower maximum and 
the telescope. In case of a stand alone telescope, this parameter is very important for 
energy estimation. 

• Miss (M): The perpendicular distance of the major axis of the image from the center 
of field of view of the camera. It is a measure of the shower orientation. 
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Fig. 5.1 : Illustration of image parameters. 

• Azwidth: The RMS angular size along a line, which is perpendicular to the line joining 
the image centroid to the center of field of view. It is a measure of both the shape and 
the orientation of the image. 

• Concentration or Frac2 (F2): Ratio of the two largest pixel signals to the sum of all 

signals. It represents the degree of light concentration and hence is a measure of the 
compactness of the image. 

• Alpha (a): The angle between the image axis ( i.e. shower axis) and the line joining 
the image centroid and the camera center (source position). It is basically the angle 
between shower axis and the optical axis of the telescope and hence is an orientation 
parameter. As events with small alpha point towards the source position, it is one of 
the most powerful parameters for the 7-hadron separation. The uncertainty in alpha is 
mainly due to the finite pixel size and the signal fluctuations (oc 1/ ^ size) [I290II . 
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• Asymmetry: It is a measure of the symmetric shape of the image. The 7-ray images 
should have tails which preferentially point away from the source position. 



5.2.2 Estimation of image parameters 

Estimation of the image parameters can be done on the basis of moment analysis. The zero- 
order moment is the size, i.e. the sum of CDC signals in all the PMTs (also called pixels) 
which are non-zero in a clean image. The first-order moments describe the position of the 
image in the focal plane of the camera and the second-order moments describe the extent of 
the image. Generally first and second order moments are sufficient for parameterization of 
the images. In this section, we summarize the formulae used to calculate the Hillas image 
parameters from the calibrated PMT signals [12961 . To calculate the image parameters, we 
first need to determine the image axis, which can be expressed as the equation of a straight 
fine with slope m and intercept c 

y = mx + c (5.1) 



The constants m and c are estimated in such a way that the resulting line minimizes the signal 
weighted sum of squares of perpendicular angular distances of the pixels i.e minimization. 
Considering Xj and yj denote the pixel coordinates of the i*^ pixel in degree with origin being 
the center of the camera and Si denotes the calibrated signal in the PMT, we can write the 
expression for as 



X 



E 

i=l 



Si iVi 



mxj 



1 + m? 



(5.2) 



where n is the total number of pixels. Assuming that both Xj and y^ have equal errors m 



and c can be obtained by differentiating the equation (5.2) The following quantities can be 
defined according to their standard definitions. 



<x>-- 



Si Xi 



En 
i=l 



<y>-- 



Era 
i=l 



(5.3) 



<x' ^-1^"^"^*') 



En 
i=l 



< xy >-- 



< y >-- 

En 
i=l 



En 
i=l 



(5.4) 



(5.5) 



(5.6) 
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< >= SklfifM < >= Sklfi^. (5.7) 

Spread of the images in different directions can then be defined in terms of moments as 

a^2 =< > - < X >^ (Ty2 ^< > - < y >^ (5.8) 



(Jxy ^< xy > - < X >< y > . (5.9) 
Some other higher order moments are 

a^3 =< > -3 < >< X > +2 < X >^ (5.10) 

Gyz ^< y^ > -3 < y^ >< y > +2 < y >^ (5.11) 

aj:2y =< x^y > -2 < xy >< X > - < x^ X y > +2 < X >^< y > (5.12) 



(7x2/2 =< xy'^ > -2 < xy >< X > - < y"^ >< x > +2 < y >^< x > . (5.13) 
Defining 

d — (Ty2 — (Ja;2 (5.14) 



z = ^Jd? + A{a,yy. (5.15) 
The constants of the image axis are then written as follows: 

m = (5.16) 

c^<y > -m <x > . (5.17) 
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The image parameters can be calculated from these moments and can be written as 



/(Jt.2 + (7„2 + Z 

length, L = (5.18) 



/ (T-r-2 + (7,,2 — Z 

width, W = J 2 ■ ^^-^^^ 

Since < x > and < y > are the coordinates of the image centroid, therefore 

distance, D = •>/< x >2 + < y >2 (5.20) 



miss, M = W -((1 + -) < X >2 +(1 - -) < y >2) ^ (5.21) 

2 z z z 



alpha, a — sin ^ ] (5.22) 



asymmetry — 



azwidth, Az = \/ (L^ sin'^a + cos'^a) (5.23) 



((Ta;3 cos^cf) + GyZ sin^cj) + 3ax2yCos'^4>sin4> + ^a^y^ coscj) sin^cj))^/^ 

L 

(5.24) 

J. J. • 7-ir> maa; max ,c r\c\ 

concentration bl— =^^^ . (5.25) 

Where sistmax and sundmax are maximum and second maximum signal respectively con- 
tained in pixels of an image. 

The image parameters described here are used for 7-hadron separation. Even after apply- 
ing 7-ray domain cuts (discussed later in this chapter) some fraction (~ 0.2%) of CR images 
is still classified as 7-rays due to the similarity of processes between shower development 
initiated by 7-ray and CR primaries. These CR events essentially become background for 
the 7-ray images. We have described here the a parameter generally used to extract the 7- 
ray signal from the background of CR images. This method of 7-hadron separation called 
"a approach" is source dependent i.e. it uses the information of the source position in the 
camera. Another source independent procedure called "^^ approach" can also be used to ex- 
tract the 7-ray signal. In this procedure source independent parameters e.g. length, width. 
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Fig. 5.2: Sketch demonstrating the Disp and 9 parameters. The Disp parameter leads to two esti- 
mated sky directions of an event (represented by solution A and B). The measurement of 
the asymmetry in the image along the major axis i.e. determining the head and tail of the 
image can provide the correct solution. Here the source is in the camera center and therefore 
solution A is the correct one. 6 is the angle between the estimated and real source positions. 
This figure has been adapted from |,129J . 



size etc. are used. While the "a approach" is useful for point like 7-ray sources, the anal- 
ysis of many extended galactic sources or on purpose off-source observation can be better 
performed by "6*^ approach". The later is often used for data taken in the Wobble mode 
(discussed later) of observation and uses the method called Disp analysis discussed here. 
However, in this thesis, we follow the "a approach". 



5.2.3 Disp procedure 

This procedure was first put forward by Fomin II152II and exploits the fact that the image 
elongation depends on the distance of the image centroid or Center of Gravity 

(CoG) of a 7-ray image from the source position in the camera plane. As depicted in Fig. |5.2[ 
the angular distance between the real and estimated source positions called 9 is obtained 
using the Disp method [I215L [T4T]| . The Disp parameter is defined as the angular distance 
between the expected source position and the CoG of the shower image along the major 
axis. In the Disp method the arrival direction of an EAS is estimated from the image shape 
parameters on an event-by-event basis assuming that the source (incoming direction of the 
primary particle) lies on the major axis of the image. This is plausible because the major axis 
is a projection of the shower axis onto the camera. Ideally the image axis coincides with the 
trajectory of a 7-ray photon II129I . 
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The elongation of the shower image and the value of the Disp parameter depend upon 
the impact parameter of the shower on the ground. For small impact parameter, the image 
should have a form close to that of a circle and be positioned very near the source position in 
the focal plane II152II . For increasing impact parmeter, it should become elongated and have 
a form close to that of an ellipse and be positioned farther from the source position in the 
focal plane. A simple form for the relationship between the elongation of an image and the 
Disp is llIBll 

/ Width \ 
Disp = d 1 - 7 r (5-26) 



Length J 

where ,^ is a scaling parameter. Generally, ^ depends on the height of the shower in the 
atmosphere, the zenith angle of observation, parameters of the model of the atmosphere and 
the energy of the primary particle. Experimentally, ^ is determined from the data so that the 
calculated shower arrival directions line up with a known source position. 

There is an important degeneracy in the Disp method: it provides two source positions 



along the major axis of the image (see Fig. 5.2 1 i.e. it does not provide information on which 
side of the shower image the arrival direction lies. Therefore a method to select the correct 
source position is needed. Images in the camera plane carry information about longitudinal 
development of the shower. The "asymmetry" charge distribution in the image contains 
the "head and tail" information of the recorded shower i.e. which edge of the image is 
closer to the source position in the camera plane. Cherenkov photons from the upper part 
of the shower produce a narrow section of the image with higher photon density (head) than 
the photons arriving from the shower tail which normally produce spread out end of the 
image [I141L 

The 9'^ distribution is expected to be flat for the background CR events and has an ex- 
ponential shape for a 7-ray signal peaking at 9'^ = 0° [I129II . By Using Disp method and a 
proper camera acceptance model a sky map of the arrival directions of the 7-ray events can 
be produced Il290l . 



5.3 Monte Carlo Simulations for TACTIC 

The development of an EAS in the atmosphere is a stochastic random process whose details 
can not be described by parametric functions or analytic expressions. Therefore Monte Carlo 
simulation are required to optimize the design and performance of an lACT. Further Monte 
Carlo simulations are helpful in understanding the reconstruction of 7-ray initiated showers 
and their segregation from much more abundant showers of hadronic origin. In order to 
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determine the effective collection area of the telescope, it is necessary to simulate its trigger 
efficiency. The energy reconstruction procedure also requires simulated 7-ray showers to 
estimate the energy of the primary particle. While 7-ray simulations agree quite well with 
the simple theory forming regular elliptical images, hadronic showers are dominated strongly 
by fluctuations in the initial nuclear interactions. 

5.3.1 CORSIKA 

While a number of Monte Carlo simulation codes e.g. CORSIKA (COsmic Ray Simula- 
tions for KAscade) Iml, MOCCA, KASCADE, ALTAI and GrISU ([[61 and references 
therein) are being used in VHE 7-ray astronomy, we have used the CORSIKA (version 
5.6211) air shower simulation code [I177H for predicting and optimizing the performance of 
the TACTIC imaging telescope. Originally developed to perform simulations for KASCADE 
experiment [I139II , CORSIKA is a detailed Monte Carlo code to study the evolution of EAS 
in the atmosphere initiated by protons, nuclei or any other particle. Its applications range 
from Cherenkov telescope experiments (~ 10^^ eV) up to the highest energies (~ 10^° eV) 
observed in EAS. The development of CORSIKA is guided by the idea to predict not only 
the correct average values of the observables but also to reproduce the correct fluctuations 
around the average value. The simulations of the electromagnetic part in an air shower in 
CORSIKA are treated with the code called EGS4, while hadronic interactions are simu- 
lated by models, GHEISHA for low energies (< 10 GeV) and VENUS, QGSJET, DPMJET, 
SIBYLL, etc. for higher energies (> 10 GeV) ( 1117711 and references therein). 

5.3.2 TACTIC simulation chain 

Generally, the complete execution of Monte Carlo simulations for an lACT is sub-divided 
into two parts. The first part comprises generation of air showers induced by different pri- 
maries and recording the relevant raw Cherenkov data (data-base generation), while folding 
in the light collector characteristics and PMT detector response is performed in the second 
part. In order to carry out the simulation study for the TACTIC telescope, we have used the 
altitude of the Mount Abu 7-ray observatory i.e. 1300 m, corresponding to air mass gram- 
mage ~ 890 g cm'"^ and the horizontal and vertical magnetic field component values of 
35.9 jiTesla and 26.6 jiTesla respectively [I281L The Cherenkov photon wavelength band 
of 300 — 450 nm have been chosen for generating the data-base. The simulated data-base for 
7-ray showers uses about 34000 showers in the energy range 0.2 — 20TeV with an impact 
parameter of 5 — 250 m. These showers have been generated at 5 different zenith angles 
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(e = 5°, 15°, 25°, 35° and 45°). A data-base of about 39000 proton initiated showers in the 
energy range 0.4 — 40 TeV, have been used for studying the 7-hadron separation capability 
of the telescope. The representative CR integral energy spectrum of the form E^^-"^ has been 
used in the data-base generation. The incidence angle of the proton showers is simulated 
by randomizing the shower directions in a field of view of 6° x 6° around the pointing di- 
rection of the telescope. A photon bunch received by a detector at the observatory level is 
characterized using the CORSIKA code by seven parameters, namely bunch size, x and y 
coordinates, two direction cosines, height of production and arrival time at the ground. The 
Cherenkov photons are ray-traced [I284II to the detector focal plane and the number of pe 
likely to be registered in a PMT pixel are inferred after folding in the atmospheric absorption 
of the photons, optical characteristics of the mirrors, the metallic compound-paraboloid light 
concentrators at the entrance window of the pixels and the photocathode spectral response. 
The Cherenkov photon data-base, consisting of number of pe registered by each pixel, is 
then subjected to noise injection, trigger condition check and image cleaning. 

The Cherenkov pulses are detected in the presence of LONS in actual experimental con- 
ditions. To estimate this contribution, measurements on a few PMTs are carried out in the 
laboratory with the help of a steady light source of a similar intensity as that of LONS at 
Mount Abu. The resulting shot noise is then superimposed on the ACR generated pe con- 
tent of each pixel. To revert it back i.e. to exclude the shot noise contribution generated 
by LONS, a standard two level "tail cut" procedure has been used to obtain the clean im- 
age H 14911 . The aim of this procedure (i.e. image cleaning) is to retain maximum pixels 
containing Cherenkov light while suppresing the pixels with signal from LONS alone. The 
procedure involves the following steps: 

(i) Pixels with less than 3. Oct signal (so called "boundary threshold") are set to zero. 

(ii) Pixels with more than 6.5a signal (known as "picture threshold") are retained. 

(iii) Pixels with signal in between 3.0cr and 6.5cr are retained provided they are adjacent to 
a pixel satisfying condition (2). 

The resulting data-bases, consisting of pe distribution in the imaging camera at various core 
distances and zenith angles are then used for estimating trigger efficiency, effective detec- 
tion area, optimum ranges of Cherenkov image parameters for discriminating between 7- 
ray and CR events, differential count rate and effective threshold energy of the telescope for 
7-ray and CR proton events and finally developing the energy reconstruction procedure for 
7-ray events. The clean Cherenkov images are characterized by calculating their standard 
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image parameters like length, width, distance, alpha, size and Frac2 as described in sec- 
tion 5.2.2 The standard Dynamic Supercuts procedure [I232II is used to separate 7-ray like 
images from the background CR. 



5.3.3 Procedure for 7-hadron separation 

ACT is quite efficient in 7-hadron separation above 100 GeV. The technique exploits the 
characteristic differences between the shower development processes initiated by these two 
primary particles and consequently the Cherenkov image generated at the telescope camera. 
Based on these differences several strategies have been developed for 7-hadron separation, 
which include supercut procedure, extended or dynamic supercut procedure and cluster anal- 
ysis [I232II . The sensitivity of an lACT strongly depends on the rejection capabilities of CR 
background events and can be improved by considering the energy and zenith angle de- 
pendence of image parameters. A comparative study has been carried out for the TACTIC 
telescope to evaluate the performance of supercuts, dynamic supercuts and zenith angle de- 
pendent dynamic supercuts. The details of this study can be found in [I202II . 

The performance of different 7-ray selection methodologies is quantified by the quality 
factor (Q) for the telescope, which is basically the ratio of fractions of 7-ray and CR retained 
events after the selection procedure is applied. In the following we summarize the supercut 
procedure since it is the basic methodology and then the dynamic supercut procedure which 
accounts for the energy dependence of image parameters. 



Supercut procedure: 7-ray images have length, width and alpha parameters, which are 
restricted to a narrow range of values as compared to proton images. This can be used to 
optimize the range of parameters to maximize the significance of the source detection. We 
have optimized the ranges of image parameters for the simulated data-base that results in 
maximum Q value i.e. the maximum acceptance of 7-rays with minimum contamination 
of background CR proton events. The range of these values for 7-ray domain are listed 



in Table 5.1 Note that the size cut is applied to eliminate events dominated by noise, while 

Table 5.1 : 7-ray domain supercuts. 



size 


S > 50 pe 


distance 


0.5° <D< 1.2° 


length 


0.11° <L< 0.33° 


width 


0.05° <W < 0.17° 


alpha 


a < 18° 
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distance parameter (rather poor discriminator) is restricted to maintain the energy resolution 
for the telescope [12811 . 

Apart from resulting in comparatively lower quahty factor, there are other disadvantages 
in using this approach directly for the extraction of spectra because length and width of the 
images increases with the energy of the primary VHE 7-rays , and this exactly is taken care 
of in dynamic supercut methodology. 



Dynamic supercut procedure: A straight forward method to incorporate the energy de- 
pendence in selection procedure is to scale the upper and lower values for each parameter 
with the energy of the primary 7-ray photon, which of course is unknown. However the 
image size is a fairly good estimate of energy and length and width are also well correlated 
with the log{size). The first order polynomial fit of these correlations {length v/s log(size) 
and width v/s log{size)) can be used to obtain the dynamic supercut ranges of parameters. 
The upper values of these parameters for the TACTIC telescope have been optimized for 
higher 7-ray acceptance compared to the supercuts, while lower values have been kept fixed 
because of insignificant effect. The ranges of dynamic supercuts for VHE 7-rays are given 
in Table m 



Table 5.2: Range of 7-ray domain dynamic supercuts. 



size 


5 > 50 pe 


distance 


0.5° <D < 1.2° 


length 


0.11° <L< (0.27527 + 0.023 ln{S)) 


width 


0.05° <W < (0.06356 + 0.026 ln{S)) 


alpha 


a < 18° 



The 7-ray retention factor as a function of 7-ray energy has been compared in Fig j5.3| 
for supercut and dynamic supercut procedures. Comparatively better 7-ray retention for 
dynamic supercuts is evident. We have used the dynamic supercut procedure in the work 
presented in this thesis. 



5.3.4 Comparison with real data 

The agreement between the expected and actual performance of the telescope can be checked 
by comparing the expected and observed image parameter distributions. Fig ]5.4 shows the 
distributions of the image parameters length, width, distance and alpha [11821 , [296[[ for 
simulated protons and for the actual Cherenkov images recorded by the telescope. The sim- 
ulated distributions of these image parameters for 7-rays have also been shown in the figure 
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for comparison. The observed image parameter distributions are found to closely match the 
distributions obtained from the simulations for proton-initiated showers, testifying to the fact 
that the event triggers are dominated by background CRs. 



5.4 Evaluation of telescope parameters 

In order to obtain the observed 7-ray flux from the detected number of 7-rays , we need to 
use several characteristic parameters of the telescope e.g trigger efficiency, effective area, 
threshold energy and CR detection rate, etc. These characteristic parameters in turn depend 
on design features of the telescope i.e. size of the light collector, pixel size, field of view of 
the camera, trigger generation scheme, single pixel threshold level and the image parameter 
cuts. We describe here some of these important parameters. 

5.4.1 Single pixel threshold and topological trigger 

The NSB strongly limits the energy threshold of an lACT since the trigger threshold is nor- 
mally set to the LONS fluctuations. However, LONS contribution can be reduced by em- 
ploying the smaller pixel size and limiting the charge integration time. The variation of SCR 
as a function of its trigger threshold can be used to select the appropriate trigger scheme, 
which also governs the trigger threshold energy of the telescope. Further, the variation of 
SCR with threshold level can be different from a pure Possionian distribution because of 
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Fig. 5.4: Comparison of observed and simulated image parameter distributions. 



the afterpulsing in PMTs [I231II . The measured variation of SCR with trigger threshold for a 



PMT is shown in Fig 5.5 Given the large number of pixels forming the trigger region, the 
telescope trigger rate turns out to be extremely high. In order to suppress these unwanted 
random triggers, we use some topological schemes for trigger generation. Several topolog- 
ical trigger configurations e.g. NNP, 3NCT and ^^^6*2 have been worked out and resulting 



CCR have been plotted for comparison in Fig 5.5 

By fixing up the tolerable CCR (say 10~^ Hz as shown in the figure) one can compare 
the single pixel thresholds for different trigger schemes. For example, in NNP trigger con- 
figuration the single pixel threshold can be set at ~ 25 pe while 3NCT configuration allows 
us to go down to ~ 8 j9e level for the same value of CCR. It is important to mention here that 
lowering the single pixel threshold can in principle be compensated by increasing the trigger 
multiplicity. However, the optimum value of single pixel threshold which finally decides the 
threshold energy of the telescope is strongly dependent on the intensity of Cherenkov pho- 
tons in the image plane and the pixel size used 028111 . Both NNP and 3NCT trigger schemes 
have been used for the TACTIC telescope from time to time, here we shall use NNP trigger 
scheme for the evaluation of other telescope parameters. 
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Fig. 5.5: Chance coincidence rates for different topological trigger schemes as a function of single 
pixel threshold. The discrimination level is converted into pe after obtaining the single 
photoelectron spectrum for the PMT. 



5.4.2 Trigger efficiency 

Trigger efficiency r]{r, E), is the probability that the telescope gets triggered by a 7-ray or 
proton shower of primary energy E, arriving at a distance r from the telescope. In simula- 
tions, it is defined as 

( E) showers triggering the telescope 

' total number of showers generated 

Apart from core distance and energy of the primary particle, trigger efficiency also depends 
on the nature of primary particle and hence on the Cherenkov photon density, trigger field of 
view of the camera, the single pixel threshold level and trigger generation scheme. The trig- 
ger efficiencies for the TACTIC imaging telescope have been evaluated for 3NCT and NNP 
trigger schemes at 8 pe and 25 pe threshold levels respectively, both for 7-ray and proton 
generated Cherenkov events. 
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5.4.3 Effective collection area 

Cherenkov telescopes have an effective collection area determined by the area of the ACR 
pool at the ground which has a typical radius of ~ 120 m. Such a huge effective collection 
area makes the Cherenkov telescopes highly efficient for the detection of VHE 7-rays where 
flux becomes significantly low due to the power law nature of most of the emissions from as- 
trophysical sources. The energy dependent effective collection area Af,ff{E) can be written 
in a general form as 

POO 

A^ff{E) =271 r X ri{r, E) dr. (5.27) 
Jo 

This integration is normally discretized in simulations and for the TACTIC imaging telescope 
the effective collection area corresponding to each primary energy value is evaluated from 
the following equation 

^.//W = ^E('-L-'-?) "'''"'/"'''-' (5.28) 

1=1 

where n = 20, ri = 5 m and r2i = 205 m. The evaluated effective collection area for 
the TACTIC imaging telescope has been shown in Fig j5.6| (a) and (b) for 7-ray and proton 
induced ACEs respectively, at two representative zenith angle values of 15° and 35°. These 
results were obtained by using the NNP trigger scheme with 11x11 pixels trigger field and 
a single pixel threshold of > 25 pe. While significant variation can arise in the calculation of 
effective area at higher primary energies, the contribution of such events is small due to the 
power law nature of 7-ray spectrum of the VHE sources. 

The effective collection area for the isotropic CR background can be evaluated from the 
following equation 

/•oo /"OO 

Aeff{E) = 2n / r X r]{r,E,n)drdn (5.29) 
Jo Jo 

where dil is the solid angle subtended by a detector element. 

5.4.4 Telescope trigger rates and threshold energies 

The differential trigger rate (^) of a telescope is defined as the number of particles with 
energies between E and E + dE, triggering the telescope per unit time. In simulations, the 
differential trigger rate for a particle with known differential spectrum (^), is obtained by 
directly multiplying the differential spectrum with the energy dependent effective collection 
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Energy {TeV) Energy (TeV) 



Fig. 5.6: Effective collection area of the telescope for 7-ray (a) and protons (b) as a function of the 
primary energy at zenith angles of 15° and 35°. Differential trigger rates for 7-ray (c) and 
proton events (d) as a function of the primary energy. These results have been obtained by 
using the NNP trigger scheme with 11x11 pixels trigger field and a single pixel threshold 

of > 25pe. 



area (^e//)- Therefore the differential trigger rate can be written as 

IE = " lE- 

The peak of the differential trigger rate curve as a function of energy determines the threshold 
energy of the telescope for a given type of primary particle. The integral trigger rate is then 



obtained by simply integrating equation (5.30) 



R{E>Eth)= / -=dE (5.31) 

where Eth is the energy threshold of the telescope for a given type of particle and Emax 
is the maximum energy of the primary particle that can be detected by the telescope. The 
differential energy spectra used in the TACTIC simulations in the energy range 0.3 — 50 TeV 
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has the following form for 7-ray PTI and CR protons [|38ll 

^ = 2.79 X 10-^^^2.59 photons m'^ s"^ TeV'^ (5.32) 
dE 

^ = 1.10 X IQ-^E-'^-'^^ particles m'^ sr'^ s'^ TeV'^ (5.33) 
dE 

respectively. Fig j5.6| (c) and (d) show the differential event rates as a function of the pri- 
mary energy for 7-ray and proton events, respectively. It is evident from these plots that 
the 7-ray trigger threshold energy of the telescope is ~ 1 TeV. The corresponding trigger 
threshold energy of the telescope for protons turns out to be ~ 2.5 TeV. The reason for 
this difference in threshold energies is the higher Cherenkov photon density (almost twice) 
produced by 7-ray photon EAS as compared to that of proton initiated EAS. Additionally, 
proton images are more diffuse, and are less likely to reach the trigger threshold. At higher 
energies (> 20 TeV), the photon density is sufficiently high for both 7-rays and protons to 
efficiently trigger the telescope and the effective area is mainly governed by the radius of 
the ACR pool. The effective area for protons, at still higher energies (> 50 TeV) increases 



as compared to 7-rays , (see Fig 5.6 (b)), due to the triggering of the telescope by proton 
showers, whose directions fall even outside the trigger field of view of the camera. 

Showers arriving at large zenith angles develop their shower maximum at relatively 
greater heights with respect to the telescope. This dependence goes approximately as Height 
oc l/cos{9), where 9 is the zenith angle of the telescope. This allows the showers arriving 
at larger core distances to remain within the trigger field of view of the camera. Thereby 
geometrically increasing the effective collection area at large zenith angles. Again, because 
of the geometry, the projection of ACR pool on the ground for a shower arriving at large 
zenith angle will be larger. The dependence of radius of ACR pool also goes approxi- 
mately as oc I /cos (9) and consequently the Cherenkov photon density goes down as nearly 
oc l/cos'^{9). At higher energies, the effective collection area increases with increasing 
zenith angle purely due to the geometrical effect as long as the photon density is sufficient 
to trigger the telescope. At lower energies, the trigger efficiency decreases nearly the same 
way as the Cherenkov photon density and results in decrease in effective collection area 
with increasing zenith angle. Correspondingly, the threshold energy of the telescope also 
increases at higher zenith angles. These effects have been taken into account while obtaining 
the energy spectrum of a 7-ray source. 

The agreement between the predictions from Monte Carlo simulations and the actual 
performance of the telescope have been checked by comparing the observed trigger rate 
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of the telescope with the predicted value. The expected prompt coincidence rate at a zenith 
angle of 15° turns out to be ~ 2.5 Hz for the NNP trigger mode. This value has been obtained 



on the basis of integrating the differential rate curve for protons Fig 5.6 (d). Reasonably 
good matching of this with the experimentally observed value of ~ 2 — 3 Hz suggests that 
the response of the telescope is very close to that predicted by simulations. 

5.4.5 Sensitivity 

As already discussed in chapter 3, the sensitivity of a 7-ray telescope is its ability to detect 
a 7-ray signal over the CR background. The sensitivity of an lACT can be expressed in two 
ways: 

(i) If the flux of the source is known (e.g. Crab Nebula like sources), we estimate the 
sensitivity in terms of the minimum time required to detect the 7-ray signal, which 
is statistically significant (usually taken as 5a). 

(ii) When the flux of the source is unknown, the source is observed for a fixed time period 
(usually 50 hours) and sensitivity is expressed as the minimum 7-ray flux required 
to detect the signal at 5cr statistical significance level. 

We shall follow the first approach here. The number of 7-ray events detected by an lACT, 
which survive 7-ray domain parameter cuts can be written as 



N,= -j^A,{E)r^,{E)TdE (5.34) 

where i]^/{E) is 7-ray retention factor and the integration is over the dynamic energy range of 
the TACTIC telescope i.e. ~ 1 .0 TeV to ~ 20.0 TeV. The lower limit of the dynamic energy 
range signifies the threshold energy while the upper limit corresponds to the energy at which 
the charge content of the ACE is sufficiently high to saturate the CDC counts. The number 
of detected events, generated by protons and again surviving after 7-ray domain parameter 
cuts, can similarly be written as 

Np= -^Ap{E)7]p{E)TndE (5.35) 
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where rip{E) is proton retention factor, Q is the solid angle of the telescope and T is the 
observation time. The statistical significance (N^) of a 7-ray signal detection is defined as 



- N, 



off 



// 



(5.36) 



where Non = Ny + Np and iVo// = Np are the number of events detected in on- source and 
off- source regions. The above equations can be solved for the minimum time required to 
obtain A^^r significance, which can be expressed as 



T 



Nt 



1 



+ 2VL- 



L 



Err 
Em 



dE 



fAp{E)7]p{E) dE 



jE:::"-^ME)v,iE)dE ' "'\j^-'-^A,iE)rj,{E)dE)\ 



(5.37) 



After substituting the appropriate values for i]^, rjp, A^, Ap, n and flux values for 7-ray and 



protons from equations (5.32) and (5.33) respectively, the minimum time required to de- 



tect a 7-ray signal from the Crab Nebula (flux is given by equation ( |5.32| )), at a statistical 
significance of 5cr turns out to be ~ 25 hours. 

This sensitivity figure has been confirmed repeatedly by analyzing the data collected 
on the Crab Nebula. Data collected on the Crab Nebula for ~101. 44 /lowrs between 10 
November 2005 - 30 January 2006 have yielded an excess of ~(839±89) 7-ray events with 
a statistical significance of ~9.64cr. 



5.4.6 Energy reconstruction procedure for 7-ray events 

Given the inherent power of Artificial Neural Network (ANN) to effectively handle the mul- 
tivariate data fitting, an ANN-based energy estimation procedure has been developed for 
determining the energy spectrum of a candidate 7-ray source. An ANN is a non-linear statis- 
tical data modeling tool which can be used to model complex relationships between inputs 
and outputs. The tool can be applied to problems like function approximation or regression 
analysis (similar to what is being attempted here), classification (pattern recognition) etc. 
Although not for energy estimation, the idea of applying ANN to imaging telescope data 
was attempted for the first time by Reynolds and Fegan II262L While a detailed description 
of ANN methodology can be found in [12621 . 12131 . we present here only the main steps which 
have been followed to ensure reliability of end results. 

The 7-ray energy reconstruction with a single imaging telescope, in general, is a func- 
tion of image size, distance and zenith angle. The procedure followed by us uses a 3:30:1 
(i.e. 3 nodes in the input layer, 30 nodes in hidden layer and 1 node in the output layer) con- 
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figuration of the ANN with resilient back propagation training algorithm [I260II to estimate 
the energy of a 7-ray like event on the basis of its imagesize, distance and zenith angle. 
The training and testing of the ANN was done in accordance with the standard procedure of 
dividing the data base into two parts so that one part could be used for the training and the 
remaining for testing. We used about 10,000 events out of a total of 34000 parameterized 
7-ray images for training. The 3 nodes in the input layer correspond to zenith angle, size 
and distance, while the 1 node in the output layer represents the expected energy (in TeV) 
of the event. In order to make ANN training easier and to smoothen inherent event to event 
fluctuations, we first calculated the < size > and < distance > of the training data sample 
of 10000 images by clubbing together showers of a particular energy in various core dis- 
tance bins with each bin having a size of 40 m. Once satisfactory training of the ANN was 
achieved, the corresponding ANN generated weight-file was then used as a part of the main 
data analysis program so that the energy of a 7-ray like event could be predicted without 
using the ANN software package. Rigorous checks were also performed to ensure that the 
network did not become "over- trained" and the configuration used was properly optimized 
with respect to number of iterations and number of nodes in the hidden layer. Optimiza- 
tion of the network was done by monitoring the RMS error while training the ANN. The 
optimized configuration (3:20:1 with 5000 iterations) yielded a final RMS error of ~ 0.027 
which reduced only marginally when the number of nodes in the hidden layer or the number 
of iterations were increased further. 

A plot of energy reconstruction error obtained for a test data sample of 24000 parame- 



terized images is shown in Fig. 5.7 a). This plot has been obtained by testing the ANN with 
individual parameterized images at different energies within 5° — 45° zenith angle range. 
Denoting the estimated energy by Eann and the actual energy used during Monte Carlo 
simulations by Emc^ Fig- |5^b) shows the frequency distribution of \n(EMcl^ANN) for all 
events shown in Fig. 5.7 ^a) alongwith a Gaussian fit. Having a value of cr(ln E) ~ 28.4% 



Fig. 5.7 ;b), directly implies that the procedure should allow us to retain ~ 84% 7-ray s in 
most of the energy bins from 1 — 20 TeV, if 6 bins per energy decade i.e. cr(ln E) of~ 40% 
are used for determining the energy spectrum of a candidate 7-ray source. Note that the 
energy resolution of ~ 40% has been used here to avoid spill over of events in between the 
neighbouring energy bins. The performance of the ANN-based energy reconstruction proce- 
dure was also compared with the results obtained from the linear least square fitting method. 
This method yielded a a(\n E) of ~ 35.4%. It is also worth mentioning here that the new 
ANN-based energy reconstruction method used here, apart from yielding a lower a(\n E) of 
~ 28.4% as compared to o"(ln E) of ~ 36% reported by the Whipple group 1123211 . has the 
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Fig. 5.7: (a) Error in the energy resolution function Xvi^EmcIEann) as a function of the energy of 
the 7-rays for the simulated showers generated within 5° — 45° zenith angle range, (b) 
Frequency distribution of \n{EMclEANN) for all events shown in Fig. 5.7 a) along with a 
Gaussian fit. 



added advantage that it considers zenith angle dependence of size and distance parameters 
as well. The procedure thus allows data collection over a much wider zenith angle range as 
against a coverage of upto 35° in case the zenith angle dependence is to be ignored. 

5.5 TACTIC data analysis 

Presently, there are two primary goals in analyzing the data recorded with the TACTIC tele- 
scope: 

(i) 7-hadron separation that allows to detect a 7-ray source. 

(ii) Estimation of primary 7-ray energy of the 7-like events, to derive an energy spectrum 
of a detected 7-ray source. 
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Fig. 5.8: Flow chart of data analysis procedure followed for 7-ray signal extraction and energy spec- 
trum determination. 



Detailed analysis of the telescope data involves several steps. Before applying the 7-ray do- 
main parameter cuts, the raw Cherenkov images recorded by the telescope need to be pro- 
cessed through the data quality checks, pedestal subtraction, noise cleaning and PMT gain 
normalization (commonly called flat fielding). All the main steps involved in the analysis 



have been shown by means of a flow chart in Fig. 5.8 Here, we describe the detailed analy- 
sis procedure which has been used in analyzing the data recorded with the TACTIC telescope 
on different TeV 7-ray sources. 



5.5.1 Data quality checks 

The quality of data recorded with the TACTIC telescope depends on a number of parameters 
hke sky conditions, dust levels, instrumentation malfunctions etc. In order to ensure that the 
quality of recorded data is of a high order, number of data quality checks have been used 
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Fig. 5.9: Distribution of arrival time differences of successive events, recorded with the TACTIC 
telescope. Red curve shows an exponential fit to the experimental data yielding a mean rate 
of ~ 2.5 Hz. Here data used have been recorded in the zenith angle range of 10° to 35° 
from the Crab Nebula direction. 



which includes the compatibility of event rate with the Poissonian statistics and its expected 
dependence on zenith angle. Some of these are described here: 

• A visual log of sky conditions is recorded through out the observations and referred 
to while handling the data. 

• CRs reaching the Earth isotropically are random both in direction and timing. It is 
therefore expected that the frequency distribution of the arrival time differences of 

these events should be exponential with the exponent determining the mean event rate 
detected by the telescope. This frequency distribution for one spell of recorded data is 



shown in Fig. 5.9 and is indeed exponential. The exponent obtained by fitting the data 
is in agreement with the observed rate of ~ 2 — 3 Hz for good spells, as confirmed 
with visual inspection of sky conditions. 
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Fig. 5.1 0: Variation of prompt coincidence rate as a function of zenith angle of observation. A fit to 
the data identifies (a) the good data run and (b) the bad data run. 



• The main data quality check, which indicates the possibility of problems during ob- 
servations is the event trigger rate. As discussed earlier, the effective area increases 
with increasing zenith angle while the photon density decreases and as a net effect the 
telescope trigger rate decreases with increasing zenith angle, following the rela- 
tionship 

Re^ Roicose)"^ (5.38) 

where 9 is the zenith angle of observation, Rq is the event rate at zenith angle 9 and 
Rq is at 6* = 0°. The value of m depends upon the sky transparency and we have used 
m between 1-3 for accepting good quality data. 
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Fig. 5.11: Distribution of sky counts for pixel no. 60. 



A significantly lower trigger rate than the expected one, generally indicates bad atmospheric 
conditions or hardware malfunctioning, while higher rates that exceed the allowed statistical 
fluctuations points to the presence of bright stars in the field of view or even could again 
be a hardware related problem, like the instability in some of the discriminator channels or 
malfunctioning of trigger generator module. A representative example of selecting good data 
spells based on this procedure has been shown in Fig. 5.10[ wherein a fit to the observed data 
based on equation ( |5.38 1 qualifies the spell shown in Fig. 5.10(a) as good quality data, while 
Fig. 5.10(b) shows a bad spell, where rate increases with increasing zenith angle and the 
value of m obtained is negative and falls beyond its allowed values for good sky. 

Onhne tests have also been performed, which demonstrate the over all status of the chan- 
nels including PMTs and associated electronics. These tests are carried out by analyzing the 
calibration data collected several times during a typical data taking run, immediately after 
their acquisition. While results of these tests are demonstrated in Figures 5.1 1[ 5.12[ 5.14 
and|5.15[ we discuss them in the following sections. 



5.5.2 Pedestal subtraction 

The pedestal of a CDC channel is the value of it's charge content converted into CDC counts 
(ICDC count = 0.125pc) which it outputs in the absence of any input signal. This is nor- 
mally set around 100 CDC counts so that small positive polarity fluctuations on the signal 
line due to LONS induced shot noise fluctuations, do not generate erroneous value of CDC 
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Fig. 5.12: Sky pedestal run, taken during an actual observation run. Plot shows the mean sky counts 
for 256 pixels with error bars of 1 cr statistical fluctuation from 2000 samples. 



counts. The pedestal value for each PMT channel is obtained by artificially triggering the 
camera and recording CDC counts. The trigger rate is kept quite high (~ 400 Hz), so that 
the probability of recording an ACE while recording pedestal data is very low. Acquiring 
the calibration data in this artificial trigger mode is called the sky pedestal run. A typical sky 
pedestal run consists of 2000 events. The mean value of CDC counts and its variance over 
2000 events is then calculated and used for pedestal subtraction and noise cleaning respec- 



tively. Fig. 5.11 shows the distribution of CDC counts for pixel number 60, while the mean 
value and associated standard deviation for all the pixels have been shown in Fig. 5.12[ The 
pedestal subtraction process involves the subtraction of the mean values of sky pedestal CDC 
counts for all the pixels from the total counts of corresponding pixels recorded in presence of 
the ACR pulses. Clearly, the aim is to subtract the mean value of LONS, which contributed 
in the signal while digitizing the ACR events. 



5.5.3 Image cleaning 

Image cleaning is one of the most crucial step in the analysis chain to be followed for an 
lACT, as the subsequent parameterization of the images uses only the signal contents of the 
pixels, without its corresponding errors. Most of the pixels in the camera contain only fluc- 
tuations due to LONS as indicated in the previous section, while the information about the 
air shower which triggered the readout is confined to smaller number of pixels. Therefore, in 
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(a) Raw image (b) Clean image 

Fig. 5.13: Example of an image of a typical shower recorded with the TACTIC telescope. Color 
code is an indication of the charge content of the pixel, (a) the raw image i.e. prior to the 
processing procedure, and (b) the same image after pedestal subtraction and application 
of picture and boundary thresholds. 



order to reconstruct the image of an EAS, the contribution of the pixels containing the LONS 
fluctuations alone have to be suppressed. Moreover, the image is discarded if any pixel in the 
image is saturated either due to presence of a bright star or even due to genuine Cherenkov 
signal, since the corresponding energy of the primary particle can not be estimated properly. 

The purpose and procedure of image cleaning have already been discussed in section 5.3.2[ 
we summarize the standard image cleaning procedure again here. A pixel is selected as part 
of the image if it has a signal above a certain threshold called picture threshold or is beside 
such a pixel and has signal above a lower threshold called boundary threshold. The pic- 
ture and boundary thresholds are multiples of the RMS sky pedestal deviation which PMT's 
signal must exceed to be considered part of the picture or boundary, respectively I1216II . 
The RMS deviations of the pixels are evaluated by using a sky pedestal run and shown in 



Fig. 5.12 The picture and boundary pixels together make up the image, while all others are 



set to zero. This image cleaning procedure is depicted in Fig. 5.13 A picture threshold of 
6.5 cr and a boundary threshold of 3.0 a have been used to select the largest number of pixels 
while at the same time limiting the inclusion of pixels with noise alone. 
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Fig. 5.14: Distribution of LED counts for pixel no. 60. 

5.5.4 Gain normalization 

This process is also called "flat fielding". In this process, we account for the differences 
in the relative gains of the PMTs. Operating in the pulsed mode, a blue LED based light 
pulser illuminates the camera uniformally (±5%) from a distance of ~ 2 m from the camera 
surface. In response to the LED pulse, the readout is triggered and corresponding CDC 
counts are recorded for all the pixels. A total number of 2000 such calibration pulses are 
recorded, by triggering the camera at a rate of ~ 400 Hz. The calibration data collected 
in this mode is called the relative calibration run. The mean and standard deviation is then 
calculated for each pixel on the basis of these 2000 events. The frequency distribution of 



CDC counts, for pixel number 60 is shown in Fig. 5.14 for a typical relative calibration 
run. The mean values and associated variances for 256 pixels are depicted in Fig. 5.15 The 
calibration data collected are used to determine the gain normalisation factors of each PMT 
by comparing their mean CDC counts to the reference pixel. Nine pixels at the center of the 
camera are kept at fixed HV during the course of observation and any one of them can be 
used as a reference pixel. 

A sky pedestal run and a relative calibration run are together referred to as a calibration 
run. In order to achieve the minimum possible trigger threshold, we allow the control system 
of TACTIC telescope to optimize the HV and therefore the gain of the PMTs. The values 
of HV are changed only for those PMTs whose SCRs have changed beyond the allowed 
SCR range in responce to the NSB light. This also ensures the safe anode current operation 
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Fig. 5.15: Relative calibration run, taken during an actual observation run. Plot shows the mean LED 
counts for 256 pixels with error bars of 1 cr statistical fluctuation from 2000 samples. 

of PMTs. We need to recalibrate the camera after each such change in the HV values. By 
examining the mean values and standard deviation of the charge content of the pixels, the 
possible problems in the hardware can be identified since these were set in an allowed range 
for a particular setting of the LED pulse. Also by comparing them from one calibration 
run to another, the online tests for overall channel health can be performed as mentioned in 
section U. 5. 1[ 



5.5.5 7-hadron separation 

The 7-hadron separation is based on image parameters obtained using the Hillas image pa- 
rameterization methodology HI 821 of the cleaned images which has been discussed in sec- 



tion 5.2.2 The derived image parameters {length, width, size, distance, alpha, frac2 etc.) 
are then used (except a, which is used separately) to segregate 7-ray events from the over- 
whelming background of hadron events by applying the 7-domain cuts, which are obtained 
from the Monte Carlo simulations. The representative examples of a 7-ray like image and 



a proton image as recorded with the TACTIC telescope have been shown in Fig. 5.16 The 
shape selected images (after applying length, width, size, distance, and frac2) are used to 



plot the a-distribution, the example of which has been depicted in Fig. 5.18 The a is also 
known as image orientation parameter and its value is expected to be smaller for 7-ray show- 
ers having directions close to the optical axis of the telescope, while a nearly flat distribution 
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(a) a 7-ray like event (b) a proton event 

Fig. 5.16: Representative example of a 7-ray like and a proton shower image, recorded with the 
TACTIC telescope. Color code is an indication of the charge content of the pixel. 



is expected for an isotropic background of CR protons. For the TACTIC telescope, the excess 
of 7-ray events from a point 7-ray source is expected to be within an a range of 0° — 18°. 

5.5.6 Statistical significance of source detection 

In the ACT, the observations on a candidate 7-ray source are normally carried out in two 
modes; (i) ON-OFF mode: in this mode an on-source run is taken with the 7-ray source 
at the center of the camera, while off- source run is taken separately covering the same re- 
gion of the sky with no probable 7-ray source in the field of view, (ii) Tracking mode: 
the tracking mode is also known as the discovery mode of observations, wherein the can- 
didate 7-ray source is continuously tracked, without taking any off-source run. This mode 
of observation is preferred to maximize the on-source observation time and to increase the 
possibility of recording the flaring activity from the source. However, the systematic er- 
rors are relatively more as compared to the ON-OFF mode of data taking. Currently a third 
method called Wobble mode is also being used, wherein the source is kept at off-center po- 
sition and a few other off-center positions are marked to carry out OFF-source observations 
simultaneously. 

Certain percentage of CR events pass through all the 7 domain cuts and essentially be- 
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come background for 7-ray events, which need to be determined in order to calculate the 
statistical significance of the detected excess of 7-rays . While the off-source run provides 
an estimate of the background events in ON-OFF mode of observation, the tracking mode 
of data taking utilizes the on-source data itself for the calculation of background events. 
To estimate the expected background for the tracking mode, we use the events which have 
passed the dynamic supercut selection criteria except the a parameter. The distribution of 
a parameter is expected to be nearly flat between 0° and 90° for the background CR proton 
events, we therefore use the events with a value between 27° to 81°, for estimation of the 
background level. While the total range of a values i.e. 0° to 90° is divided into 10 bins 
with bin size of 9°, the on-source region is defined by the first two bins (0° to 18°). The third 
and last bins are not included for estimation of the background events, due to possible spill 
over of on-source events and comparatively less stability of events (possibly due to image 
truncation) from run to run, respectively. 

The significance of the signal can be estimated using either of the following two formulae 
proposed by Li and Ma [12171 (equation 5 and 17 of the paper) 



S 



vNoff 



or 



S = V2l Non In 



l + r] 
V 



Non + No 



^Non + r]moff 



+ Noff In 



ffJl 



(1 + ^) 



off 



Non + NoffJ_ 



(5.39) 
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(5.40) 

where Non is the number of events in the on-source region of a distribution (events in 0° to 
18°), Noff is the number of background events (events in 27° to 81°) and i] is the ratio of 
number of bins in on-source to number of bins in off-source region. In our case 77 = 3. 



Equation (5.39) is derived by estimating the standard deviation of the observed signal and it 
simply follows from the Poisson law of the counts i.e. Non and Noff. For the case ~ 1 



(0.5 < 77 < 1.5), the equation (5.39) can be used to evaluate the significance. The other 
equation ( 5.40| ), is derived by applying the method of statistical hypothesis test and can be 
applied in a general case of r/ 56 1 and also with fewer observed counts. 

While evaluating the significance of a 7-ray signal in our data analysis, we have used 



both the equations and found that equation (5.40) provides only marginal improvement in 



significance compared to evaluation by equation (5.39). 
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5.5.7 Energy spectrum estimation 

As already discussed, we use an ANN based technique to estimate the energy of the primary 
7-rays detected with the TACTIC telescope. The energy spectrum of the 7-ray source is then 
determined by binning these 7-ray photons in terms of energy. The ANN is trained using a 
Monte Carlo 7-ray sample and generates a weight file after achieving a satisfactory training. 
The ANN generated weight file is used in the analysis program, for determining the energy 
of the primary 7-ray photons, based on the size, distance and zenith angle parameters. 
The energy is determined only for those images which pass through all the 7-ray domain 
parameter cuts except a. Having an estimate of energy resolution of the telescope and also 
some idea about its threshold energy, the energy bins i.e. the energy bin width and first 
energy are selected, accordingly. The energy bin width is defined as 

AE^lnf (5.41) 

where is the middle energy of i*'* energy bin. The middle energy of a particular bin 
is related to lower and upper bounds of the same bin and to the middle energy of the other 
energy bins in the following manner, which follows from the above equation; 

E,^^E^^-e^ ■ Eu,^E^^-e^ (5.42) 

Erm+i — Efrii ' ^ ' Efni+i — -^mi ' i ' C . (5.43) 

After determining the energy, events are put in corresponding energy bins. The number of 
excess events for each energy bin are then determined by performing the frequency distribu- 
tion of their a parameter. 

The differential photon flux per energy bin has been computed using the formula 




(5.44) 



where SNi and d^{Ei)/dE are the number of events and the differential flux at energy Ei, 
measured in the i*^ energy bin of width SEi and over the zenith angle range of 0° — 45°, 
respectively. 7} is the observation time in the zenith angle bin with corresponding energy- 
dependent effective area (Aij) and 7-ray acceptance (rjij). The 5 zenith angle bins (j — 1 — 5) 
used are 0° - 10°, 10° - 20°, 20° - 30°, 30° -40° and 40° - 50° with simulation data available 
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Fig. 5.1 7: Light curve for Crab Nebula as recorded by the TACTIC imaging telescope during 2005- 
06 observation. Error bars shown are for statistical errors only. 



at 5°, 15°, 25°, 35° and 45°. The number of 7-ray events (5N.i) in a particular energy bin is 
calculated by subtracting the expected number of background events, from the 7-ray domain 
events. 



5.5.8 Light curve 

When the flux is determined with respect to time bins, the so-called light curve is generated. 
Conceptually, its determination is similar to that of the flux. Instead of determining the 
number of excess events in the bins of energy, it is determined in the bins of time and energy. 
If the statistics allow, this can be done fully differentially, otherwise integral flux can be 
calculated. Generally, in lACT astronomy, the time bin covers the whole observation night. 
In case of TACTIC, we calculate the number of 7-ray like events from the source direction 
for the whole observation night and generate a light curve in terms of photons / hour as a 



function of Modified Julian Day (MJD). Fig. 5.17 shows the light curve for the Crab Nebula 
as recorded by the TACTIC telescope during 2005-06 observations. The horizontal line in 
the figure shows the average 7-ray rate. 
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Fig. 5.1 8: Distribution of a parameter from the Crab Nebula direction (solid line and left hand scale) 
observed for ~101.44 /lonrs during 10 November 2005 - 30 January 2006. Horizontal 
line represents the background level per 9° bin derived using reasonably flat a region 
from 27° < a < 81°. Distribution of a parameter is also shown (dotted line and right- 
hand scale) from the Crab Nebula off-source direction observed for ~12.5 hours during 
2006. Error bars shown are for statistical errors only. 



5.6 Validation of analysis procedure 

The complete data analysis chain has been validated by analyzing the data, collected on the 
standard candle Crab Nebula, with the TACTIC imaging telescope, for ~ 101.44 hours dur- 
ing 10 November 2005 - 30 January 2006. The excess number of 7-ray like events selected 
after using the Dynamic Supercuts procedure was determined to be ~(839 ± 89) events with 



a statistical significance of ~ 9.64cr. Fig. 5.18 shows the frequency distribution of the a 



parameter. The corresponding average 7-ray rate turns out to be ~(8.27 ± 0.88) / hour. 
In order to generate the Crab Unit (CU) of flux in terms of the detected 7-ray like events, 
we analyzed the same data sample of Crab Nebula again by restricting the zenith angle of 
the observations from 15° to 45°. This zenith angle range is normally covered by Mrk421, 
during its observations with TACTIC telescope. The resulting 7-ray rate from Mrk421 (dis- 
cussed in next chapter) can then be converted into reference CU, while interpreting it's data. 
This re-analysis yielded an excess of ~(598 ± 69) 7-ray events in an observation time of 
~63.33 hours with a corresponding 7-ray rate of ~(9.44 ± 1.09) / hour, thus leading to 
the conversion: ICU = (9.44 ± 1.09) / hour. While one would have expected this rate 
to decrease because of increase in threshold energy, the reason behind this increase in 7- 
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Fig. 5.19: Spectrum of Crab Nebula as measured by the TACTIC telescope during 2005-06. Error 
bars shown are for statistical errors only. 



ray rate is the superior 7-ray acceptance of Dynamic Supercuts at higher zenith angles which 
over-compensates the decrease in the rate due to the increase of the threshold energy of the 
telescope. 

In order to test the validity of the energy estimation procedure, we have used the entire 
data i.e. for ~101. 44 /lowrs, collected on Crab Nebula. The 7-ray differential spectrum 
obtained after applying the Dynamic Supercuts and appropriate values of effective collection 
area and 7-ray acceptance efficiency (along with their energy and zenith angle dependence) is 



shown in Fig. 5.19 While determining the energy spectrum, we used the excess noise factor 
method [|97l for converting the image size in CDC counts to number of photoelectrons. The 
analysis of RGC data yields a value of Ipe = (6.5 ± 1.2) CDC when an average value of 
~ 1.7 is used for excess noise factor of the photomultiplier tubes. The differential energy 
spectrum of the Crab Nebula shown in Fig. 5.19 is a power law fit {d^/dE = f^E^^) 
with /o = (2.74 ± 0.19) x lO-^cm-^s-^Tel/"^ and V = 2.65 ± 0.06. The fit has a 
X^/dof = 0.53/6 (here and afterwards the denominator represents the value for the degrees 
of freedom ) with a corresponding probability of 0.997. The errors in the flux constant and 
the spectral index are standard errors. Excellent matching of this spectrum with that obtained 
by the Whipple and HEGRA groups [11851 , |50l reassures that the procedure followed by us 
for obtaining the energy spectrum of a 7-ray source is quite reliable. 



Chapter 6 



TeV observations of Markarian 421 



6.1 Introduction 

Mrk421 is a nearby blazar (z=0.030) and one of the brightest BL Lac objects located at a 
distance of ~ 134.1 Mpc (Ho = 71 km s'^ Mpc-\ Qm = 0.27, Qx = 0.73) [fT43ll . Blazars 
are the most extreme and powerful variable sources of photons with energies ranging from 
the radio to the 7-ray regimes. These objects are a subclass of AGN and are believed to have 
their jets more aligned toward the Earth as compared to any other class of radio loud AGN. 
They display high luminosity, irregular amplitude variability in all accessible spectral bands 
and have a core-dominated radio morphology with flat radio spectra, which join smoothly to 
the infra-red (IR), optical and ultra-violet spectra. In all these bands the flux exhibits high 
and variable polarization. With more observations, variety of features of such objects are be- 
coming evident and many different models can explain their observed properties [|106l . l265l . 
Their high variability and broadband emission make long-term observations of blazars very 
important for understanding their emission mechanisms and other related properties [I188L 

The blazar SED has a double-peaked structure in the vF,, versus frequency plot. Both 
peaks are found to vary, often both in intensity and peak frequency as the activity level of 
the blazar changes. The first peak is usually referred to as the synchrotron peak, both in 
leptonic and hadronic models for blazar emission. Further, it is generally believed to be the 
result of incoherent synchrotron emission from relativistic electrons and positrons, which are 
conjectured to be present in the magnetic fields of the jet. The origin of the second peak, 
usually referred to as the IC peak, is less well determined. 

In SSC models [|T05l . [T94l . 12271 . 12681 . [T63]| it is assumed that the synchrotron photons are 
up-scattered to higher energies by the electrons while in EC models HlOll , 11371 , 11611 , 12921 , 
I269l . ll36ll . these seed photons can come from the accretion disk, the broad-line region, the 
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torus, the local infrared background, the cosmic microwave background, the ambient photons 
from the central accretion flow or some combination of these sources. Hadronic models have 
also been invoked to explain the broadband spectra of blazars ll36]. 12241 . 1236II . It is proposed 
that the X-ray to 7-ray emission is synchrotron radiation from protons accelerated in highly 
magnetized compact regions of the jet. Other authors propose that proton-proton collisions, 
either within the jet itself or between the jet and ambient clouds, give rise to neutral pions 
which then decay to 7-rays [HI , [BO] , EH . 

Mrk421 was the first extragalactic source detected at TeV energies in 1992 by Whipple 
observatory 1125 5 II and subsequently by HEGRA group 11247 L using lACTs. The observations 
which led to the detection of Mrk421 at TeV energies were initiated in response to the 
detection of AGNs by EGRET. The initial detection indicated a 6a excess and the flux above 
500 GeV was approximately 30% of the Crab Nebula at these energies. The source has been 
regularly monitored by different groups since then Il305] , |25l] , HI] , [lOT] , [H , |69] , |270] , 
12581 . 13011 . It has been seen that the TeV 7-ray emission from Mrk421 is highly variable 
with variations of more than one order of magnitude and ocassional flaring doubling time 
of as short as 15 minutes II156I . l45l. The first clear detection of the flaring activity in the 
VHE emission of an AGN came in 1994 observations of Mrk421 by the Whipple telescope 
where a 10-fold increase in the flux, from an average level that year of approximately 15% 
of the Crab flux to approximately 150% of the Crab flux, was observed. The observations 
of Mrk421 in 1995 UllOU revealed several distinct episodes of flaring activity as in previous 
observations. Perhaps more importantly, they indicated that the VHE emission from Mrk421 
was best characterized by a succession of day scale or shorter flares with a baseline emission 
level below the sensitivity limit of the Whipple detector. 

The hypothesis that the VHE emission from Mrk421 could flare on sub-day time scales 
was borne out in a spectacular fashion in 1996, with the observation of two short flares II156L 
In the first flare which was observed on 07 May 1996, the flux increased monotonically dur- 
ing the course of ~ 2 hours of observations. The doubling time of the flare was ~ 1 hour. 
The next night the flux had dropped to a flux level of ~ 30% of the Crab Nebula flux, 
implying a decay time scale of < 1 day. The second flare, observed on 15 May 1996, 
although weaker, was remarkable for its very short duration: the entire flare lasted approx- 
imatey 30 minutes with a doubling and decay time of less than 15 minutes. Further in 
2001, Mrk421 underwent an extraordinary period of activity during which it was consis- 
tently brighter than the Crab Nebula over a three-month period II295II . 

Since its detection in the TeV energy range, Mrk421 has also been the target of sev- 
eral multiwavelength observation campaigns UE , [US , MB , llM , 1275] , 1276] , [IMl. One 



Chapter 6. TeV observations of Markarian 421 



141 



of the earliest multi-wavelength compaigns was organized in 1995 to measure the multi- 
wavelength properties of Mrk421. This compaign revealed, for the first time, correlations 
between VHE 7-rays and X-rays. The source also exhibited a large amplitude flare in VHE 
7-rays which was also clearly seen in the ASCA (X-ray) and Extreme Ultaviolet Explorer 
(EUVE) obsevations. The X-rays and 7-rays appear to vary together, limited by the one-day 
resolution of the VHE observations. Observations in April 1998 at TeV energies with the 
Whipple telescope and at X-ray wavelengths with the BeppoSAX satellite, established an 
hour scale correlations between the X-rays and 7-rays in a blazar [I295L 

Measurements of the VHE spectra are important for a variety of reasons. While the shape 
of the high energy spectrum is a key input parameter of AGN emission models, the way in 
which the spectrum varies with flux compared to lower energy observations provide further 
tests for the models. Several groups have determined the energy spectrum of Mrk421, both 
at low average flux levels of < 1 CU and from intense flares of > 2 CU. The recent results 
of these studies [1251 1 . 149] . l54l . 12101 suggest that the spectrum is compatible with a power law 
combined with an exponential cutoff of the form dN/ dE = K exp{—E/ Eq). It has also 
been reported that the spectrum hardens as the flux increases [|49l , l210ll . either because of an 
increase in the cutoff energy or a change in the spectral index itself. Differences in the energy 
spectrum of Mrk421 and Mrk501 have also been addressed to understand the 7-ray produc- 
tion mechanisms of these objects and absorption effects at the source or in the intergalactic 
medium due to interaction of VHE 7-rays with the EBL photons H 14411 . The energy spectra 
of Mrk501 and Mrk421 measured by HEGRA telescope system during their high states in 
1997 and 2001, respectively, are described by T = 1.92 and Eq = 6.2 TeV for Mrk501 and 
r = 2.19 and Eq = 3.6 TeV for Mrk421. The exponential cutoff in the spectrum of Mrk421 
starts earlier. Since both the sources are located at approximatey the same distance, the dif- 
ferences in the cutoff energies can be interpreted as an indication against the hypothesis that 
attributes the cutoffs to the pure intergalactic absorption effect. The earlier flux measurement 
of Mrk421 with the TACTIC telescope during 2004 observations have revealed a power law 
spectrum with a photon index of 2.4 ± 0.2 in the energy range of 2 — 9 Tel^ II258II . 

In this chapter, we present the TACTIC results obtained during 2005-2007 observation 
period. In addition, we also compare the TACTIC TeV light curves with those obtained 
with the RXTE/ASM ^ and Swift/BA^I [11 J in the X-ray energies for the contemporary 
periods. 
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Table 6.1 : Details of Mrk421 ON-source observations during the period 2005-06. 



Spell 


Observing period 


No. of nights 


Observation time(hours) 


I 


07 Dec 2005 - 1 1 Dec 2005 


5 


09.24 


II 


27 Dec 2005 - 09 Jan 2006 


12 


35.71 


III 


23 Jan 2006 - 07 Feb 2006 


16 


61.53 


IV 


19 Feb 2006 - 02 Mar 2006 


10 


34.54 


V 


19 Mar 2006-30 Mar 2006 


11 


31.14 


VI 


18 Apr 2006 -30 Apr 2006 


12 


29.55 


VII 


16 May 2006-21 May 2006 


05 


06.21 


II+III 


27 Dec 2005 - 07 Feb 2006 


28 


97.24 


Total 


07 Dec 2005 - 21 May 2006 


71 


207.92 



Table 6.2: Details of Mrk421 ON-source observations during the period 2006-07. 



Spell 


Observing period 


No. of nights 


Observation time (hours) 


I 


18 Dec 2006- 29 Dec 2006 


6 


8.96 


II 


17 Jan 2007 - 27 Jan 2007 


11 


26.54 


ni 


1 1 Feb 2007 - 22 Feb 2007 


11 


16.19 


IV 


09 Mar 2007 - 22 Mar 2007 


12 


26.76 


V 


06 Apr 2007 - 18 Apr 2007 


7 


4.61 


Total 


18 Dec 2006 - 18 Apr 2007 


57 


83.06 



6.2 TACTIC observations of Mrk421 (2005-07) 

We have used the imaging element of the TACTIC telescope at Mount Abu to observe 
Mrk421 during the period 07 December 2005 to 18 April 2007. These data were recorded 
using two trigger configurations NNP and 3NCT for the observation periods from 07 Decem- 
ber 2005 to 21 May 2006 and 18 December 2006 to 18 April 2007 respectively. The single 
pixel threshold was set to > 8/25 pe for 3NCT/NNP topological trigger logic. The telescope 
was operated at an individual SCR of 3 kHz, ensuring the safe operation of PMTs with typi- 
cal anode currents of < 3/i A. Mrk421 was observed for ~ 208 hours between 07 December 
2005 to 21 May 2006 and for ~ 123hours between 18 December 2006 to 18 April 2007. 
However, ~ 202 hours and ~ 83 hours of good quality data collected during 2005-06 and 
2006-07 respectively were used for further analysis. The observations were carried out in 
the tracking mode, where the source is tracked continuously without taking off- source data 
II256L This mode improves the chances of recording possible flaring activity from the candi- 
date source direction. Nearly 30 hours of the off-source data were recorded during 2005-06, 
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Table 6.3: Dynamic Supercuts selection criteria used for analyzing the TACTIC data. 



Parameter 


Cut Values 


length (L) 


0.11° <L< (0.235 + 0.0265 x In 3)° 


width (W) 


0.06° <W < (0.085 + 0.0120 x In 3)° 


distance (D) 


0.52° <D< 1.27°cos°-^^e ;(^=zenith ang.) 


size (S) 


3 > AbOd.c ;(6.5 digital counts=1.0 pe ) 


alpha (a) 


a < 18° 


frac2 (F2) 


F2 > 0.35 



in order to cross-check the background estimation. The total on-source data for 2005-07 
have been divided into different spells, where each spell corresponds to one lunation period. 
The zenith angle of the observations was <45°. Details of the these observations are given 
in Table [OandlOfor 2005-06 and 2006-07 respectively. 



6.3 Data analysis and results 

The general quality of the recorded data was checked by referring to the sky condition log 
and compatibility of the PGR and CCR with Poissonian statistics, as already explained in 
this thesis while discussing the data analysis procedure. The imaging data recorded by the 
telescope were corrected for inter-pixel gain variation and then subjected to the standard 
two-level "image cleaning" procedure II198II with picture and boundary thresholds of 6.5a 
and 3. Oct, respectively. The image cleaning threshold levels were first optimized on the Crab 
data and then applied to the Mrk421 data. The clean Cherenkov images were characterized 
by calculating their standard image parameters like length, width, distance, alpha, size 
and frac2 HI 821 , 12961 , 1198II . The standard dynamic supercuts II232II procedure was then 
used to segregate 7-ray like images from the huge background of CRs. The 7-ray selection 



criteria (Table 6.3 1 used in the analysis have been obtained on the basis of dedicated Monte 
Carlo simulations carried out for the TACTIC telescope. The cuts have also been validated 
further by applying them to the TACTIC Crab Nebula data analysis. 

A well established procedure to extract the 7-ray signal from the CR background using 
a stand alone imaging telescope is to plot the frequency distribution of the a parameter of 
shape and distance selected events. This distribution is expected to be flat for the isotropic 
background of CR events. For 7-rays , coming from a point source, the distribution is ex- 
pected to show a peak at smaller a values. Defining 0° < a < 18° as the 7-ray domain 
and 27° < a < 81° as the background region, the number of 7-ray events is then calculated 
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Table 6.4: Detailed spell wise analysis of 2005-06 data. 



Spell 


Obs. time 
(hours) 


7-ray 
events 


7-ray 
rate (hour^^) 


Signifi- 
cance((T) 


/dof 
(27° <a<8r) 


Prob. 


I 


9.24 


9 ±25 


0.97 ±2.67 


0.37 


2.97 /5 


0.705 


II 


35.71 


275±49 


7.70 ±1.37 


5.79 


5.77 /5 


0.329 


III 


61.53 


676 ±66 


10.99 ±1.07 


10.64 


8.57 /5 


0.127 


IV 


34.54 


91 ±47 


2.64±1.37 


1.94 


5.10/5 


0.404 


V 


31.14 


61±38 


1.96 ±1.23 


1.61 


2.44 /5 


0.785 


VI 


29.55 


123±33 


4.16 ±1.11 


3.86 


3.01 /5 


0.698 


All data 


201.72 


1236±110 


6.13 ±0.55 


11.49 


4.57 /5 


0.471 


II +III 


97.24 


951±82 


9.78 ±0.84 


12.00 


5.00 /5 


0.416 



by subtracting the expected number of background events (calculated on the basis of back- 
ground region) from the 7-ray domain events. Estimating the expected background level in 
the 7-domain by following this approach is well known HI 1411 and has been used quite exten- 
sively by other groups when equal amount of off-source data is not available. However, we 
have also validated this method for the TACTIC telescope by using separate off-source data 
on a regular basis and the a distribution of this data in range a < 81° is in good agreement 
with the expected flat distribution. The significance of the excess events has been finally 
calculated by using the maximum likelihood ratio method of Li & Ma [I2171I . 



6.3.1 Results of 2005-06 data analysis 

For the data recorded during 2005-06, we have first analyzed the data for each spell sepa- 
rately. Furthermore, it has also been ensured that these data are subjected to exactly similar 
analysis procedures as we did in the case of Crab Nebula to avoid any source dependent bias 
while determining the energy spectrum of Mrk421. The results of the spell wise analysis are 



presented in Table 6.4 An examination of this table clearly indicates that Mrk421 was most 
active during spell II and spell III observations. Fig. 6.1| ^a) gives the a-distribution when all 
the data collected for ~201.72 hours (07 December 2005 to 30 April 2006, spell I to spell 
VI) are analysed together. The total data yields an excess of ~(1236±1 10) 7-ray like events 
with a statistical significance of ~1 1.49cr. Fig. 6.1| ^b) gives the corresponding a-distribution 
for 29.65 hours of off-source data. The procedure for estimating the expected background 
in the 7-domain by using the background region is also consistent with the results from the 
off-source alpha plot (Fig. |6.1 ^b)), when its 7-ray domain events (0° < a < 18° events) 
are appropriately scaled up (to account for the difference in the on-source and off-source 



observation time) and compared with the 7-ray domain events of Fig. 6.1 a). The results 
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Fig . 6.1: (a) On-source a plot for Mrk42 1 when all the data collected between 07 December 2005 and 
30 April 2006 for ~ 201.72 hours are analysed, (b) Off-source a plot for ~ 29.65 hours 
of observation time. The horizontal lines in these figures indicate the expected background 
in the 7-domain obtained by using the background region (27° < a < 81°). 



of this calculation yield a value of ~(8320±238) events for the background level which is 
in close agreement with the value of ~(8171±52) events obtained when background region 



of Fig. 6.1 a) is used itself. It is worth mentioning that % /dof of the background region 



(indicated by column 6 of Table 6.4) is also consistent with the assumption that the back- 
ground region is flat and thus can be reliably used for estimating the background level in the 
7-domain. 

It is clear that the source was in a high state during 27 December 2005 to 09 January 2006 
and 23 January 2006 to 07 February 2006 at a combined average flux of ~(1.04±0.14)CU. 
It is worth mentioning here that the results of the Whipple group [[T3I. I188II also indicate that 
the source was in a high state during the period of their observations from 27 December 2005 
to 07 January 2006 and 23 January 2006 to 03 February 2006. Furthermore, our results for 
spell 5 (19 March 2006 to 30 March 2006) are also consistent with the Whipple observations 
from 23 March 2006 to 04 April 2006 when the source was seen to flare at a lower flux level. 
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Fig. 6.2: (a) Light curve for Mrk421 as recorded by the TACTIC imaging telescope between 07 
December 2005 and 30 April 2006 for a total observation time of ~ 201.72 hours, (b) Daily 
plot of the average count rate in the 2 — 10 keV energy range as determined by RXTE/ASM 
X-ray detector, (c) Daily plot of the average counts / cm? / s in the 15 — b()keV energy 
range as determined by SwjJf/BAT detector. Error bars shown are for statistical errors only. 



These observations clearly indicate that, despite differences in the observation times between 
the two telescope systems, Mrk421 was observed to be in a high state by both the systems 
for a prolonged duration, somewhat similar to the flaring episodes during January-February 
2001 and April-May 2004. 

6.3.2 Light curves 
TeV light curve 

We have examined day-to-day flux variability in the TeV light curve of Mrk421 for spell II 
and spell III data. The results indicate that we cannot claim any statistically significant flux 
variation on a nightly basis because of rather large error bars. The light curve of Mrk421 
as recorded by the TACTIC imaging telescope is shown in Fig. |6.2| ^a), where the average 
hourly 7-ray rates on daily basis have been plotted as a function of the modified Julian date. 
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The TeV light curve is characterised with a reduced value of 115.1/70 with respect to the 
mean level of 4.97 ± 0.49 7-ray photon events, with corresponding probability of 0.000554 
which is consistent with the variability hypothesis. 



ASM X-ray light curve 

The Rossi X-ray Timing Explorer (RXTE) satellite observes the time structure of astromoni- 
cal X-ray sources in the energy range 2— 200 keV. It is designed for studying known sources, 
detecting transient events. X-ray bursts, and periodic fluctuations in X-ray emissions. The 
RXTE has three instruments flU : 
(/) AU Sky Monitor (ASM): 2-10 keV 
(ii) Proportional Counter Array (PC A): 2 — 60 keV and 
(///) High-Energy X-ray Timing Experiment (HEXTE): 15 — 200 keV 
The ASM was built by the Centre for Space Research (CSR) at the Massachusetts Institute 
of Technology (MIT). It rotates in such a way as to scan most of the sky (~ 80%) every 
1.5 hours monitoring the long-term behavior of a number of the brightest X-ray sources, 
and giving observers an opportunity to spot any new phenomenon quickly. The ASM con- 
sists of three wide-angle shadow cameras equipped with position- sensitive Xenon propor- 
tional counters with a total collecting area of 90 cm^. The detector has a sensitivity of 
~ 20 mCrab d. 

The web page of RXTE ||51 makes available all ASM data products, including light 
curves data in ASCII format. Each raw data point represents the fitted source flux from 
one 90 second dwell and the Crab nebula flux is about 75 ASM counts / sec. Each "one- 
day average" data point represents the one-day average of the fitted source fluxes from a 
number (typically 5-10) of individual ASM dwells. The error value used in plotting the one- 
day averages is the quadrature average of the estimated errors on the individual dwells from 
that day Corresponding to the TACTIC TeV light curve discussed earlier, RXTE/ASM 



||5l light curve is shown in Fig. 6.2 ^b). This light curve is characterised with a reduced 
value of 574.4/70 with respect to the mean level of 2.12 ± 0.04 counts/sec. The probability 
obtained is very low which is consistent with the variability hypothesis. 



Swift/BAT light curve 

Swift is a multi-wavelength observatory which deploys three instruments [fT2]| 
(/) UV/Optical Telescope (UVOT): 170 - 600nm 
(//) X-ray Telescope (XRT): 0.3 - lOkeV and 
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Fig. 6.3: Differential energy spectrum of Mrk421 for the data collected between 27 December 2005 
- 07 February 2006 when the source was in a high state. 



(///) Burst Alert Telescope (BAT): 15 - 150 keV 

BAT is a large field of view (1.4 steradian) X-ray telescope with imaging capabilities. It has 
a coded aperture mask with 0.52 m? CdZnTe detectors. The BAT typically observes 50% to 
80% of the sky each day and has accumulated light curves of many non-GRB sources. In 
order to compare the TACTIC light curve of the source obtained during the observation 
period with that at 15 — 50 keV we have used Swift/B AT [11 J archived data of Mrk421 to 



plot its Ught curve which is shown in Fig. 6.2 c). It is characterised with a reduced value 
of 661/69 (very low probability, consistent with the variability hypothesis) with respect to 
the mean level of 0.0025 ib 0.000077 counts cm-^ s'K 



6.3.3 Energy spectrum of Mrk421 

We have estimated TeV energy spectrum of Mrk421 using the excess of 7-ray like events 
obtained during the high state of this source. The relevant details regarding the estimation of 
energy spectrum have already been discussed in the last chapter. The resulting differential 
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Table 6.5: Differential energy spectrum of Mrk421 derived using 2005-06 observations with the 
TACTIC telescope. 



Energy 


Energy bin width 


Differential flux 


Error in flux 


(TeV) 


(TeV) 


(photons cm~'^ s"^ TeV"^) 


(photons cm^^ s^^ TeV^^) 


1.0000 


0.403 


4.1399x10-11 


1.2975x10-11 


1.4918 


0.601 


1.2261x10-11 


1.7846x10-12 


2.2255 


0.896 


4.3461x10-12 


4.5502x10-13 


3.3201 


1.337 


1.1012x10-12 


1.8538x10-13 


4.9530 


1.994 


2.9385x10-13 


8.9087x10-14 


7.3891 


2.975 


8.1301x10-1^ 


4.2371x10-14 


11.0232 


4.439 


1.0425x10-14 


2.1035x10-14 



Table 6.6: Estimated intrinsic differential energy spectrum fitted parameters for three models. 



r 


Error in T 


fo 

7 cm-2 s-iTe\/-i 


Error in /o 

7 cm-2 s-17-ey-i 


Model 


2.99 


0.101 


6.69X 10-11 


6.38 X 10-12 


Primack [254J 


2.88 


0.108 


6.53 X 10-11 


6.19 X 10-12 


baseline[274| 


2.84 


0.093 


7.84 X 10-11 


6.20 X 10-12 


fast evolution [|274|| 



flux values for various energy bins are given in Table |6.5[ Since we know that the source 



was in a high state during the spell II and spell III of 2005-06 observations, we have used 



this data alone to determine the time averaged energy spectrum of Mrk421. Fig. 6.3 shows 
the differential energy spectrum after applying the Dynamic Supercuts to the combined data 
of spell II and spell III. A power law fit of the form 

d^/dE = (4.66 ± 0.46) x 10"^^ ^(-s.iiio.ii) ^^-2 ^-1 j-^y^i (g j) 

in the energy range 1 — 11 TeV with a x^/dof = 2.45/5 (probability = 0.78). 
A power law with an exponential cutoff 

d^/dE = foE-^exp{-E/Eo) (6.2) 

was also tried and the result yields the following parameters of the fit fo = (4.88 ± 0.38) x 
10-11 cm-2 s-i TeV-i, r = 2.51 ±0.26 and = (4.7±2.1)Te1/ with a xV^o/ = 0.88/4 
(probability = 0.93). The errors in the flux constant, the spectral index and cutoff energy 
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Fig. 6.4: TACTIC observed differential energy spectrum of Mrk421 derived using the data collected 
between 27 December 2005 - 07 February 2006 when the source was in a high state (bot- 
tom base filled triangle). De-absorbed source spectrum estmated by using the models (i) 
Primack given in II254II (hollow circle), (ii) baseline model given in II274II (top base filled 
triangle) and (iii) fast evolution model given in li2741 (filled circle). 



are again standard errors. While work on understanding the telescope systematics is still in 
progress, our preliminary estimates for the Crab Nebula spectrum indicate that the systematic 
errors in flux and the spectral index are < ± 40 % and < ± 0.42, respectively. 

The energy spectrum of Mrk421 as measured by the TACTIC imaging telescope in the 
energy region 1 — 11 TeV is compatible with both a pure power law fit and a power law with 
an exponential cutoff. However, a systematic deviation from a pure power law, although 



not statistically very significant, is evident in Fig. 6.3 at the energy values of 7.4 TeV and 
11.0 TeV. Difficulties like limited 7-ray event statistics coupled with rather large error bars 
do not allow us to claim the cutoff feature at a high confidence level. The cutoff energy 
Eq = (4:.7±2.1)TeV inferred from our observations is fairly consistent with the cutoff values 
of 3.6(+0.4 - 0.3),tai(+0.9 - 0.8) sysTeV and (4.3 ± 0.3)TeV reported by the HEGRA g5\ 
and the VERITAS [I210II groups, respectively. Furthermore, the results of the HESS group 
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Fig. 6.5: Distributions of image parameter a for the data, taken during 2006-07 observations. Hori- 
zontal line represents the mean background level per 9° bin derived by using resonably flat 
a region of 27° < a < 81°. Error bars shown are for statistical errors only. 



115411 . based on 9 nights in April and May 2004, also indicate that the time averaged energy 
spectrum of Mrk 42 1 is well described by a power law with an index T = 2 . 1 ± . 1 ^tat ± . S^j^s 
and an exponential cutoff at 3.1(+0.5 — OA) stat^O-^sysT eV . The HESS results also indicate 
that the cutoff signature in the energy spectrum is intrinsic to the source. 

We have also attempted to estimate the source intrinsic 7-ray spectrum by using the SED 
of the EBL provided by three models given in 025411 and 027411 (Primack, baseline and fast 



evolution models). The resulting intrinsic spectra are shown in Figure 6.4 for these three 
models. A power law fit to the estimated intrinsic data of the type d^/dE = foE~^ in the 



energy range 1 — 11 TeV yields values of F and /o listed in Table 6.6 As is clear from 



the Fig. 6.4, there is no significant difference between the measured source spectrum and 
the three de-absorbed spectra, which is not surprising as the source is associated with a low 
redshift value. 



6.3.4 Results of 2006-07 data analysis 

In this subsection, we present data analysis results of 2006-07 observations. When all the 



data are analyzed together, the corresponding results obtained are shown in Fig. 6.5 In this 
figure the histogram of the a parameter has been shown after having applied shape and orien- 



tation related imaging cuts given in Table 6.3 This distribution is almost flat and the number 
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Table 6.7: Monthly spell wise analysis of Mrk421 2006-07 data with only statistical errors. 



Spell 


7-ray 
photons detected 


Significance (cr) 


I 


1.67 ± 20 


0.08 


II 


49.33 ± 23 


2.11 


III 


8.33 ± 17.43 


0.48 


IV 


-56.00 ±21.00 


2.70 


V 


-13.00 ± 8.00 


-1.64 


I+II+III+IV+V 


-lO.OOi 42.00 


-0.25 



of 7-ray like events are —10 ± 42 with a statistical significane of -0.25o". Thereby indicating 
that the source was possibly in a low TeV emission state (below TACTIC sensitivity level) 
during the period of these observations. We have also carried out an exercise in which the 
2006-07 data are divided into five monthly spells I, II, III, IV and V of December 2006, 
January, February, March and April 2007 observations respectively and analyzed by using 
the same data analysis procedure. The number of 7-ray like events obtained are 1.67±0.2, 
49.33± 23.0, 8.33± 17.43, -56.0± 21.0 and -13.0± 8.0 for the respective spells, indicating 
that the source TeV 7-ray signal level has remained below TACTIC sensitivity during the 



monthly spells of 2006-07 observations. These results have been given in the Table 6.7 Fur- 
ther, these data have also been analyzed on nightly basis in order to explore the possibility 
of very strong episodic TeV emission. The results so obtained are depicted in Fig. 6.6^ , 
which shows the night to night variations of the 7-ray rate (7-rays /hour) for Mrk421 during 
2006-07 observations. This light curve is characterised with a reduced value of 37.59/28 
with respect to the mean level of -0.50±0.47 photon events, with corresponding probability 
of 0.10 which is consistent with the no variability hypothesis. The magnitude of an excess 
or deficit recorded on different nights is within ±2 a level for these observations and hence 
indicates the absence of a statistically significant episodic TeV 7-ray signal from the source 
direction. 



Next we compare the observed TeV light curve which is shown in the Fig. 6.61 with 
those of the source with the RXTE/ASM O (in the energy range 2-10 keV) and Swift/BKT 
ifTTll (in the energy range 15 — 50 keV) detectors. First we start with the RXTE/ASM con- 



temporary light curve shown in Fig. 6.6 5 which has been plotted by using the daily average 
count rates of the ASM from its archived data Q. The ASM light curve is characterised with 
a reduced value of 300.5/120 (probability obtained is very low which is consistent with 
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Fig. 6.6: Mrk421 light curves for 2006-07 observations. 



the variability hypothesis) with respect to the mean level of 0.58 ± 0.03 counts / sec. Next 
the Swift BAT contemporary light curve of the source obtained from its archived data [flTII 



is shown in Fig. 6.6 z. It is characterised with a reduced value of 158/116 (probability = 
0.0058 which is consistent with the variability hypothesis) with respect to the mean level of 
0.00042 ± 0.00000071 counts cm'^ s-\ 

We do not find any evidence for the presence of a statistically significant VHE 7-ray 
signal either in the overall data or when the data are analyzed on a month to month basis 
or a night to night basis during 2006-07 observations. We derive an upper limit of 1.44 x 
10'^"^ photons cm'"^ s~^ on the VHE 7-ray emission at 3a confidence level on the steady 
state flux from the source direction which is about 9% of the TACTIC detected integrated 
flux of the Crab Nebula above > ITeV by using the method of Helene [|180L This 
result suggests that the BL Lac object Mrk421 was possibly in a quiescent TeV state (below 
TACTIC sensitivity level) during 2006-07 observations. 



Chapter 7 



TeV observations of H1426+428 



7.1 Introduction 

H1426+428 is an established member of the class of extreme synchrotron BL Lac objects 
detected at TeV energies. Generally, objects of this class are characterized by a flat X-ray 
spectrum (F^ oc with F < 1) extending upto hard X-ray band. In the framework of SSC 
models, these objects are interesting for observations at TeV energies, as the IC peak of the 
SED is expected to be at higher energies. The work targeted to measure the IC peak position 
is hampered by the effect of pair production of VHE photons with low energy (optical and 
infrared) photons of the EBL. This source is located at a redshift of z = 0.129 for which the 
optical depth exceeds unity even for energies of a few 100 GeV. 

As already discussed in the last chapter, the SED of these objects extend from radio to 
7-ray energies with the double humped shape feature in a uFi, representation [[86l , 12881 , 
I289L The broadband SED of the BL Lac objects identified in X-ray surveys are significantly 
different compared to those identified in radio surveys and have been classified into HBL and 
LBL, depending on their X-ray to radio flux densities [I243II . The first hump in SED of these 
objects is expected at X-ray energies and generally accepted as originating from synchrotron 
radiation of relativistic electrons in the magnetic field around the object. The second peak is 
expected at 7-ray energies and in the framework of leptonic models [11051 , 11941 , 12271 , 12681 , 
IIOII . 11371 . 1161II is thought to be stemming from the IC scattering of low energy photons to 7- 
ray energies by the same population of relativistic electrons which produces the synchrotron 
radiation. Whereas, the hadron based models [|85]. ll30l .l40l. l235ll attribute the second hump 
and X-ray emission to some extent to processes involving protons which are also accelerated 
along with the leptons in the relativistic jets. These objects are extremely variable and are 
characterized by low emission quiescent states with occasional flaring when the flux can 
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increase by several orders of magnitude. For example, the HESS group detected the activity 
of the object PKS2155 -304 in July 2006 when its emission reached a flux level more than 
two orders of magnitude higher than its quiescent flux [|87l . 

H1426+428 was discovered in the 2-6 keV energy band by HEAO 1 II299II and was 
classified as a BL Lac object in 1989 [I261II . In an observation campaign performed by 
BeppoSAX in 1998-99, this source was identified as an extreme HBL because of its high 
synchrotron peak frequency, therefore a potential source of VHE 7-rays [I123II . 

The Whipple group reported the first VHE detection of this source in 2002, using the data 
recorded during 1999 -2001, at a statistical significance of about 6a [I187II . Subsequently it 
was confirmed by the CAT [I138II and HEGRA groups [j44] , |48l at statistical significance of 
5.2(7 and 7.5a respectively. The CELESTE and STACEE groups, both with non-imaging 
detectors have also observed the source and reported upper limits during the periods 2002- 
2004 [|270ll and 2001- 2002 [12011 respectively. However, the Crimean group has reported 
the detection of VHE 7-rays from this source using their GT-48 system in April 2004 with 
a statistical significance of 5.8cr [|151ll . More recently, the VERITAS group observed this 
source for 12.5 hours using their stereoscopic telescope system at Mt. Hopkins. Preliminary 
results obtained by them indicate a marginal excess of VHE photons at 3.2 a statistical level 
and they have placed an upper limit of about 3% of the Crab Nebula flux [I206II . In order 
to draw unbiased conclusions concerning the spectral shape of radiation relevant to the well 
defined states of the source, it is necessary to derive the energy spectrum in the desired energy 
range by using a single instrument during a relatively short observation period. The low 
fluxes and significant variability of such sources make this a rather difficult task. Hence, care 
should be taken while drawing conclusions about the energy spectrum based on information 
of fluxes obtained by different groups at different observation epochs. 

VHE 7-ray photons are believed to be absorbed by their interactions with the EBL lead- 
ing to an energy dependent horizon for viewing of VHE 7-ray sources. As a result of this 
absorption the energy spectrum of such a source may show a cut-off feature or steepening of 
the spectrum, but at the same time such features can also be attributed to the source itself. 
The interpretation of such features can be used as a probe for EBL measurements [I273II . In- 
terestingly, the recent detections of the AGN 3C 279 (z = 0.536) by the MAGIC group (EH, 
and H2356-309 (z = 0.165) and lES 1101-232 (z = 0.186) by the HESS group [[53 have 
infused a lot of interest in understanding of the 7-ray horizon in the VHE region. 

The limits on EBL density can be derived from two approaches: (?) assuming a limit to 
the hardness of the intrinsic spectrum of the AGN and deriving limits which assume very lit- 
tle about the EBL spectrum [|98l and (zz) assuming some shape of the EBL spectrum, based 
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Fig. 7.1 : Energy spectra of H1426-I-428. While the filled symbols indicate the measured fluxes re- 
ported by Whipple [248], HEGRA |48| and CAT 1 138] groups, the open symbols represent 
intrinsic fluxes after EBL correction based on three models. This figure has been adapted 
from BU. 



on theoretical modeling, and adjusting the normalization of the EBL density to match the 
VHE spectra [I134II . The latter can be more stringent but is necessarily more model depen- 
dent [12951 ■ Being at very large redshift, the TeV radiation from this extreme BL Lac object 
suffers severe intergalactic absorption and the measured spectrum can be used to derive lim- 
its on EBL. For any resonable EBL model, the absorption effect is at least a factor of 10 
or, more likely, factor of 100 [37J. This implies not only significant deformation of the 
primary spectrum, but also very large TeV luminosity of the source. It can be seen from 
Fig. |7.1[ the observed 7-ray spectra of the source, reported by Whipple [[248L HEGRA [i48l[ 



and CAT [1138B have a quite unusual shape i.e. the spectra are steep below 1 TeV but become 
flat above 1.5 TeV. However, the source intrinsic spectra after correction for the intergalactic 
absorption looks significantly different as compared to those observed. While for the EBL 
model suggested by Primack et al. [|253l . the intrinsic 7-ray spectrum has a power law form 
with a photon index ~ 1.5 as shown in Fig. |7.1| ^a), for other two models shown in Fig. |7.1| ^b) 
and (c) the reconstructed intrinsic source spectra have unusual forms with very sharp posi- 
tive slopes (E'^dN/dE oc E^) below 1 TeV for a model discussed in [441 or above 1 TeV 
suggested in [|22T]| . 

The deformation of the primary TeV 7-ray spectrum due to the intergalactic absorp- 
tion is characterised by the factor e^""*^^'^), where t{z,E) is the energy-dependent optical 
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depth for a 7-ray photon of energy E emitted by a source at a redshift z. The spectral 
deformation factors calculated for the above mentioned 3 different EBL models show that 
above 300 GeV the intergalactic absorption leads to a strong steepening of the spectrum of 
7-rays from H1426+428, but from ITeV to several TeV the spectrum is deformed only 
slightly, although the suppression of the absolute flux may be as large as a factor of 100. 
Therefore the intrinsic source spectrum of the 7-rays i.e. spectrum after correction for the 
intergalactic absorption 

J,,(E) = J,,,(E)e^(^'^) (7.1) 

is expected to be quite different from the observed spectrum Jobs{E). 

The sharp positive slopes of the 7-ray spectra (E'^dN/dE oc E^) in Fig. |7.l| ^b) and (c) 
can not be easily explained by conventional models invoking either a leptonic or a hadronic 
origin of 7-ray emission of BL Lac objects. Even a decent power law intrinsic TeV spec- 
trum of H1426-I-428 {E'^dN/dE oc E^-^), obtained for a model of EBL spectrum shown in 
Fig.jTTTJ^a), does not yet guarantee an easy explanation of the broad-band SED of this source. 
It may require a significant revision of the current conceptual view according to which the 
synchrotron (X-ray) peak in the SED dominates over the IC (TeV) peak. In the case of 
H1426-I-428, the detected TeV flux after correcting for intergalactic absorption becomes at 
least an order of magnitude larger than the X-ray energy flux [|124[| . Thus, the TeV lumonos- 
ity significantly exceeds the X-ray luminosity and it is clear that this should have an impact 
on the TeV blazar models [i37il . Generally, the two bumps in the SED of TeV blazars are 
always called "synchrotron" and "IC" peaks which tacitly implies that the leptonic origin of 
the TeV radiation is accepted as an undisputed fact, and all possible difficulties with the in- 
terpretation of the data within this concept can be blamed on the imperfectness of the existing 
models. However, the TeV observations of sources like HI 426-1-428, indicate that it is quite 
important to develop in parallel with improvement of leptonic models, alternative scenarios, 
attributing, for example, the TeV emission to the synchrotron radiation of protons, and the 
X-ray emission to the synchrotron radiation of primary (directly accelerated) or secondary 
(pair produced) electrons fl37l . 

The apparent luminosity of 7-rays after correcting for the intergalactic absorption appears 
at least an order of magnitude larger than the synchrotron luminosity. Within the one-zone 
SSC model, this would imply a very low magnetic field and a very large density of relativistic 
electrons. The EC models and the models of electromagnetic cascade in radiation-dominated 
environments allow more flexibility in modeling the unusual SED of this source. Alterna- 
tively, the hard intrinsic TeV 7-ray spectrum of H1426-I-428 and high L^/Lx ratio (> 10) 
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can be explained by the proton-synchrotron model (assuming ~ fraction of the proton 
synchrotron radiation suffers internal absorption). This scenario couples X-rays to TeV 7- 
rays through synchrotron radiation of secondary electrons, therefore we must expect very 
tight X-ray s/Te\^ 7-ray correlations. The secondary electrons are promptly cooled in the 
strong magnetic field, so one may assume that the X-ray will arrive almost simultaeously 
with the bulk of the TeV 7-rays . However, more complicated correlations can not be ex- 
cluded, especially if the magnetic field drops significantly at the outskirts of the emitting 
region, so an essential part of the secondary electrons are produced and cooled in the low 
field environment. This unfortunately makes the predictions about the X-ray/TeV 7-ray cor- 
relations less robust. 

In this chapter, we present TACTIC results obtained from our observation compaign on 
this source during the period 2004-07. In addition, we also compare the TACTIC TeV light 
curves with those obtained with the RXTE/ASM ^ at X-ray energies for the contemporary 
periods. 



7.2 TACTIC observations of H1426+428 (2004-07) 

We have observed the BL Lac object HI 426-1-428 using the imaging element of the TAC- 
TIC telescope, during the period 2004-07. The imaging camera was configured with inner 
225 pixels for recording ACR images, while innermost 121 pixels were used for the trigger 
generation. The trigger is based on the 3NCT configuration for 2004 and 2007 data while 
2006 data were collected with NNP. The single pixel threshold was set to > 8/25 pe for 
3NCT/NNP topological trigger logic. 

In order to maximize the on-source observation time and to increase the possibility of 
recording flaring activity from the source, we have made these observations in the continuous 
source tracking mode. In this mode, the source was always placed at the centre of the camera 
while observing it. The zenith angle range covered during these on-source observations was 



from 18° to 45°. Details of these observations are given in Table 7.1 



7.3 Data analysis and results 

While analysing, the data were subdivided into 3 spells each corresponding to the observa- 
tions made in a particular year. Standard data quality checks, e.g. compatibility of the PCR 
and CCR with Poissonian statistics, have been performed in order to evaluate the system be- 
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Table 7.1 : Details of TACTIC observations on H1426+428. 



Spell 




Months 


Ohsprvatinn 
Dates 


Total data 

(JlOUT S) 


Data selected 

(JlOUT S) 


I 


2004 


March 
April 


22-26, 28,29,31 
12-15, 18,20-22, 24-26 


71.0 


44.8 


II 


2006 


March 


26-31 










April 


1,2, 4,21,23-30 


34.9 


18.4 






May 


01,02, 19,21,23-24, 29-30 






III 


2007 


March 


16-27 










April 


12-15, 18-20, 22-26 


137.7 


102.5 






May 


06-15, 17-19, 21-22 










June 


07-11 







havior and the general quality of the recorded data. As a result of this exercise, 77.9 hours of 
data were rejected and the rest of the data were used for detailed analysis. While performing 
detailed analysis of the data, firstly all the data were subjected to pedestal removal by us- 
ing the standard two-level "image cleaning" procedure [I198II . The image cleaning threshold 
(picture and boundary) levels were first optimized on the Crab Nebula data HI 6611 and then 
applied to the data presented here. The digital counts of each pixel in the clean images were 
then corrected for inter-pixel gain variation. The clean Cherenkov images were character- 
ized by calculating their standard image parameters like length, width, distance, alpha, 
size and frac2 [|182M1981 , |296L The standard extended supercuts [I232II procedure was then 
used to separate 7-ray like images from those due to the huge background of CRs. The 7-ray 



selection criteria based on the imaging cuts, given in Table 7.2 used in the analysis have been 
obtained on the basis of Monte Carlo simulations carried out for the TACTIC telescope. It 



may be noted that slightly different 7-ray domain cuts (Table 7.2 ) have been used for the data 



collected on H1426-I-428 as compared to those used for Mrk421 (Table 6.3 1. This is done be- 
cause HI 426-1-428 is normaly observed between March to June (see Table 1 .1) during which 
period the telescope operates at slightly higher threshold as compared to its winter (October 
to Februbary) operation due to variations in sky transparency between the two seasons. 

The next step of the data analysis deals with estimating the number of 7- ray events 
from the source direction in the presence of an overwhelming background of CR events. 
This is done by plotting the frequency distribution of the a parameter of length, width, 
size and distance selected events. The 7-rays from a point source have smaller values of 
the a parameter as mentioned earlier in this thesis, whereas the CR events, because of their 
isotropic nature, are uniformly distributed in all a values ranging from 0° to 90°. For the 
TACTIC system we find that the a range for 7-ray events is found to be from 0° to < 18°. 
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Table 7.2: Extended Supercuts selection criteria used for analyzing the TACTIC data. 



Parameter 


Cut Values 


length (L) 


0.11° <L< (0.235 + 0.0265 x In S)° 


width (W) 


0.065° <W < (0.085 + 0.012 x In 5)° 


distance (D) 


0.5° <D< 1.27° 


size (S) 


S > 350d.c; (6.5 digital counts = 1.0 pe ) 


alpha (a) 


a < 18° 


frac2 (F2) 


F2 > 0.38 



The flat distribution of a in the range 27° < a < 81° has been used to determine the 
background events present in the signal region BOlll . The number of 7-ray events is then 
calculated by subtracting the expected number of background events (calculated on the basis 
of the background region) from the 7-ray domain events and the significance of the excess 
events has been finally calculated by using the maximum likelihood ratio method of Li & 



7.3.1 TeV results 

When all the data recorded during the year 2004 are analysed together, the corresponding re- 



sults obtained are shown in Fig. 7.2 1, wherein the histogram of the a parameter has been plot- 



ted after having applied length, width, size and distance imaging cuts given in Table |7.2 
As is clear from this figure, the distribution is almost flat and the number of 7-ray like events 
within the 7-ray domain of the distribution are 19±48, thereby indicating that the source was 
possibly in a low TeV emission state (below TACTIC sensitivity level) during the period of 
these observations. Further the same data were also analysed on a nightly basis to explore the 
possibility of very strong episodic TeV emissions. The corresponding results obtained are 



depicted in Fig. 7.3 1, which shows the day-to-day variations of the 7-ray rate {'j -rays/ hour) 
for 2004 observations. This light curve is characterised with a reduced value of 19.48/12 
with respect to the zero degree polynomial fitted constant value of 0.28±0.97 photon events, 
with corresponding probability of 0.1 1 which is consistent with the no- variability hypothe- 
sis. The magnitude of an excess or deficit recorded on different nights is within ±2 a level 
for 2004 observations and hence indicates the absence of a statistically significant episodic 
TeV gamma-ray signal from the HI 426-1-428 direction. 

We have followed the same data analysis methodology while analysing the data recorded 



during the years 2006 and 2007 and the corresponding a histograms are shown in Fig. 7.2 5 



and 7.2; respectively. As is clear from these figures the distributions are almost flat and 
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Fig. 7.2: (a) Distribution of image parameter a after applying the cuts given in Table 7.2 for 
ii.S hours, ISAhours and 102.5 hours of data taken on H1426+428 during (a) 2004, 
(b) 2006 and (c) 2007 observations respectively. The horizontal lines in the figure indi- 
cate the expected background in the 7-domain obtained by using the background region 
(27° < a < 81°). The error bars shown are for statistical errors only. 



the number of 7-ray like events within the 7-ray domain of the distributions are 68±35 and 
-54±80 for the 2006 and 2007 observations respectively, thereby indicating that the source 
was possibly in a low TeV emission state (below TACTIC sensitivity level) during the pe- 
riod of these observations as well. Further these data were also analysed on a nightly basis 
and the corresponding results obtained are depicted in Fig. |7.3| j and |7.3| : respectively, again 
showing the day-to-day variations of the hourly 7-ray rate with time. Light curves are char- 
acterised with the reduced values of 12.66 /13 and 31.31 /34 respectively, with respect 
to the zero degree polynomial fitted constant value of 2.22 ±1.79 and -0.49±0.74 photon 
events with corresponding probabilities of 0.47 and 0.69 respectively which are consistent 
with the no-variability hypothesis. The magnitude of an excess or deficit recorded on differ- 
ent nights is also within ± 2 a level for 2006 and 2007 observations thereby again indicating 
the absence of a statistically significant episodic TeV 7-ray signal during these observations 
as well. 
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Fig. 7.3: TACTIC light curves for (a) 2004, (b) 2006 and (c) 2007 observations and corresponding 
ASM light curves are shown in (d), (e) and (f) respectively. 



In Fig. 7.4[ we have shown the a histogram obtained when the data for all the years 
were combined together which is also consistent with the flat distribution. There is thus no 
evidence for the presence of a statistically significant 7-ray signal from the source direction 
during the period of our observations. We have placed an upper limit of < 1.18 x 10^^^ 
photons cm'"^ s^^ at 3a level (13.1% of the TACTIC detected Crab Nebula integrated 7-ray 
flux) on the integrated TeV 7-ray flux above 1 TeV from the source direction, using the 
method of Helene II180B . In Table [73| we have compiled all results obtained on a yearly 
spell basis as well as for the case when all the data are considered together. 



7.3.2 ASM X-ray light curves 

In order to examine any correlation between VHE 7-ray and X-ray (2 — 10 keV) source 
emissions, the Whipple group who discovered this source has done correlation studies for 
this source and no evidence for significant correlation was found, rather their 2001 data 
show some evidence for an anticorrelation between the average monthly rates L187J. In the 
present work we also compare the TACTIC light curves with those of RXTE/ASM X-ray 
observations in the energy range of 2 — 10 keV. For this we have used the daily average count 
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Fig. 7.4: Distribution of image parameter a after applying cuts given in Table 

165.7 hours of data on H1426+428. The horizontal line in the figure indicates the expected 
background in the 7-domain obtained by using the background region (27° < a < 81°). 
The error bars shown are for statistical errors only. 



Table 7.3: Yearly spell wise and consolidated data analysis results obtained on H1426+428 with the 
TACTIC telescope. 



Year 


Obs. time 


Excess/Deficit of events 


Significance 


Upper Limit at 3cr level 




(h) 


from source direction 


(^) 


photons cm^^ s^^ 


2004 


44.8 


19±48 


0.39 


< 2.13 X 10-12 


2006 


18.4 


68±35 


1.95 


< 3.72 X 10-12 


2007 


102.5 


-54±80 


-0.68 


< 1.54 X 10-12 


Total 


165.7 


33±100 


0.33 


< 1.18 X 10-12 



rate of ASM from its archived dataflSl to obtain the light curves for the contemporary periods. 
These light curves are shown in Fig. |7.3| i, 7.3^ and |7.3| f for the corresponding TACTIC 
observation periods of 2004, 2006 and 2007 respectively. ASM light curves are characterised 
with the reduced values of 5.27/12, 1 1 .03/13 and 43.78/34 with respect to the zero degree 
polynomial fitted constant values of -0.2±0.088, 0.095±0.082, and 0.152±0.046 counts and 
corresponding probabilities of 0.96, 0.60 and 0.17 respectively, thereby indicating that these 
light curves are consistent with the constant flux hypothesis and no variability is observed in 
a time scale of a day or more in the RXTE/ASM energy band also. We have not done any 
correlation studies between the TACTIC and RXTE/ASM observations, mainly because of 
non-detection of the source in the present work. It may be noted here that we have not used 
the EBL absorption effect while deriving the upper limit. We have used ASM data [|5|| in 
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the presented work just to inspect the contemporary RXTE source light curves in the energy 
band of 2 — 10 keV and to compare with those obtained with the TACTIC observations. 
In the earlier reported VHE detections, no evidence for the significant correlation between 
X-ray and VHE 7-ray was found HI 8711 as has been mentioned earlier. 



Chapter 8 



Discussion and conclusions 



The aim of this work was to study the VHE 7-ray emission from two extragalactic sources 
(of type blazar) Mrk421 and H1426+428. While Mrk421 is the closest {z = 0.030) known 
VHE 7-ray emitter, H1426+428 was the farthest (z = 0.129) known source at the time of 
start of this work. However, the most distant source as of now is the FSRQ 3C279 (z = 
0.5362) detected by the MAGIC telescope flM . The detection of this source, signifies the 
importance of having a detector with lower energy threshold as the universe is comparatively 
more transparent to low energy photons. 

Blazars are the most powerful variable sources of electromagnetic emission ranging from 
the radio to the 7-ray regimes. Their SED show a two peak structure and both of them are 
found to vary in intensity as well as in peak frequency with the emission level of the blazars. 
The first peak is commonly believed to originate from synchrotron emission both in leptonic 
and hadronic models for blazar emission. The origin of the second peak (usually called 
IC peak) is less well determined and is still being debated. Presently, the leptonic models 
represent the most preferred concept of 7-ray emission from blazars and are divided into 
the following two sub-classes, (i) In SSC models M5\ , , 1227] , |268] , MM the seed 
photons are synchrotron photons emitted by electrons which are then upscattered with the 
same electron population via IC process, (ii) In EC models [fTOT] . [TT71 . [T6T] . 12921 . 1269) . [T36l . 
the seed photons are ambient infrared or optical photons, photons of CMBR or thermal 
radiations either directly from the accretion disk or first scattered by the surrounding gas and 
dust clouds. On the other hand, hadronic models assume that the observed 7-ray emission 
is from the accelerated protons interacting with ambient matter [I252L photon fields [12241 . 
|225],|226l, magnetic fields [136J or both magnetic and photon fields [I235] . l236ll . 

In order to understand the properties of blazars and the underlying emission mechanisms 
responsible for observed variability over the entire electromagnetic spectrum, it is necessary 
to detect more and more sources. Till date 42 extragalactic sources have been detected (T^. 
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Most of these detected extragalactic VHE sources belong to the blazars category of AGN 
except two radio galaxies, two starbust galaxies and three unidentified sources. Further, 
majority of the detected blazars belongs to HBL type, while LBL and FSRQ count only 
three sources each. Four recently detected sources belong to type IBL. 

8.1 Mrk421 

Using TACTIC 7-ray telescope, we have observed Mrk421 during two spells: (?) 07 Decem- 
ber 2005 to 30 April 2006 and (//) 18 December 2006 to 18 April 2007, with a total on-source 
observation time of ~202 hours and ~123 hours respectively. The source was detected in a 
bright state with a statistical significance of ~12.0(j during the period 27 December 2005 to 
09 January 2006 and 23 January 2006 to 07 February 2006 in about ~97 hours of observa- 
tion time. Total number of TeV photons detected were 951 ± 82 during these observations. 
The combined average flux was measured to be ~(1.04±0.14)CU. The Whipple group has 
also detected this object in a high state during following periods: (/) 27 December 2005 - 07 
January 2006 and (//) 23 January 2006 - 03 February 2006 [I188L Furthermore, our results 
for 19 March 2006 to 30 March 2006 are also consistent with the Whipple observations from 
23 March 2006 to 04 April 2006 when the source was observed at a lower flux level. These 
observations clearly indicate that, despite difference in the observation time between the two 
telescope systems, Mrk421 was observed to be in a high state by both the systems for a pro- 
longed duration, somewhat similar to the flaring episodes during January - February, 2001 
and April - May, 2004 |f T04l . 

We have also determined the observed time-averaged differential 7-ray spectrum of the 
source in the energy range 1 — 11 TeV. This spectrum when fitted with a pure power law 
of the form d^{E)/dE = K R-^ yields V = 3.11 ± 0.11 and K = (4.66 ± 0.46) x 
10^^^cm^'^s''^TeV^^ with (x^/dof = 2.45/5 ). In addition, we have also tried another 
power law function with an exponential cutoff of the form d(p{E) / dE = K E~^exp{—E / Eq) 
which yields K = (4.88 ± 0.38) x 10-^^cm-^s-^TeV-\ T = 2.51 ± 0.26 and Eq = 
(4.7 ± 2.1)TeV with a (x^/dof = 0.88/4). We find that the later function fits the observed 
spectrum well with lower value of reduced x^- The cutoff energy Eq = (4.7 ± 2.1)Te\^ 
inferred from our observations is fairly consistent with the cutoff values of 3.6(+0.4 — 
0.3),t„t(+0.9 - 0.8)sysTeV and (4.3 ±0.3)TeV reported by the HEGRA [|45J and the VERI- 
TAS [I210II groups, respectively. Furthermore, the results of the HESS group [[54|| , based on 9 
night observations in April and May 2004, also indicate that the time averaged energy spec- 
trum of Mrk421 is well described by a power law with index F = 2.1 ± O.lgtat ± O-'^sys with 
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Fig. 8.1 : Comparison of the observed spectral energy distribution of Mrk 421 in the VHE range 
obtained by the HESS, MAGIC, VERITAS, Whipple, HEGRA, CAT and TACTIC |54 , 
|69] , |305] , [32| , [39] , |45] , |25l| , [Ml- The TACTIC observed differential energy spectrum 
during 2005-06 observations is in the energy range 1 — 11 TeV. The arrow represents an 
upper limit obtained at a 3(j confidence level on the integrated 7-ray flux above ITeV from 
this source direction during 2006-07 observations. 



an exponential cutoff at 3.1 (+0.5 — OA)stat ± 0.9sysTeV. The HESS results also indicate 
that the cutoff signature in the energy spectrum is intrinsic to the source. However, during 
the observation period from 18 December 2006 to 18 April 2007, the source was found to be 
in a low emission state and we have placed an upper limit of 1.44 x 10^^^ photons cm~^ s^^ 
at 3cr level on the integrated TeV 7-ray flux above 1 TeV from the source direction. 

In addition, we have compared the TACTIC TeV light curves with those obtained with 
the RXTE/ASM ^ and Swift/BAI [|TT| satellite-based experiments which are sensitive 
in the energy ranges 2 — lOkeV and 15 — 150 keV respectively. Due to the large er- 
ror bars associated with TACTIC light curves coupled with non-simultaneous observations 
with RXTE/ASM and SwiftfBAI, we could not obtain any statistically significant correlation 
among these data. 
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From the TACTIC observations of Mrk421 presented in this thesis, we conclude that 
during 27 December 2005 to 09 January 2006 and 23 January 2006 to 07 February 2006 the 
source was in a high emission state characterized by a spectral index of F = 2.51 ± 0.26 and 
cutoff energy of Eq = (4.7 ± 2.1)TeV. During rest of the observation period the source was 
possibly in a TeV quiescent state. 

Further, we have obtained the SED of Mrk421 in the VHE range using TACTIC tele- 
scope ll30T1l and compared it with those obtained by the HESS, MAGIC, VERITAS, Whip- 
ple, HEGRA and CAT telescopes [[51 , M , [305] , |32] , |39] , |45] , |25ll. This comparision is 
shown in Fig. |8.1[ The source 2006 observations (28 April to 14 June) reported in [1761 have 
confirmed the indications of rising peak energy with increasing flux, as expected in leptonic 
acceleration models. It may however be noted that the TACTIC observations for the year 
2006-07 reported in this thesis do not overlap with those reported in the f76J. Therefore 
there is no conflict between these two observations. 

If we compare the earher high states of this source the flare of May 1996 detected by 
the Whipple group remains the most intense one so far at the level of about 10 times of 
the Crab Nebula flux. Therefore, in order to increase the probability of recording a similar 
intense flare more consistent multiwavelength observations are needed with the present age 
detectors. It seems that such flares are very rare in nature but have a great potential to provide 
a deeper insight into the astrophysics of AGN. 



8.2 H1426+428 

The BL Lac object HI 426-1-428 is an established source of TeV 7-rays and detections of 
these photons from this object also have important implications for estimating the EBL, in 
addition to the understanding of the particle acceleration and 7-ray production mechanisms 
in the AGN jets. We have observed this source for about 2AAhours in 2004, 2006 and 
2007 with the TACTIC 7-ray telescope. Detailed analysis of these data do not indicate the 
presence of any statistically significant TeV 7-ray signal from this source in the total data set 
or when the data set is divided into monthly subsets. Accordingly, we have placed an upper 
limit of < 1.18 X 10^^'^ photons cm^'^ s^^ on the integrated 7-ray flux at 3a significance 
level [I302]| . 

In Fig. |8.2[ we show the available light curves of the source from 1999 to 2004 ob- 
served by the Whipple, CAT, CELESTE and GT-48 ^W] , [I38] , |270] , \T5^ groups. As is 
clear from this figure, the GT-48 telescope has recorded the source at more than 200% of 
the GT-48 detected Crab level on April 19, 2004 (MJD 53115) at 3.9 a level ifBTI.fBOl. 
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Fig. 8.2: H1426+428 TeV light curves obtained by various groups, namely Whipple II187II . 
CATini, CELESTEUva and GT-48II5I1. 



It may be noted here that the TACTIC 2004 results discussed above on this source are in 



agreement with those obtained with the CELESTE system 1270] around 100 GeV during the 
same period. However, our results are in conflict with those of GT-48 telescope above 1 TeV 
during the period from 15 April to 25 April, 2004, wherein a 7-ray signal at 5.8cr level was 
reported ifTSTIl . 

In addition, we have compared RXTE/ASM [|5l X-ray light curves (2 — 10 keV) with 
those obtained with the TACTIC observations. We do not find any enhanced activity of 
the source at X-ray energies too. Further, we have also compared the available TeV light 
curves of the source from 1999 to 2004 observed by the Whipple, CAT, CELESTE and GT- 
48 groups [|I87| , UMI , I22Q] , USB. It is interesting to point out here that the TACTIC 2004 
results are in close agreement with those obtained by the CELESTE group above 100 GeV 
during the same period. We conclude that the emission level of H1426-I-428 was below the 



TACTIC sensitivity with the estimated flux level ~ 13.1% of the Crab Nebula. In Fig. 8.3 



we have compared the TACTIC upper limit obtained during 2004-2007 with the source VHE 
spectra results obtained during different epochs by the Whipple [I187L HEGRA [48 J and 
CAT alongwith the upper hmits provided by the CELESTE [12701, STACEE [20| and 
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Fig. 8.3: TACTIC upper limit on H1426+428 derived using 165.7 hours of data has been shown 
along with results obtained by Whipple EUl, HEGRA El, CAT CELESTE |[270l . 
STACEE O and VERITAS 12061 groups. 



VERITAS II206II collaborations. Further, those groups who have reported VHE 7-ray de- 
tections have derived source spectra which are consistent with the expected signature of 
absorption of TeV 7-rays by EBL [|248] , |44| , |48] , [HI . 

We feel that the long term multi- wavelength observations of the source are required to 
possibly record more flaring episodes of the source by using the technologically advanced 
telescopes. These observations in the VHE regime would also possibly record a better statis- 
tics of VHE photons to derive the source related observed and intrinsic spectra corresponding 
to various high states. Indeed this would also help in understanding the particle acceleration 
and 7-ray production mechanisms in the AGN jets, in addition to the VHE photon horizon 
and EBL related quest. 



8.3 Future outlook 



Exploring the 7-ray sky in the VHE energy range is expected to lead to a potentially rich har- 
vest of astrophysical discoveries. This has been amply demonstrated in the last few years by 
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ground based lACTs such as HESS, MAGIC and VERITAS at 7-ray energies > 100 GeV. 
Further, at still lower energies, the new generation of satellite based 7-ray telescopes partic- 
ularly Fermi which was launched in 2008 has heralded a new era in 7-ray astronomy [fTTl . 
The Fermi LAT which is sensitive to 7-rays with energies in the range from 20 MeV to 
100 GeV has to date detected about 1500 objects. However, due to low flux, instruments for 
detection of VHE 7-rays of more than few tens of GeV, require a large effective detection 
area which is not possible with space based instruments. On the other hand, ground based 
instruments allow much larger detection area. The recent rapid growth in the field of VHE 
7-ray astronomy was therefore primarily driven by lACTs which reach a sensitivity of ~ 1% 
of the Crab Nebula flux for typical observation time of 25 hours. The sensitivity degrades to- 
wards lower energies due to threshold effects and towards higher energies due to the limited 
detection area [I126L 

For many unanswered questions in the field of 7-ray astronomy, low energy threshold 
instruments are expected to be promising. For example, VHE 7-ray photons from distant 
sources interact in the intergalactic space with infrared and other low energy photon fields 
leading to an energy dependent 7-ray horizon for vewing VHE 7-ray sources. The obser- 
vations at lower energy threshold give access to 7-ray sources at high redshifts. The low 
threshold energy can be attained by increasing the light collector area of the telescopes and 
installing them at higher altitudes where the photon density of the atmospheric Cherenkov 
events is higher ll42l . 

As part of augmentation of our activities in 7-ray astronomy we are setting up a large 
area imaging 7-ray telescope MACE (Major Atmospheric Cherenkov Experiment) at Hanle 
(32.8° N, 78.9° E, 4200 m asl) in the Ladakh region of India. This high altitude site offers an 
average of about 260 uniformly distributed spectroscopic nights per year which is a major 
advantage in terms of sky coverage for source observations. Located closer to the shower 
maximum the Cherenkov photon density at Hanle is substantially high as compared to that 
at sea level [I125II . The higher photon density along with the low background light level at 
this site helps in lowering the energy threshold of the Cherenkov telescope being setup there. 

The MACE telescope with high resolution imaging camera is designed to operate in the 
suh-TeV energy range. The altitude- azimuth mounted telescope will deploy a 21m diameter 
parabolic light collector made of 356 panels of 984 mm x 984 mm size with each panel 
consisting of 4 spherical mirror facets of 488 mm x 488 mm size. Each facet is diamond 
turned to a mirror finish yielding a reflectivity of > 85% in the visible band. The telescope 
with the F/D of ~ 1.2 will use graded focal length (25 m to 26.2 m) mirrors in order to 
reduce the Dgo spot size (defined as the diameter of the circle within which 80% of the 
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reflected rays lie) of the light collector to ~ 15 mm for on-axis incidence. Each mirror panel 
will be equipped with motorized orientation controllers for aligning them to form a single 
parabolic light collector. 

The focal plane instrumentation will have a photomultiplier tube based imaging cam- 
era covering a field of view of 4°x4°. The imaging camera will comprise of 1088 pixels 
arranged in a square matrix with uniform pixel resolution of 0.125°. The inner 576 pixels 
with a field of view of 2.4° x 2.4° will be used for generating the event trigger. The PMTs 
will be provided with acrylic front- aluminized light cones for enhancing the light collection 
efficiency of the camera. The signal processing instrumentation will also be housed within 
the camera and the acquired data will be sent to the control room over the computer network 
for processing and archiving. Detailed Monte Carlo simulation studies have been carried out 
using CORSIKA HI 7711 code and the results suggest that using a pixel threshold of >4pe and 
a 4 nearest neighbour pixel trigger, gamma-ray energy threshold of ~30 GeV is achievable 
by the MACE telescope. When operational by 2013 the MACE telescope will have the dis- 
tinction of being the largest high altitude 7-ray telescope in the world. Operating at a flux 
sensitivity of a few milli Crab units the MACE will augment international efforts in the field 
of ground based 7-ray astronomy. 
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